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HAUPTAUFSÄTZE 


Beitragzur Mechanik desbildsamen Verhaltensvon Flußstahl'). 
Von K. Hohenemser und W. Prager, Göttingen. 


1. Übersicht. Die vorliegende Arbeit bringt zunächst in kurzen Zügen einen systematischen 
Aufbau der Mechanik des bildsamen Verhaltens der Metalle unter der Annahme, daß diese 
als quasiisotrop angesehen werden können. Es werden dann Versuchsergebnisse mitgeteilt, 
die zum Teil in Übereinstimmung mit den bisherigen Ansätzen der Plastizitätsmechanik 
stehen, zum Teil ihnen widersprechen. Die Versuchsergebnisse werden von einem neuen 
allgemeineren Ansatz, welcher die bisherigen in gewisser Weise als Sonderfälle enthält, 
wesentlich besser wiedergegeben, weisen aber im übrigen darauf hin, daß eine vollständige 
Beschreibung des plastischen Verhaltens von Stahl mit Hilfe der Mechanik isotroper Kontinua 
nicht möglich zu sein scheint. 


2. Die Fließgrenze.. Die beim Fließen der Metalle möglichen Spannungszustände unter- 
liegen einer gewissen Beschränkung, der Fließbedingung: Eine Funktion der Spannungen 
behält während des Fließens konstanten Wert. Die Frage nach der Form dieser Fließbedingung 
wurde für Stahl bereits mit praktisch hinreichender Genauigkeit gelöst. Für nicht zu große 
mittlere Normalspannung 
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ist die Fließbedingung von o unabhängig”). Man spaltet daher vom Spannungstensor 
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I), Der Notgemeinschaft der Deutsenen Wissenschaft, dureh deren Unterstützung dem erstgenannten Verfasser 


das Arbeiten im Göttinger Institut für angewandte Mechanik ermöglicht wurde, sei aueh an dieser Stelle gedankt. 
>, Vgl. z. B. F. Sehleiceher, Z. f. ang. Math. u. Mech. 6 (1926), S. 199. 
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bezeichnet. Für isotropes Material kommen dann als Fließbedingungen nur Funktionen der 
Invarianten des Deviators in Frage. Setzt man «den Betrag der größten Differenz der Haupt- 
werte von X, gleich einer Konstanten, so hat man die Schubspannungsfließbedingung, setzt 
man die quadratische Invariante des Spannungsdeviators gleich einer Konstanten, so hat man 
die Misessche Fließbedinzung ?) 


Be 


wo % die Fließgrenze bei reiner Schubbeanspruchung ist. Für den Eintritt des Fließens 
bei Stahl ist in guter Annäherung die Schubspannungsbedingung maßgebend’), während bei 
fortgeschrittenen plastischen Verformungen besser die Misessche Fließbedingung erfüllt ist °). 
Es hängt dies wohl damit zusammen, daß zu Beginn ein ungleichmäßiges, örtlich be- 
schränktes Fließen auftritt, und erst mit beginnender Verfestigung eine gleichmäßige plastische 
Verformung möglich ist®). Wie die Lodeschen Versuche zeigen, genügen die Spannungen 
auch im Verfestirungsgebiet angenähert der Gl. (1), man hat lediglich die Konstante der 
reehten Seite nach Maßgabe der jeweiligen Verfestigung zu verändern ’). 

Bei Stahl ist die Zunahme von I,” mit wachsender Verformung so gering, daß diese 
(‚röße als konstant angesehen werden kann, wenn sich die Formänderungen nicht über ein 
allzu eroßes Gebiet erstrecken. 


3. Statisch bestimmte Aufgaben der Plastizitätsmechanik. Statisch bestimmt sollen solche 
Aufeaben der Plastizitätsmechanik genannt werden, die gelöst werden können, ohne auf den 
Zusammenhang der Spannungen mit den Verformungen einzugehen. Als Randbedingungen 
kommen natürlich nur Spannungsrandbedingungen in Frage. Damit eine Aufgabe statisch 
bestimmt wird, müssen mindestens drei voneinander unabhängige Gleichungen zwischen den 
Komponenten des Spannungstensors I oder denen des Deviators 2, bestehen. Das Typische 
der Differentialeleichungen statisch bestimmter Aufgaben der Plastizitätsmechanik ist darın 
zu sehen, daß sie sich von einem begrenzten Teil des Randes aus innerhalb der von den 
Randendpunkten ausgehenden Charakteristiken in den Bereich hinein integrieren lassen, 
während bei dem elastischen Körper die Spannungen im Innern eines Bereiches immer von 
den Randspannungen an «allen Punkten des Randes abhängen. 

Beim Torsionsproblem für einen prismatischen Stab aus plastischem Material sind die 
Charakteristiken die Orthogonalen zu den Schubspannungslinien, d. h. die Fallinien der 
Spannungsfläche. Auf dem Rande des Querschnittes ist die Normalkomponente der Schub- 
spannung null, die Tangentialkomponente %k. Die Spannungsfläche läßt sich dann vom Rande 
aus als Böschungstläche aufbauen. Die Ermittlung der Spannungsverteilung in einem tor- 
dierten Stab, der mit dem größten tragbaren Drehmoment belastet, also überall plastisch 
eeworden ist, gehört somit zu den statisch bestimmten Aufgaben der Pla.tizitätsmechanik. 
(Ganz anders geartet ist die gleiche Aufgabe für einen Stab, bei dem nur in einzelnen Gebieten 
der Fließzustand erreicht wird, während in den übrigen das Material sich elastisch verhält. 
Von dem im Plastischen liegenden Teil der Querschnittsberandung aus läßt sich die Spannungs- 
läche auch hier als Böschungsfläche aufbauen, doch wird das dureh ein Randstück beherrschte 
(Yebiet begrenzt einerseits durch die von den Endpunkten dieses Stückes ausgehenden Fall- 
linien der Spannungsfläche, andererseits durch die Grenze gegen das noch elastisch eebliebene 
(sebiet. Diese Grenze selbst kann nur mit Hilfe von besonderen Annahmen ermittelt werden, 
z. B. der Forderung, daß auch die Tangentialkomponente der Spannungen beim Überschreiten 
der Grenze zwischen elastischem und plastischem Gebiet stetig verläuft. 

Statisch bestimmte Fälle können weiter auftreten bei ebenen Problemen, die dadurch 
zekennzeichnet sind, daß die eine Hauptachse des Spannungstensors eine vorgegebene unver- 
änderliehe Riehtung hat, «die wir weiterhin als z-Riehtung bezeichnen. Zur Festlegung des 


Spannungszustandes genügt dann außer dem Betrag der Hauptspannung in der z-Richtung 


die Angabe des „ebenen“ Tensors mit drei Komponenten für die Ebene der beiden anderen 
Hauptachsen. Setzt man die Hauptspannnng in der z-Riehtung gleich null, so erhält man 


3, Mit 2,2 (%,-%,) ist das skalare Produkt von &, mit sieh selbst gemeint, also die Summe der Quadrate aller 
Komponenten von £ , entsprechend der Vorschrift über skalare Multiplikation von Tensoren. 
1) Verl. z. B.UC.A.M. Smith, Proe. Institut. Mech. Engg. 1909, 4, S. 1237. 
W, Lode. VDI-Forsehnngsheft 305 (1928) und Z. f. Physik 36 (1026). S. d12. 
Vel. z.B. H. Hort, VDI-Forscehungsheft 41 (1907). 
Die Frage, in weleher Weise die Verfestigung von den vorangegangenen Deformationen abhängt, wird zur 


(ref rer Institnt Fir aı ve wandte \lechanik experimentell henrbeitet. 
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die Spannungsverteilung, wie sie angenähert in einer dünnen Scheibe entsteht, die nur durch 
Kräfte in ihrer Ebene beansprucht wird. Eine wichtigere Gruppe statisch bestimmter Auf: 
eaben der Plastizitätsmechanik ist unter denjenigen ebenen Problemen enthalten, bei welchen 
die Hauptspannung in der z-Riehtung gleich dem arithmetischen Mittel der beiden anderen 
vesetzt wird, bei denen also der Hauptwert des Spannungsdeviators I, in der z-Richtung 
verschwindet. Wir werden später beim Eingehen auf die plastischen Formänderungen sehen, 
daß ein solcher Spannungszustand angenähert bei einer ebenen plastischen Deformation auf- 
tritt. Aus den Gleichgewichtsbedingungen und der Fließbedingung 


(a. - 0 Hi =4M, 


die in diesem Falle auch die Konstanz der größten Schubspannung bedeutet, erhält man eine 
partielle hyperbolische Differentialgleichung zweiter Ordnung für r, deren Charakteristiken 
eegeben sind durch \ 
dy, nn dy, k+r 
d.r k+r’ dx k—ı 


und somit zwei orthogonale Scharen bilden. Man überzeugt sich leicht, daß die Charakte- 
ristiken an jeder Stelle die Hauptschubrichtungen haben und daher indentisch mit den 
Hencky-Prandtlschen Gleitlinien sind °). Gelöste statisch bestimmte Aufgaben dieser Gruppe 
sind die Ermittlung des Spannungszustandes in einem ganz plastisch gewordenen Rohr unter 
Innendruck und in einem ganz plastisch gewordenen Material zwischen zwei reibenden einander 
parallelen oder gegeneinander geneigten Platten®). Nicht statisch bestimmt sind natürlich 
solche Probleme wie das Eindringen eines starren Stempels in einen elastisch-plastischen 


Körper. 


4. Der Zusammenhang der plastischen Formänderungen mit den Spannungen. Die Be- 
handlung allgemeinerer Aufgaben der Plastizitätsmechanik erfordert das Eingehen auf den 
Zusammenhang zwischen Spannungen und Formänderungen. Bei der Bildung von Ansätzen 
macht man von der Tatsache Gebrauch, daß es einen großen Bereich von Deformations- 
zeschwindigkeiten gibt, bei denen die Erscheinungen der inneren Reibung und der Relaxation 
vernachlässigt werden können. Wie Versuche von H. Hannemann und R. Yamada) zeigen, 
beträgt die plastische Volumdehnung von Stahl nur einen Bruchteil der elastischen. Wir 
nehmen daher den Zusammenhang der mittleren Dehnung 
| 
.) 
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(Ex -r ey + &2) 
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mit der mittleren Normalspannung auch im plastischen Zustand als den gleichen an wie im 
elastischen Zustand, also 
m —2 
+ @; 
2m +-1)G 


wo G den Gleitmodul und m die Poissonsche Zahl bedeutet. Der gesuchte Zusammen- 
ten) 


hang zwischen Spannungstensor 2 und Formänderungtensor FE reduziert sich so auf einen 
Zusammenhang zwischen den entsprechenden Deviatoren I, und 
| | 
Ex € .) x 7) .) Yx: 
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Ehe wir spezielle Zusammenhänge zwischen 2, und E, diskutieren, wollen wir einige 
erundsätzliche Bemerkungen über Tensorgleichungen vorausschieken. Bezeichnen A, B,C... 


Tensoren, so bedeutet die Gleichung 
FA,BEC..)=®, 


das zwischen den für die gleichen Grund- und Koordinatenrichtungen genommenen Kompo- 
nenten von A,B,C... neun entsprechende Gleichungen bestehen. Der funktionale Zu- 


x, I... Prandtl, Z. f. ang. Math. u. Mech. 1 (1921). S. 15, und 3 (1923), 8. 401. H. Heneky, Z. f. ang. Math. 
1. Mech. 3 (1923), S. 241. 

», Vel. z.B. A. Nadai, Der bildsame Zustand der Werkstoffe, Berlin 1927, S. 1311. 

10, 4. Hannemann und R. Yamada, Arch. f. Eisenhüttenwesen 4 (19301), S. 399. 


| * 
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sammenhang f (A,B,C...) 0) muß invarlant sein gegenüber Transformationen der Grund- 
und Koordinatenriehtungen. Eın soleher Zusammenhang zwischen Tensoren ist z. B. gereben 
dureh die lineare Beziehung 


BAFLBFEIF GER: 5 ar ee as 
Die Größen €. €, €, . . . können Funktionen von skalaren Invarianten der Tensoren A,B, 
(!... sein, also z. B. der skalaren Produkte je zweier Tensoren. Trägt man in einem neun- 


dimensionalen Raum vom Ursprung eines Koordinatensystems aus den Vektor auf, dessen 
Komponenten gleich den neun Komponenten eines gegebenen Tensors sind, so kann man 
diesen Vektor unter gewissen Einschränkungen als Bild des Tensors auffassen. Im folgenden 
werden wir alle Tensoren auf ein und dasselbe feste Koordinatensystem beziehen und können 
dann von der obigen Darstellung Gebrauch machen. Der Linearform (2) entspricht im Bild- 
raum eine Addition von Bildvektoren. Bestehen zwischen den Komponenten eines Tensors n 
voneinander unabhäneize lineare Gleichungen, so liegt der Bildvektor in einem (9 — n) 
dimensionalen Bildraum. Diese Tatsache wird uns später von Nutzen sein. 

Wir denken uns jetzt «die Deviatoren E, und I, von einem Parameter f abhängig, der 
z. B. die Zeit darstellen kann. Im elastischen Zustand ist die Änderung des Spannungs- 
deviators proportional der Änderung des Deformationsdeviators und unabhängig von den je- 
weiligen Spannungen und Deformationen : 


ar. a2, 
dt dt 


T 


Wir wollen zunächst annehmen, daß im plastischen Zustand die Änderung des Spannungs- 
deviators nieht nur von der Änderung «des Deformationsdeviators, sondern auch noch von den 
augenblieklichen Werten von E, und 2, abhängt. Als einfachsten Ansatz wählen wir die 
lineare Beziehung 


d E, N #2, 


I |  % | .) 
cET+t%2,+ 6; + BA a Enz 3). 
| 0 9 0 | ; lt I N dt (+) 
Die Größen e, bis e, können noch Funktionen der skalaren Invarianten der vier Tensoren 
ge 1} und jy zein, jedoch in der Weise, daß die Beziehung voraussetzungsgemäß 
/ ( 


zeitunabhängie wird, d.h. unabhängig davon, mit welcher Geschwindigkeit die Deformationen 
vorgenommen werden. Haben die Größen e, bis e, konstante Werte, dann ist die geforderte 
Zeitunabhängiekeit nicht erfüllt. Wir weisen nochmals ausdrücklich auf diesen erund- 
sätzlicehen Unterschied ın den Ansätzen der Plastizitätsmechanik und denen für innere 
Reibung und Relaxation hin. da hier Unklarheiten zu bestehen scheinen ''). Außer der Gl. (3) 
ist noch die Fließbedineung zu erfüllen: 


A 1? (4) 

i . . . Bü \ . ii C Cr er 
überdies sind noch zwei weitere Gleichungen zur Festleeung der drei Größen , 
ne 

4 4 K} 


erforderlich. 
5. Der Saint-Venant-Misessche Körper. Wir setzen in Gl. (3) e,=e,=Ö und haben 


d BE, 
Gy -—& ), 


3 dt 
Durch skalare Multiplikation mit I, und Berücksichtigung von (4) erhält man 


6: R dt 


A?- | 
e. (dE,:2,) 
Es ergibt sich dıe zeitunabhängize Beziehung 


2.-4E) 2, 


E, A A 


Die Bildvektoren von dE, und 2, sınd also gleichzerichtet, und es ist 


(2, dE,) 
l 


veE,, 


Vel. z. B. B. Gutenberg und H. Schlechtweg, Phys. Z. 31 (1931), S. 745, die den „Fließvorgang* mit 


Hilfe der Gleichungen eines elastischen Körpers mit innerer Reibung und Relaxation beschreiben wollen. 
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wir haben daher 
dE +2 


’0 


) d FE, “ | 


(N). 


Dem Bildvektor d FE, einer infinitesimalen Deformationsänderung entspricht eindeutig 
ein Spannungsbildvektor 2, von gleicher Richtung wie dE, und mit dem absoluten Be- 
trage A. Bei Vorgabe von &, bleibt jedoch der Betrag von dE, unbestimmt. Vernachlässigt 
man die Volumdehnung. was bei genügend großen plastischen Formänderungen stets zulässig 
ıst, dann sind die Gl. (5) zusammen mit der Bedingung 


E 0 


identisch mit den Misesschen Gleiechungen'?). Die hier gegebene Formulierung erscheint 
uns glücklicher als die übliche, da sie die Zeitunabhängiekeit dieser Gleichungen klarer her- 
vortreten läßt. Ändert sieh , während der ganzen Deformation nicht, so kann man 
(il. (9) auch schreiben: 


(6). 


Die Gültiekeit dieser Gleiehung ist dureh Versuche von W. Lode'*) erwiesen worden und 
wird auch durch unsere später mitzuteilenden Versuche bestätigt. Man kann also bei un- 
verändertem Spannungsdeviator I, nur solche Deformationen E, vornehmen, bei denen der 
Bildvektor von FE, gleiche Riehtung wie der von 2,, im übrigen aber einen beliebigen ab- 
soluten Betrag hat. Wir wollen derartige Formänderungen „frei“ nennen, da sie ohne Ände- 
rung der äußeren Kräfte vorgenommen werden können. Die zu einem bestimmten Spannungs- 
zustand gehörende freie Deformation wird z. B. dann auftreten, wenn man die Spannungen 
in beliebiger Weise «die Fließerenze um einen unendlich kleinen Betrag überschreiten läßt. 
Der Saint-Venant-Misessche Ansatz gibt also die einem bestimmten Spannungszustand 
entsprechenden freien Formänderungen in Stahl richtig wieder. Derartige Deformationen 
können auftreten bei den statisch bestimmten Problemen der Plastizitätsmechanik im Sinne 
der Definition in 3. Nachdem die Spannungen im Innern aus den Randspannungen berechnet 
worden sind, ist die Formänderungsaufgabe, wie aus Gl. (6) hervorgeht, noch in keiner Weise 
eindeutie lösbar, vielmehr ergeben sich für die Formänderungen partielle Differential- 
eleichunzen, zu deren Lösung noch die Angabe der Randwerte der Formänderungen erforder- 
lich ist. Bei denjenigen ebenen Problemen, bei welehen die Hauptspannung in der z-Riehtung 
eleich dem arıthmetischen Mittel der beiden anderen ist, folgt aus Gl. (6) 


€, e=U). . 


oder, da «die Volumdehnung & gegenüber den plastischen Dehnungen meist vernachlässigt 
werden kann, &.—=0, d. h. die Verformung ist eben. Die Charakteristiken der Differential- 
oleiehuneen für die Deformationen sind in diesem Falle identisch mit den in 3 erwähnten 
Prandtl-Henekyschen Gleitlinien '*), wodurch die Bezeichnung der Hauptschubspannungs- 
linien als Gleitlinien erst gerechtfertigt wird, da Unstetigkeiten der Verschiebungsrandwerte sich 
längs dieser Linien in das Innere des Körpers fortpflanzen. Die Randwertaufgaben der 
Plastizitätstheorie, die hinsichtlich der Spannungen und Verformungen vollständig gelöst sind, 
sind die unendlich ausgebreitete Masse zwischen zwei parallelen reibenden Druckplatten und 
das Rohr unter Innendruck '’). 


Wir diskutieren weiter das Verhalten des Saint-Venant-Misesschen Körpers bei 
beliebigen Formänderungen. Aus Gl. (5) folgt, daß bei einer unstetigen Riehtungsänderung 
von dE, auch I, die gleiche unstetige Richtungsänderung mitmacht. Dies würde z. B. be- 
leuten, daß in einem dünnwandigen Rohr, welches zuerst plastisch verdreht und dann gedehnt 
wird, die Torsionsspannungen im Augenblick des Deformationswechsels verschwinden und 
dafür die volle Fließzugspannung entsteht. Damit hängt zusammen, daß in dem Saint- 
Venant-Misesschen Körper, an welchem gerade keine Deformationsänderungen vor- 
genommen werden (4 F, =), entsprechend Gl. (5) ein beliebiger Spannungszustand an der 
Fließgerenze herrschen kann, während bei den Metallen auch ın diesem Falle ein dureh die 
Deformationen eindeutig bestimmter Spannungszustand auftritt. Infolge Niehtberücksichtigung 
des elastischen Anteils der Formänderung gibt der Saint-Venant-Misessche Ansatz das 
wirkliche Verhalten der Metalle in diesem Falle nieht richtig wieder. Tatsächlich haben auch 

12, R.v. Mises, Göttinger Nachrichten math. phys. Kl. 1913, S. >82. 

13) W. Lode, vg 


14) 1. v. Mises, Z. f. ang. Math. u. Mech. 5 (1925), 8. 147. 
15, 11. Pollaezek-Geiringer, Verhdle. III. int. K. f. tech. Mech, Stockholm 1950, II, S. 185, 











Versuche des erstgenannten Verfassers '") gezeigt, daß in dem obigen Falle des erst plastisch 
verdrehten und dann gedehnten Rohres aus Stahl die Torsionsspannung mit wachsender 
Dehnung lanzsam und stetige absinkt und noch bei Dehnungen, die ein Vielfaches der elastischen 
Dehnung an der Fließgrenze ausmachen, einen von Null verschiedenen Betrag hat. Im gleichen 
Maße und unter recht guter Einhaltung der Fließbedingung (4) steigt die Zugspannung all- 
mählich an. Bei nicht statisch bestimmten Aufgaben (z. B. gemischten Randwertaufeaben) 
dürfen die Saint-Venant-Misesschen Gleichungen daher nicht ohne weiteres anzewandt 
werden, da sich dann ım allgemeinen die Richtung von 3, während der Deformation 
ändern wird. 


6. Der Prandtl-Reußsche Körper. Da der Saıint-Venant-Misessche Ansatz die freien 
"ormänderungen richtig wiedergibt, die elastischen dagegen vernachlässiet, kann man daran 
denken, jede infinitesimale Formänderung dE, aus einer freien Deformationskomponente 


2.62) 2, 


- u | 
und einer elastischen Komponente 


Ad 


UE () 


dE, 


‘0 = I(7 
zusammenzusetzen. Wir haben dann 


IE (2, ° d F,) 3 “2, (Y 
A "A786 1 
Dazu kommt noch die Beziehung zwischen der mittleren Normalspannung o und der Volum- 


dehnune 


u 4 


wor, 
I(m-+-1)G 


Die Beziehung (9) folgt auch aus Gl. (3), wenn man e, identisch null setzt und die Fließ- 
bedingung (4) berücksichtigt. Es ergibt sıch 

ö- (2,:dE,) 
C. A’dt 


r 
\us der Forderung, daß auch ' nicht von E, abhängen soll, läßt sich weiter folgern 
96 


wenn man nämlich eine Formänderung länes der Fläche E, = ım #,-Bildraum vornimmt. 


2G 
Die Gl. (9 ist wieder eine, wie uns scheint, glücklichere Formulierung der Reußschen 
Gleichungen '”), da sie die Zeitunabhängigkeit der Be- 
ziehung klarer hervortreten läßt'‘). In Gl. (9) ist 


oe 
(2,48. 2, 





‘0 
l 1 

die Komponente von dE, in der hRiehtung von 2,, 
d2, steht wegen I,’ = 4? oder (2, : d&,)=0 senkrecht 
auf &,. IstdE, und 2, gegeben, so erhält man nach 

= 
Gl.) ,,, dadurch, daß man die Komponente von 

ud LT 


dE, in Riehtungz senkreeht zu 2, nimmt (vgl. auch 


Abb. 1). Ist umgekehrt d 2, zegeben, dann erhält man 
von dE, nur die Komponente in Richtung von d2,, 
während die Komponente ın Richtung von 2, unbestimmt bleibt. Selbstverständlich gilt die 
(1. (Od) nur dann, wenn die Komponente von d E, in der Richtung von 5. positiv ist, da andern- 


falls der Körper nieht mehr im plastischen Zustand verbleibt. 





\bhh. l 


16), X. Hohenemser, Z. f. ane. Math. u. Mech. 11 (1931). S. 1». 


17), X. Reuss, Z. f. ang. Math. u. Mech. 10 (1930), S. 267. 
l,.. Prandt! teilte den wleiehen Ansatz wie A. Reuß Ffir den Sonderfall des ebenen Formänderungszu 
standes. allerdines mit anderer Fließbedineune, schon 1924 in Delft mit, Verhdlg. des I. int. K. f. tech. Mech., 


Delft 1024, 8. 45. 
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In Abb. I sei AB ein Formänderungeswee im E,-Bildraum. An jeden Punkt des Form- 
u 
änderungsweges soll der Vektor der elastischen Formänderung IT anzetraegen werden, der 
wo) UT 


die konstante Länge IT: hat. Die Zerlegung eines jeden Elementes des Formänderungsweges 
u ' 


in einen freien Teil 
(2, dE,) 3 


| A 
und einen elastischen Teil 
AS, 
2G 
d2 & 
ist ın Abb. I durchgeführt, und die YA sind zu den jeweiligen 50 addiert worden. Man 
. \T . oT 
sieht, daß im Grenzfall der Endpunkt 0 von 56 auf einer Schleppkurve liegt zum Form- 


änderungsweg, den der Endpunkt des Vektors E, beschreibt. Die Schleppkurve zu einer 
Kurve wird dabei dadurch gekennzeichnet, daß je zwei entsprechende Punkte von Kurve und 
Schleppkurve einen konstanten Abstand voneinander haben, und daß ihre Verbindungslinie 
Tangente an die Schleppkurve ist. Die Schleppkurve zu einer gegebenen Kurve wird be- 
stimmt durch Angabe eines Paares entsprechender Punkte, 

Man erkennt den Fortschritt der Beziehung (9) gegenüber der Gl. (5), denn einer un- 
stetigen Riehtungsänderung von dE, entspricht nunmehr eine stetige Änderung der Riehtung 
von 2. Für den Fall, daß die Bildvektoren der Gl. (9) alle in einem zweidimensionalen Raum 
liegen, ist die Integration der Reußschen Gleichungen mechanisch auf Grund der Schlepp- 
kurvenbeziehung mit Hilfe des Beilplanimeters möglich. Man erhält so die Spannungen am 
Ende eines vorgezrebenen Formänderungsweges eindeutig, wenn dieselben am Anfang des 
Weges gegeben sind. Die umgekehrte Aufgabe, aus dem Verlauf der Spannungen denjenigen 
der Formänderungen zu bestimmen, ist nieht eindeutig lösbar, entsprechend dem tatsächlichen 
Verhalten der plastischen Metalle. Die kennzeichnende Eigenschaft des Prandtl-keußschen 
Körpers, nämlich die Unabhängt!gkeit der Spannungs-Formänderungs-Beziehungen von der 
Anfangsdeformation FE, äußert sich in unserer Darstellung der Abb. I «darın, daß man die 
Figur 9ABb im E,-Bildraum beliebig verschieben darf, sie gibt immer die Spannungsänderung 
auf dem entsprechenden Wege wieder. 

Die mechanische Integration der Reußschen Gleiehungen gelingt bei denjenigen ebenen 
Problemen, bei welchen die Hauptspannung für die z-Riehtung den Wert 

| 


OÖ. Ö .) 6x 0 


Mr 


annimmt, vgl. auch 3. Die Gl. (9) heißt für die e-Komponente 


(2, - d FE) d(o. 6) 
die. e) FE (6. 6) >; 
und man erhält wegen o- 6--0 hieraus 
€, E GONSE. 
Weil o, Hauptspannung ist, gilt weiter 7... 7,- 0 und es folgt aus den Gl. (9 
"xz  eonst und  y,- = const. 


Da zu Beginn des Fließens Formänderungen vorhanden sein sollen, wie sie nach dem Elasti- 
zitätsgesetz den Spannungen entsprechen, ist auch während des Fließens &. 2 ey: =yy- =. 
Aus der Kleinheit der Volumdehnung & folgt, daß angzenähert ebene Deformation vorliegt. Der 
Spannungs- und Formänderungszustand ist folglich bestimmt durch die ebenen Tensoren 














| 
Öx T Ei =, 
\, . 
> und # 
| 
T ou ») Y ey 
Die entsprechenden ebenen Deviatoren heißen F 
Ox Oy x Ey Y 
.) T .) -) 
vr A 
ar und E, 
Öx Oy y { £ 7} 
T .) .) .) 
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Zwischen je zwei entsprechenden Komponenten der beiden Deviatoren besteht die gleiche 
invariante Beziehung, die Bildvektoren 2, und 2G E, liegen daher in einer zweidimensionalen 


Ebene des vierdimensionalen Bildraums und haben die Komponenten 


= 5”, ryY2 und 20: 5’, 20 75. 
Die Fließbedinzung wird hier 
we (6, m 37: 3, 
Die Gleichung (9) gilt also auch für die ebenen Deviatoren I, und FE, und wir können die 


in Abb. I angedeutete Konstruktion mit Hilfe des Beilplanimeters ausführen'®). Ist der 
Verlauf der Deformationen gegeben, dann findet man mit Hilfe der Schleppkurvenkonstruktion 
die Spannungsgrößen o, o, und r. Außerdem bestehen noch die Gleichungen 


.) x 
m — 2 0x +0, 
F £ Bu : 
In +-1)G > 
j Öx jr. Öy r PR j (0) 
0; > ‚ Tı,=T'Ty; ET Du 7 de 


so daß die Spannungen eindeutig bestimmt sind. Setzt man den Körper als inkompressibel 
voraus, d.h. m =2, e=0, so ist die ebene Deformation streng realisiert, dagegen bleibt jetzt 
der Wert von (9, + 9,)} unbestimmt. 

Bei demjenigen ebenen Zustand, für welchen die Hauptspannung in der z-Richtung 
verschwindet und außerdem für eine dazu senkrechte Richtung die Normaispannung null 
ıst, also 
0> 0, 


OÖ, 


läßt sieh die mechanische Integration der Reußschen Gleichungen gleichfalls durchführen. 
Aus Gl. (9) folet wie oben, daß während der Formänderung 


Yz, = const, Yy, > const 


4 


ist. Da zu Beginn der Bewegung der Körper wiederum Deformationen haben soll, wie sie 


nach dem Elastizitätsgesetz den Spannungen entsprechen, folgt aus 7.,-=rT,z;=0 auch y.-= 
„-=V. Aus 0,=6.=0 ergibt sich ©, e., und die Deviatoren I, und E, heißen: 
-) .) nr 
-0, I 0 = W, Ey) = 0 
) +) i _ 
| | 
a er 
N } — 6, v und E, = > (£,—£,) 0 
+.) dd +) 
(0) 0) — (. 0 0 3 (Ex €) 
+.) e 














Die Bildvektoren von I, und 2 GE, liegen für beliebige Spannungen und Deformationen in 


ein und derselben Ebene des Bildraums und haben ın dieser Ebene die Komponenten 


.) ») ar 
Öx >, y2 und 2G (er, —e,) u, UT: 
+) ö +.) ] _ 
Bei «(dem zuletzt besprochenen ebenen Zustand o,= 0.==0 ist auch die Form 6°-+ e? 7? 


mit beliebigem e eine Invariante des Spannungsdeviators 2. sie kann daher als Fließbedingung 
verwandt werden. Bestimmt man ın den Reußschen Gleichungen 


’ 4 x f "EM 
d E, Aa 3, 
die Größe 7 mit Hilfe der Fließbedingung 
rettete, 


dann erhält man Differentialgleichungen, die, wie man sich leicht überzeugt, mit Hilfe des 
Beilplanimeters integriert werden können, wenn man in der Bildebene 


G,, er und 2Gle,—e,), Gcy 


‘ 


aufträgt. 


Die Darstellung eines ebenen Deviators dureh einen ebenen Vektor und die Angabe der Schleppkurvenkon 
struktion ım obigen Falle stammt von L. Prandtl, siehe a. a. ©. 
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7. Versuchsergebnisse. Die Versuche wurden an dünnwandigen Stahlrohren von 34 mm 
mittlerem Durchmesser und 3 mm Wandstärke angestellt, die in einer Spezialmaschine 
eleichzeitig gedehnt und verdreht werden konnten’). Als Fließbedingung konnte mit guter 
Annäherung 


u a ee a 
verwendet werden. wobei e etwa zwischen 1.70 und 155 schwankte. kKın Teil dieser 


Schwankungen rührt von der Verfestigung während des Versuches her, die sieh im 0°  r’ 
Diagramm (Abb. 2) dadurch bemerkbar macht, daß die Fließkurve bei während des Versuchs 
fallender Schubspannung flacher, bei steigender Schubspannung 
steiler verläuft. Doch auch wenn man den Einfluß der Ver- 
festigung durch aufeinanderfolgende Versuche mit fallender und 
steigender Schubspannung zu eliminieren versucht, bleiben 
Schwankungen der e-Werte von etwa + 3vH. Da diese Schwan- | 
kungen in systematischer Weise von den vorausgegangenen Form- | er 
änderungen des Rohres abhängen, liegt der Schluß nahe, dab Mm 0, 
die bei großen plastischen Formänderungen nachweisbaren Abb. 2. 
Anıisotropieeigenschaften (Zieh- und Walztextur) auch schon 

bei den kleinen plastischen Verformungen unserer Versuche Einfluß haben (vgl. auch 
Anmerkung 7). 


. 


Wir können die Versuchsergebnisse in der in 6 beschriebenen Weise darstellen und 
damit außer einer guten Übersicht über den Verlauf der Spannungen und Formänderungen 
auch eine Prüfung der Versuchsergebnisse mit dem Beilplanimeter ermöglichen. Der Form- 
änderungsweg ist in diesem Fall also in einem Koordinatensystem darzustellen, in welehenn 


2G(e,—e,) und Gecy 


aufzutragen sind. Zu jedem Formänderungspunkt des Weges gehört der Spannungsbildvektor 
mit den Komponenten o, under. Die Länge dieses Vektors ist entsprechend der angenäherten 
Gültiekeit von Gl. (10) angzenähert konstant. 


Um den Vergleich von Versuchen zu ermöglichen, die in verschiedenen Verfestigungs- 
stadien vorgenommen sind, machen wir die aufzutragenden- Größen dimensionslos, indem wir 
sie alle durch die Fließspannung bei reinem Zug 077, dividieren. Unter Benutzung der 
Beziehung 


ETyy Of] 
und der Definitionen 


| 71 . Tyı 
ei — E z / A (} 


4 


wo E und @G Elastızitäts- und Sechubmodul sind. erhält man: 


an Ey 2(r Ex € lıcy y OÖ, 1 T 
au — _. =e, y und 8, / 
of l- e/ı of ‚I of or Try 


Berücksichtigt man noch den Zusammenhang der mittleren Normalspannung mit der Volum- 
dehnung: 


im y 


i > Di 
&xt&ry Ox, 


3m -+-1)G 


dann erhält man für e die Beziehung 


B) m €, m 2 
== S. 
2 (ml) ey, Im—+|]) 


Die Dehnungen und Verdrehungen konnten mit Spiegelgeräten und Fernrohr gemessen werden, 
dla die während eines Versuches auftretenden Formänderungen klein waren. Die Auswertung 
der Versuche erfolgte folgendermaßen: Die Zugekräfte P und Drehmomente M wurden in 
a M\? 

einem Diagramm mit den Koordinaten 7° und | aufgetragen, wo r „, der mittlere Rohr- 


l in 


radıus ist. Die Punkte liegen dann sehr angenähert auf einer Geraden, vgl. Abb. 7 und 8, 


20, Eine Beschreibung der von der Helmholtz-Gesellschaft zur Förderung physikalisch-teehnischer Forschung 
zur Verfügung gestellten Maschine findet sieh bei K. Hohenemser, siehe a. a. ©, 
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welehe zu den Versuchen der Abb. 5 gehören. Die Gerade schneidet auf den Achsen die 


Myr\ 
Größen PP’; und | |) ab. Man erhält dann 


‚m 


O, pP r MW 
r 07 Pyı' sy  Myr 


Aus einem Versuch im elastischen Bereich, für welchen die Beziehungen 
I’ ce ex b) \/ 9) 


velten, wobei #, und y Dehnung und Schiebung in Skalenteilen der Meßgeräte sind, erhält 
man ec, und e, und daraus 


P’;, My, 
fl . m 
e- e,. 
10 
une weiter mit _ 
+.) 
5 
e 2... —=s und g9=- 
13 ef] 13 F vs 


Die Verfestigung während der Versuche war gering, im Maximum 4 vH. Eine Berücksichtigung 
der Verfestigung wurde bei der Auswertung einiger Versuche durebgeführt, ergab jedoch nur 
sehr eerinefüriee Korrekturen. 


> 


In Abb. 3 ıst ein Fließversuch bei unver- 
änderten Spannungen dargestellt, um noch einmal 
die Gültigkeit der Gl. (6) zu erhärten. Die Span- 
nungsbildvektoren sind der Übersichtliehkeit halber 
um 90" gedreht gezeichnet. In Abb. 4 sind zwei 
Versuche eingetragen, bei denen nach einem ver- 
schieden langen Stück Fließen unter unveränder- 
lichen Spannungen die Dehnung konstant gehalten 
und das Rohr lediglich verdreht wurde. Die Ver- 
suche lassen sich dadurch zur Deckung bringen, 
daß man die Formänderungskurven zusammen mit 
den Spannungsvektoren parallel von A nach B ver- 
schiebt. d. h. ein Fließen unter unveränderlichen 
Spannungen hat keinen Einfluß auf den weiteren 
Verlauf des Zusammenhangs zwischen Spannungen 
und  Verformungen. Die Grundannahme des 
Prandtl-Reußschen Ansatzes, Unabhängigkeit 
des Spannungs-Verformungs-Zusammenhangs von der vorausgegangenen Verformung, ist 
also in diesem Falle erfüllt.  @Quantitativ allerdings bestehen beträchtliche Abweichungen 
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\bb. 4. 
der Versuchsergebnisse von dem Prandtl-Reußschen Ansatz, wie man an den ge- "= 
strichelten Spannungsvektoren sieht, die mit Hilfe des Beilplanimeters bestimmt wurden *'). 2 


\rl. auch die Versuche von K. Hohenemser, siehe a. a. Ö. 
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Ks erhebt sieh nun die Frage, ob nicht der Prandtl-Reußsche Ansatz, trotz 
quantitativer Abweichungen, das bildsame Verhalten von Stahl im Prinzip richtig wiedergibt, 
insbesondere ob der Spannungs-Formänderungs-Zusammenhang wirklich völlig unabhängig 
ist von den vorauszerangeenen Formänderungeen. Hierüber geben die zwei in Abb. 5 dar- 
eestellten Versuche Aufschluß. Die Formänderuneswege sind in beiden Fällen bis auf den 
kleinen „Haken“ bei A identisch. Der Haken im Formänderungsweg bedeutet, daß ein kurzes 
"ließen unter unveränderten Spannungen stattgefunden hat. Man erkennt die einschneidende 
Wirkung, die diese kleine freie Formänderung auf den weiteren Verlauf der Spannungen aus- 
übt. Bis zum Punkte A stimmen beide Versuche überein, von dort ab dagegen zeigt sich 
eine wachsende Abweichung. Die striehpunktierten Vektoren in Abb. 5 oben geben die 
Spannungen des Versuchs Abb. 5 unten wieder. Abb. 6 zeigt einen Versuch, bei welchem 
nach anfänglichem Fließen unter konstanten Spannungen, wobei der Spannungsbildvektor 
die gleiche Riehtung hat wie derjenige in den Versuchen der Abb. 5 im Punkte 4, im 
weiteren der eleieche Formänderungsweez wie in Abb 5 eingeschlagen wurde. Der Verlauf 
der Spannungen bei diesem Versuch ıst der gleiche wie der in Abb. 5 oben vom Punkte 4A 
an. Man erkennt aus diesen sehr bemerkenswerten Versuchen, daß der Prandt]l-Reußsche 
Ansatz nicht nur quantitativ von den Versuchsergebnissen abweicht, sondern die Fließvorgänge 
in Stahl erundsätzlich nicht richtig beschreibt. 


8. Ein neuer Ansatz. Aus den in 7 mitgeteilten Versuchen entnehmen wir, daß eine 
beliebiee Deformationsänderung sıch nicht einfach aus einer elastischen und einer freien 
Formänderungskomponente zusammensetzen läßt. Der Fließmechanismus bei beliebigen 
plastischen Formänederuneen ıst offenbar wesentlich komplizierter, trotzdem fällt den freien 
Deformationen eine ausgezeichnete Rolle zu. Aus Abb. 3 folgt, daß sie auf alle späteren 
Voreänee im Stab keinen Einfluß haben. aus Abb. 5 und 6 folet, daß bereits eine kleine 
(reie Formänderung genügt, um gewissermaßen das Gedächtnis des Materials an alle voraus- 
eerangenen Vorgänge auszulöschen. Denn der Hasen im Formänderungsweg bewirkt, daß 
das Material sich im weiteren so verhält, als hätte die Formänderunge von B nach A, die in 
dem Versuch der Abb. 5 unten einen wesentlichen Einfluß auf den weiteren Verlauf der 
Spannungen hat, nicht stattgefunden. Wollen wir diese Tatsache durch einen Ansatz der 
Kontinuumsmeehanik wiedergeben, so muß offenbar ın diesem Ansatz die seit der letzten 
freien Formänderung stattgefundene Deformation eingehen. Im Bildraum des Formänderungs- 
deviators (Abb. 9) sei I1BC ein Formänderuneswee und B der Punkt. an dem zum letzten 


Male von € aus gesehen eine freie Formänderung stattzefunden, d.h. an welchem zum letzten 
kn 


Male der Bildvektor . ., die Riehtung der Tangente an den Formänderuneswege eehabt hat. 
-) (1 te es 


x 


Wir bezeiehnen mit E," die Deformationsänderung zwischen B und €, mit E, : 5 die elastische 
' 


"ormänderune ın DB. mit 
A E+-HK" 


die Summe dieser beiden Formänderungen. Wir knüpfen nun an Gl. (3) an, in die wir statt 
F, die Größe E,* einsetzen. Im übrigen fordern wir e,=c,=0, wie es die Versuche nahe- 


A 
leeen. Mit Berücksiehtieunge der Fließbedineunge (4) erhalten wir °°) 
F, 2 11 
: N". rn te ea ee A en ; 
VE, .| | 
Außerdem gilt noch 
u 24 12 
CE er eh Wa Wie ee er ES 
>(m-+1)G 


Der neue Körper, der dureh die Gl. (1b) und (12) beschrieben wird, hat zwei kenn- 
zeiehnende Eigenschaften. Die erste ist darin zu sehen, daß ım Gegensatz zu dem Prandtl- 
Reußschen Körper verschiedene Formänderungswege nach dem gleichen Formänderungs- 
zustand die gleichen Spannungen ergeben, wenn nur freie Formänderungen auf allen Ver- 
eleiehswegen vermieden werden. Im E,-Bildraum gibt es also ein Feld für die Spannungs- 


>), Die von H. Henekv. Z. f. angew. Math. u. Mech. 4 (1924), S. 325 aufgestellten Gleichungen stimmen formal 
mit den Gl. (1b) und (12) überein, allerdings mit dem wesentlichen Untersehied, daß Heneky die Gesamtdeformation 
einsetzt statt E,* und seine Gleiehnngen deshalb im Widerspruch zu den Beobachtungen an Metallen im plastischen 
Zustand stehen. MH. Henekv leitet seine Gleiehungen aus einem Variationsprinzip her, dessen Gültigkeit schon 
\. Haar und Th. v. Karman,. Gött. Naehr. math. phyvs. Kl. 1909, S. 204 postuliert haben, doch ist die Herleitung auch 
dann nieht zwangesläufie, wenn man die wahrscheinlich zutreffende Gültigkeit des verwendeten Spannungsvarliations 
prinzips voranssetzt, enthält vielmehr die willkürlicehe Annahme der Identität der in dem Variationspreblem auf 
tretenden Laerangesehen Fakteren zu den Gleiehgew ichtsnebenbedingungen mit den Verschiebungsgrößen. 
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bildvektoren ra diese haben die Riehtung der Vektoren vom Ursprung E,"=0 nach den 
nd 1 

betreffenden Formänderungspunkten und den dureh die Fließbedingung vorgeschriebenen 


Betrag Man nennt solehe Verformungen vielleicht zweckmäßig quasielastisch, da ja 


2@ 
hier wie im elastischen Gebiet auch jedem Formänderungszustand eindeutig ein Spannungs- 
zustand zugeordnet ist. 

Die zweite kennzeichnende Eigenschaft des neuen Körpers ist darin zu sehen, daß der 
quasielastische Zustand labil ist. Eine beliebig kleine Störung im Formänderunesweg, die mit 
einer kurzen freien Formänderung verbunden ist, bewirkt eine Verlagerung des Nullpunktes 
für E,*. UÜberlagert man dem Fließvorgang Schwingungen, die im Verhältnis zur Fließ- 
eeschwindiekeit genügend rasch sind und beliebig kleine Amplitude haben, so wird der labile 
quasielastische Zustand dauernd gestört, und der neue Körper verhält sich in diesem Falle, 
wie man leieht einsieht, wie der Prandtl-Reußsche Körper. Der neue Ansatz enthält also 
als Spezialfälle sowohl den Saint-Venant-Misesschen Körper, wenn nämlich die Form- 
änderung bei unveränderten Spannungen erfolgt, wie auch den Prandtl-Reußschen Körper, 
wenn der labile quasielastische Zustand dauernd gestört wird. 


9. Weitere Versuchsergebnisse. Die Versuche in 7 werden alle von dem neuen Ansatz 
gut wiedergegeben, d. h. die Spannungsbildvektoren gehen recht genau durch den Punkt 


EX=0. Die gestrichelten Vektoren stellen immer die Spannungen nach Prandtl-Reuß 
dar. ‚Jeder Versuch wurde mit einer freien De- 


) : N 
[formation eingeleitet, und zwar aus folgendem | 8 


Grunde: Zwischen dem eigentlichen Fließgebiet 
und «dem elastischen Zustand liegt ein Über- 
eangsgebiet (AB in Abb. 10), in diesem Über- 
gangsgebiet herrscht nicht nur eine andere Ge- 
setzmäßigkeit als im Fließgebiet, sondern der 
Verlauf der Spannungen auch jenseits von B 
hängt davon ab, in welcher Weise man 
Übergangsgebiet (lurehsehritten hat. Um einen 
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einheitlichen Ausgangszustand für die Versuche [RArzzm) Verformung 
zu erhalten, erwies es sich daher als not- Abb. 10. 
ge 
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Abb. 12. 


wendig, nach Durchschreiten des Übergangszustandes freie Formänderuneen vorzunehmen. 
durch welche der Einfluß der vorausgeganzenen Voreänge auszelöscht wird. Obwohl Versuche, 
die hier nicht wiedergegeben werden sollen, darauf hinweisen, daß auch eine Entlastung in 
das elastische Gebiet hinein die gleiche Wirkung hat wie eine freie Deformation, und die Er- 
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innerung an die vorangegangenen Vorgänge auslöscht, verhindern die undurchsichtigen Vor- 
vänge im Übergangsgebiet, diese Tatsache mit der wünschenswerten Exaktheit zu zeigen. Wir 
haben daher davon abgesehen, im neuen Ansatz auch die Labilität des quasielastischen Zu- 
standes gegenüber Entlastungen unter die Fließerenze zu fordern, da diese Erscheinung ver- 
suchsmäßig nicht einwandfrei realisiert werden konnte. Überhaupt sei ausdrücklich bemerkt, 
daß der neue Ansatz zwar die plastischen Eigenschaften von Stahl in großen Zügen richtig 
wiedergibt. aber gezenüber dem tatsächlichen Verhalten eine starke Idealisierung darstellt. 

Dies zeigt sich z. B. in dem Umstand, daß bei 








: , allen Versuchen ım fortgeschrittenen Stadium, 

hs : wenn der Winkel zwischen der Tangente an 

4 den Formänderungsweg und dem Spannungs- 

ce 4 , n vektor klein geworden ist, die Spannungen sich 
u wi f so verhalten, als ob freie Deformationen statt- 
| E e gefunden hätten, obwohl keine meßbaren Störun- 
. “ gen im Formänderungsweg aufgetreten sind. 

% / . In Abb. 11 ist der Versuch von Abb. 6 vollständig 
wiedergegeben. Man sieht, daß im späteren 

| UP, Verlauf Abweichungen vom  quasielastischen 
lo Verhalten in Richtung des Prandtl-Reuß- 
5 schen Verhaltens auftreten. Die strichpunk- 

ge" tierten Vektoren stellen wieder die Spannungen 

bei ungestörtem «uasielastischen Verhalten dar, 

e die gestrichelten diejenigen beim Prandtl]- 


Rkeußsehen Verhalten. Daß es sich hierbei 


Eu offenbar um eine typische Labilitätserscheinung 








— 05 109 handelt und nicht um eine gesetzmäßige Ab- 
mm ee weichung von dem quasielastischen Verhalten, 

" eeht aus dem bisher vorliegenden sehr umfang- 

\bh. 15 © 


reichen Versuchsmaterial hervor. Versuche mit 
eleichem Formänderungsweg stimmen zunächst gut miteinander und mit den quasielastischen 
Gleichungen überein, in einem späteren Stadium weichen sie voneinander ab, und zwar in dem 
Sinne, daß ein Mitschleppen des E,"-Bezugspunktes bei manchen Versuchen früher, beı 
manchen später einsetzt. 


Kin Fall des mehrfach gestörten quasielastischen Zustandes ist in Abb. 12 wiedergegeben. 
Man sieht. daß hier das Prandtl-Reußsche Verhalten besser erreicht ist. In einem ge- 
nügend kleinen Formänderungsbereich ıst die Feldeigenschaft der Spannungsvektoren gut er- 
füllt, wie Abb. 13 zeigt, in der verschiedene Versuche übereinandergetragen sind, so daß die 


FE -Bezugspunkte zusammenfallen. Die gestrichelten Vektoren in den Schnittpunkten zweier 


Formänderungswege stellen wieder die zu den betreffenden Wegen gehörenden Spannungen 
nach dem Prandtl-Reußschen Ansatz dar. Im späteren Stadium all dieser Versuche er- 
eben sich Abweichungen vom quasielastischen Verbalten in Richtung nach dem Prandtl- 
keußschen Verhalten hin. Wenn auch die Frage nach einem Kriterium für diese Ab- 
weichungen, nach einem Maß der Labilität, noch durchaus ungeklärt ist, kann man doch als 
durch die Versuche feststehende Tatsache ansehen, daß die quasielastischen und diePrandtl- 
Reußschen Gleichungen «die äußersten Grenzen darstellen, zwischen denen das wirkliche 
Verhalten von Stahl ın allen Fällen liegt. 


10. Zusammenfassung. Wir fassen die Ergebnisse unserer Untersuchungen noch einmal 
zusammen. Die Misessche Fließbedingung ist bei veränderlichen kombinierten Zug- und 
Schubbeanspruchungen von Stahl nicht so gut erfüllt, wie es W. Lode bei seinen Versuchen 
mit festgehaltenen Hauptachsenriehtungen fand. Auch wenn man den Einfluß der Verfestigung 
zu eliminieren versucht, bleiben merkliche Sehwankungen in den beobachteten Fließ- 
bedingungen, die auf Anisotropien infolge der vorangegangenen Verformungen hinzuweisen 
scheinen. Trotzdem wird der Zusammenhang zwischen Spannungen und Verformungen durch 
einen neuen Ansatz für ein isotropes Kontinuum annähernd wiedergegeben. Dieser Ansatz 
enthält als Sonderfälle die Ansätze von Saint-Venant-Mises und Prandtl-Reuß. Es 
sei noeh einmal darauf hingewiesen, daß der neue Ansatz lediglich eine formale Beschreibung 
der Versuchstatsachen in idealisierter Form gibt. Eine Erklärung der von uns gefundenen 
Erscheinungen müßte natürlich an den Aufbau des Materials aus anisotropen, verschieden 
orientierten Elementen anknüpfen. Es wäre zunächst zu untersuchen, ob die hier mitgeteilten 
Beobachtungen auch bei anderen Metallen und Einkristallen zutreffen. 172 
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Fließpotential oder Gleitebenen? 


Von A. Reuss in Budapest. 


1. Verallgemeinerung der Fließpotentialtheorie. Die Flielspotentialtheorie von R. v. Mises'). 
welche mit den Versuchen gut übereinstimmende Resultate ergibt, läßt sich so verallgemeinern. 
daß sie neben den plastischen Formänderungen auch die elastischen berücksichtigt. 


Die Grundlagen der Fließpotentialtheorie lassen sich in der verallgemeinerten Fassung 
wie folgt aussprechen: Die Deformationsgeschwindigkeiten des plastischen 
Körpers unter Einfluß eines Spannungszustandes, der an der Fließgrenze 
liegt, regeln sich derart, daß sie bei einer unendlich kleinen Variation der 
Spannungen innerhalb der Fließgrenze keine zusätzliche plastische Form- 
änderungsarbeit leisten und das Volumen des Körpers nur elastisch ver- 
ändert wird. Die Berücksichtigung der elastischen Formänderung neben der plastischen 
macht es notwendig, daß die zusätzliche Arbeit genauer als in der R. v. Misesschen Ab- 
handlung, als die plastische Formänderungsarbeit umschrieben werde, ferner, daß die 
Forderung nach der Art der Volumenänderung aufgestellt werde. Eine ähnliche Forderung 
war auch in der ursprünglichen Misesschen Fassung enthalten, aber nur implizite, in 
den Bedingungen, welchen die Fließgrenze zu genügen hat. Dadurch, daß wir sie in den 
Grundsatz übernehmen, kann der Einfluß eines überlagerten hydrostatischen Druckes auf 
die Fließerenze, auf den bereits R. v. Mises hingewiesen hat, leicht berücksichtigt werden. 


Die Fließbedingung können wir ganz allgemein als eine Gleichung 
PIE ee De Pi ee 


zwischen den Spannungskomponenten 6%,0y..:.. ’. annehmen, von der bloß gefordert wird, 
daß sie gegenüber 'Transformationen des Koortdıinatensysteims, welche «das Kristallsystem zu- 
läßt, invarıant sei. 


Die Komponenten der elastischen Formänderung berechnet man aus dem elastischen 
Potential 





| 
| .) S 10x 20x Oy 513 0x 9 ra So 0x T 
| [77 
.) S,0y AR Öy Ö r 85,,0y T 
I) 
| i rw (2), 
= # Öö- -.. .. + N, (7) T 
ER 
| 
p ») arT I | 
wo die s;7; 8x; Materialkonstanten sind, zu 
OU ‚ BU OU 
Ki . ey e. eerret Y (3), 
00, 00, Zu ÖT, 
Es &ys se) bezeichnen die elastischen, &,, &,, ...y- die gesamten Formänderungskom. 


ponenten. Die Bedingung, daß die Volumenänderung nur elastischen Ursprungs sei, wird: 
et Br tt Bat 85) Oy TFT I Bu tat, (4) 


Die Leistung der plastischen Formänderungsarbeit soll als die Differenz der Leistung 
der gesamten und der elastischen Formänderungsarbeit angeschrieben werden: 


de, dex 


dt ct | 9 


| ] n (“& de, dy, dy, 

4 4 » I a ER Yo 

2 9 " At dt BZ Fl «At 
Wird nun laut unserem Satze L, unter Berücksichtigung der Nebenbedingungen (1) 

und (4) variiert, so ergeben sich unsere Grundgleichungen der Fließpotentialtheorie: 


ıı RB. v. Mises, Meehanik der plastischen Formänderunge von Kristallen. Diese Zeitschrift. Bd. 8. (1028), 
S, 161 bis IN.». 
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de, de, Oo F | 
/. MIE.. 8s..— 8...) v) 
lt lt 06, ni .. 
de, de’ .231 
/. ut (Ss S - 8,..) (0) 
lt At O6, di . 
/ ee a 
dd dv. 04 
/ mt |(s -8..+8.)=0U 
dt lt Or ul ' aid 


wo z und « von den Spannungen unabhängige Multiplikatoren sind. 


ls soll ausdrücklich hervorgehoben werden, daß bei der Variation sowohl die gesamten, 


wie «die elastischen Formänderungsgeschwindigkeiten unvarliert bleiben. 
Aus den Gl. (1), (3), (4), (6) lassen sich bei gegebenen Formänderungsgeschwindigkeiten 
(die elastischen und plastischen Anteile derselben, ferner die Spannungen berechnen. 


Die Zeit # kommt ın den Gl. (6) als eın Parameter vor. welcher. da er bei der Be- 


rechnung herausfällt, ohne den Sinn der Gleichungen zu ändern, durch eine beliebige Zeit- 


funktion ersetzt werden kann. Die Grundgleiehungen gelten übrigens streng genommen nur 


für unendlich lanesame Formänderune. 
Wenn man die drei ersten Gl. (6) addıert und mit (4) vergleicht, so entsteht 


‚öF dF OF 5 
A|l-- P3 ”.\ Hul(s,, PB. Tr 8) 7a lb +8s,,+8,)]: 0. 


f} { Oy (0 


( „N 


Drückt man ın Fo,. 6,, 6. dureh 


‘) (2) 0 
\ 7 . , an 5 
pP’ > - Or O0, P» Oy OÖ, y 
aus, so wird 
’ h ’ P ‚ Pe 
Or Or pP O0 Oy PP» 0. Ox Oy rp 
Of OF Of oJ 
‘) 
Op 06, (00, Rz ( 
und 
‚of 
S Op | 
u ee er a ET 
8,; 8. 7 8,., 7818... 7 8,, 8 


Man sieht also unmittelbar, daß, falls die elastische Formänderung neben der plastischen 
vernachlässigt werden kann und die Fließgrenze vom hydrostatischen Druck p unabhängig 
ist, die Gl. (6) in die R. v. Misesschen Gleichungen 

de, ‚oF de, ‚of dy- ‚OF 
/ / NR / 
lt 00 lt 06, lt OT. 


\ 
übereehen. 


Für ısotrope Körper mit der Fließbedineune 


(6, 6.) + (0 0) +10, — 0,” +6(1,.’ +7, + )=flD - » «: » .:..). 


wo «die reehte Seite allein von p abhängt, erhält man 


ds: \ dp d z | dyx ‚Ar, ’ 
(Ad » (6, P) mATO, p): 5 / T ATx | 
dt lt dt 2 dt lt ' (10). 
rar ) 
j a m — 2 l 0. Ai 
die bereits bekannten Gleichungen), wo _, ‚a und ,„,„‚=Pp gesetzt wurde. Die Kom- 
>(m+1)G 2G 


ponenten der plastischen Formänderungsgeschwindigkeit können auch als Ableitungen der 
in Gl. (9) links stehenden quadratischen Form aufgefaßt werden, so daß in diesem Falle das 
Fließpotential von der Fließbedingung in der rechten Seite von Gl. (9) abweichen würde. 
Auf die Möglichkeit, daß das Fließpotential und die Fließbedingung voneinander ver- 


\. Reuss, Berteksiehtigenng der elastischen Formänderung in der Plastizitätstheorie. Diese Zeitsehr. Bd. 10. 
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schieden seien, hat bereits R. v. Mıses ın seiner Abhandlung hingewiesen. In der ver- 
allgemeinerten Fließpotentialtheorie scheint es jedoch zweckmäßig, das Fließpotential mit 
der Fließbedingung identisch anzunehmen. 


Da zur plastischen Formänderung unbedingt Arbeit aufgewendet werden muß, hat man 
für das Auftreten der plastischen Formänderung außer Gl. (1) noch die Ungleichung 


ln 0 . . . . . . . . . . . . . . . . . . ° . . . . . (11). 
Multiplizieren wir die Gl. (6) der Reihe nach mit 0x, O4 .:... r. und addieren, so wird mit 
Rücksicht auf (5) und (8) 
of, L 
(° F | A) F oO F | h) F Op (Ex Ey € ) 0 ) 
de / u 06,7 0... T- z " . =). 
2 06x ” 00, “ 06; OT, Su trat 8; + < (St 85 7 812) 


Für isotrope Körper mit der Fließbedingung (9) wird speziell, da die linke Seite eine 
uadratische Form ist, 


ARE # % 5: , Me PR Ta N Eee Er re ee 54° 
Die linke Seite von (9) muß positiv sein, f(p) ist also ebenfalls positiv, infolgedessen muß 


te 
sem. 


2. Die Gleitebenentheorie. Die plastische Verformung wird bei. vielen Autoren als ein 
(leiten entlang ganz bestimmter Ebenen aufgefaßt. Diese Ansicht ist für Kristalle in vielen 
Fällen durch Versuche bestätigt worden. O0. Mohr’) hat für isotrope Körper die Existenz 
von Gleitebenen postuliert und die Hypothese aufgestellt, daß das Gleiten in denjenigen 
Ebenen auftritt, in welchen die Schubspannung 7 zuerst einen von der Normalspannung o 
abhängigen Maximalwert = F(o) erreicht. Aus dieser Hypothese folgt nun die bekannte 
Mohrsche Fließbedingung, welche in den Hauptspannungen o, > o, > o, ausgedrückt 


Beer ee re ern Kern 
von der mittleren Hauptspannung o, unabhängig ist. 


Es ergeben sich zwei Gleitebenen, die zur Achse 3 parallel und zu den Achsen I und 
2 symmetrisch liegen. Aus den Gleitungen entlang dieser Ebenen resultiert eine Dehnung 
in der Richtung 1 und eine Verkürzung in der Richtung 2. In der Richtung 3 tritt keine 
plastische Längenänderung auf. 


Auf Grund der Gleitebenentheorie hat L. Prandtl bereits in seinem Vortrag anı 
Delfter Kongreß für angewandte Mechanik *) im Falle des ebenen homogenen Problems die 
plastisch-elastische Formänderung erschöpfend dargestellt. Im folgenden will ich seine Aus- 
führungen bloß durch die Berücksiehtigung der mittleren Hauptspannung ergänzen. 


Die Richtung 3 falle während der ganzen Formänderung immer mit der Riehtung z 
zusammen. Auf die Spannungshauptsachen bezogen sind dann die Spannungen o,, 0,, 0, 0, 
0, © und der plastische Anteil der Formänderungsgeschwindigkeit 4, - 4, 0, 0, 0, 0. Bringt 
man die Achsen 1, 2 durch eine Drehung vom Winkel a um die Achse 3 mit den Achsen «., 
y zur Deckung, so werden in diesem Koordinatensystem die Spannungen 


=6, cos’a E 6,81" a, 6, G, sın? a = 0,608’ a, Ö; Ö. .=T,= v0. T, (0, 0,) sın a cosa 


Ö 33 


A 


und der plastische Anteil der Formänderungsgeschwindigkeit 


&, —=ÄAcosBa, &y = —AcosBa, 8" =P0, y," Yy VÖ, y; =2Asin2a, 
welche mit Rücksicht auf 0, 0,= (6,  0,)cos2a und nach Unterdrückung des gemein- 
samen Faktors 1:(o, — o,) in 

"Pie == / (6x Öy)s Pt / (6x Öy), Pu 0, Pe ei), vu 0, y '— | N F, 


übergehen. 


3) 0, Mohr, Abhandlungen aus dem Gebiete der technischen Mechanik. 1906, Abhandlung V. 

4), L. Prandtl, Spannungsverteilung in plastischen Körpern. Proceedings of the first International Congress 
for Applied Mechanies. Delft, 1924, S. 43 bis 46. Dieser Vortrag war mir beim Erscheinen meiner Abhandlung über 
die „Berüeksichtigung der elastischen Formänderung in der Plastizitätstheorie‘ nieht bekannt. 








Ztsehr. f. angew. 
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Setzen wir die gesamte Formänderungsgeschwindigkeit gleich der Summe von elastischer 
und plastischer Formänderungsgeschwindigkeit, so treten an die Stelle unserer Gl. (10) ın 
der Gleitebenentheorie die Gleichungen 


de. dp ‚d 


I 


dar trPail'& pP) + 4x 0) 


de, dp ‚d 


} 





‚lt d dir » dt (0, pP) / (Ox O,,) | 
(16). 

de. dy er; 

7er Ta a7 TR zu 

dy, Zr d (; 

"ill dr 7 

Die Fließbedingung muß gleichzeitig durch 
(a, II SFT: 7 3 5 Eee ehren (17) 


ersetzt werden. 


Gehen wir von dieser Gleichung als Fließpotential aus und wenden wir die Gl, (10), 
(7) und (S) an, so entsteht mit Unterdrückung eines belanglosen Faktors 2 in 4 


de, dp d . / df(ox + 6,,) 

dt ar PT pP Tran) TE lo, -+06,) 

de, dyp ‚d eu | ), df(ox + 6,) 

dt ” dt "dt u u > NO; ey 6 d(o,—+9,) | 

(19). 

de, dy ' d ı df(ox + 9,) 

dt mr dt rPat (9 P) s 3 d(6.—+ 0,) 
dy, 2 d 4: 

dt 0 u 





Wie wir sehen, stimmen die Gl]. (16) und (15) im allgemeinen nicht überein, die Flie&- 
potentialtheorie widerspricht also, wenn f(o,—+ 0,) nicht konstant ist, der Gleitebenentheorie. 

Ist o, einer der beiden anderen Hauptspannungen, z. B. o, gleich, so wird dieser Zu- 
stand nach Th. v. Kärmän im Gegensatz zum bereits behandelten „halbplastischen“ „voll- 
plastisch“ genannt und die Formänderungsgeschwindigkeiten sind nach der Mohrschen Auf- 
fassunz unbestimmt. Die Gleichungen verlieren ihre Gültigkeit, da bei der Variation von (5) 
außer (17) nun auch die Nebenbedingung eo, — o, bzw. 


u n a | 2 ( 
6x6, (öx—+ 06,) 0: = Ö- T. 0. . i . . . . . j r (19) 


berücksichtigt werden muß. Der rechten Seite von (15) muß man also der Reihe nach noch 
die Glieder # (0, — 0.), #(0x — 0.), — #(0x + 0, —- 20,), 2z# Tr, zufügen. Wie man sieht, bleiben 
die Formänderungsgeschwindigkeiten in diesem Falle auch nach der Fließpotentialtheorie un- 
bestimmt. 


3. Zugversuch mit einem tordierten Rohre. Die vorher entwickelten Theorien wenden 
wir auf den Zugversuch mit einem bis zur Fließgrenze verdrehten unendlich dünnwandigen 
Rohre an. Die „-Achse legen wir in die Zylindererzeugende, die y-Achse in die Tangente. 
die -Achse in die Normale. Da 7.=r,=0 ist, bleibt o, während des ganzen Versuches 
Hauptspannung, und wie es auf Grund der Endergebnisse leicht bestätigt werden kann, immer 
die mittlere. 

ls stehen zwei Fließbedingzungen einander gegenüber, die Mısessche, Gl. (9), und 
die Mohrsche, Gl. (17.) Die rechtsstehenden Funktionen in denselben denken wir uns in 
Potenzreihen entwickelt und nach dem quadratischen Gliede abgebrochen. Die ersten drei 


Koeffizienten können durch die Fließgrenzen für Zug +, Druck o. und Schub r, aus: 
vedrückt werden. Man erhält so für (9) 
3»\| 3» 
,— +... +6 dert | || I + | -1sp . u DR DEN 
5 OA oO 
und für (IN) 
R i Or —+ 0,,\ Ort 0,\., £ 
(0 4 7% kz,11 - 1 . — ”\ r (0; U 5 (21). 
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In unserem Beispiele ist o,= 0, —=0), 7. = Tr, =0, ye = Yyy=0, yz = konst., ferner schreiben 
wir zur Abkürzung o,=0o und r,=r. Mit diesen Werten geht sowohl (20) wie (21) in 


z £ m 0, 
E = r:|1 = I + . ea 5 1. 
“ ( 


über. 
Betreffend der Formänderung stehen drei "Theorien einander gegenüber, und zwar 
1. die Fließpotentialtheorie mit der Misesschen Fließbedingung, 
2. die Fließpotentialtneorie mit der Mohr schen Fließbedingung, 


3. die Gleitebenentheorie mit der Mohrschen Fließbedingung. 


Für das behandelte Beispiel wird laut der Fließpotentialtheorie mit der Misesschen 
Fließbedingung aus Gl. (10) 


de, a+2pdo % 2 . 
a 3 5“ 
de, de, a—ßdo |1, ” 
7 F Audi” > Muemnder "aan: >’ Yamada EHRE, BE Sa BE SEE See See Be a Se SE Eu | 
dr E 
2 
VB rt 





laut der Fließpotentialtheorie mit der Mohrschen Fließbedingung aus Gl. (18) 








de, a+2Ppdo , 5 u ER 
: o+— N ı, +20 — /0 
di EBRTHFTTTETEN I 
de, a—ßPdo , 2 Am Ä e E , 
; A Ö oO -ö0 > 0 
dt 3 dt Dee ee ee 37 
| (24) 
ee, a-Paa Bir, ER 2, 
er Em Inner TI TFT 
dr 
0 — +24 
P dt e 
endlich laut der Gleitebenentheorie aus (16) 
dezz a+2ßdo ,„ de, a—pPdo 
dt 3 dt dt 3 di i 
(25). 
de, a—-ßBdo ,„ 0 we. i; 
dt ar ut a T Te 
Differenziert man die Gl. (22) logarithmisch nach der Zeit, so wird 
I dr l / En - Il \do Br 
=: geek a ae (26). 


Drückt man A aus den letzten Gleichungen einer jeden der Gleiehungsgruppen (23), (24), (25) 
. . . ı / f} .. ® a do .. « 
aus und berücksichtigt die Gl. (22) und (26), so läßt sich A durch o und lt ausdrücken. Setzt 
man diesen Wert von 4 in die ersten drei Gleichungen der Gleichungsgruppen (23), (24), (25) ein, 
so können die so erhaltenen 9 Gleichungen ohne weiteres integriert werden. Die Konstanten 
bestimmen wir so, daß für 6=0 +. =£:,=e,=0 werde Führen wir die Abkürzungen 





0 \o_ 
Me ı- 
«= log- —— 7 
Fe) O T 
1-2) 
Oı 
(27) 
0 
2 l+ 
1) Ö- 
v=5 (6 ,-0-)log bo 
- 0, 3 Ö 
| 
O0; 
>” 
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ein. so ist das Resultat der Integration in der folgenden Tabelle zusammeneefaßt: 


Fließpotentialtheorie mit 
Misesscher Mohrscher Gleitebenentheorie 
Fließbedingung 
7 aD 5 at=sp n ataöp | 
\ | z [77 ‚pP t | rn O a (n ") » Ö IL u 
F +.) +.) - u 
) ) ) 2 Im 
(——.£ lo u >" \o+>(-2u+r) | _. o - u ER 
2 > gr 3 a Del, 3 a 7. 
a B | a pP pP 7 pP 
- 0 u ; Z (" 2v) | > 
| | | , ! | n | ) 5 
ıdex| a--2p9 a+-2B B (0 04) _ a+2p 
fi 0) “ Fa . 
atic—a 3 3 6 014°0_ 3 





Das erste Glied ist in Jedem Ausdrucke die elastische, das zweite die plastische Form- 
änderung. In der letzten Zeile ist noch die Tangente im Anfangspunkte angegeben. r kann 
aus (22) als Funktion von o berechnet werden. 


Für den speziellen Fall, daß in den Fließbedingungen (20), (21) die rechten Seiten als 
konstant angenommen werden, daß also 6. : o_=yV3r,=k bw. =2r,=k, wird 


=. 
" Zr—hkAUNG, re | ° 


Die Fließpotentialtheorie mit Mohrscher Fließbedingung und die Gleitebenentheorie werden 
identisch und Tabelle (28) geht in 


Misessche Fließbedingung Mohrsche Fließbedingung 
ao 2 Mr 3 o a6 E. | Ar 3 ö 
€ 2 ‚ ) 3 Tzıc . a ” DO+-T DK P z zZ 
Sb u %, I  Tarerr I, 
a0 | Mr 3 0 a0 | | kr 3 7 GE 5; 
E, : r AT = Tu P0 Hk AU2g- 
HJ 3 3 p g IR 3 6 p : I / g IK 
(LO | Oo AO | 
) » \) * > Zu. 
3 z Pr Ag, 3 Er 





über. 


Die Gl. (25) eignen sich für eine experimentelle Prüfung der oben entwickelten Theorien, 


Verkürzung des Zylinderdurehmessers gleich der der Wanddicke wird, wogegen laut der 
(leitebenentheorie die Wanddieke sich bleibend überhaupt nieht vermindert. 


4. Genaue Berechnung für endliche Wanddicken im Falle der Gleitebenentheorie. Die im 
vorigen erhaltenen Resultate sind mit Ausnahme der Fließpotentialtheorie mit Misesscher 
Fliesbedingung streng nur für unendlich dünne Rohre gültig. Um einen Anhaltspunkt für 
die Brauchbarkeit der hergeleiteten Gleichungen auch für endliche Rohrdieken zu haben, 
habe ich für die Gleitebenentheorie auch die genaue Berechnung durchgeführt. 


Wir führen Zylinderkoordinaten ein. Es bestehen die Gleichungen 


(rer) Oo (r o,) 

€ . Of _ 

or or 

oder wenn man den Operator 

Q Ö 

I (r)=1-+r (31) 

I ö ° . ° ° . . . ° . . . . . 

hi " 0 


einführt, die Gleichungen 


da laut der Fließpotentialtheorie mit Misesscher Fließbedingung die bleibende spezifische 
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Das Rohr denken wir vor dem Zugversuch so weit verdreht, bis die Fließerenze im 
Rohr überall erreicht wird. Die Gleiehungen, welche den Formänderungstensor mit dem 
Spannungstensor verknüpfen, lauten nach «16) 





ÖE: Op, „oO Ä | 
7 —- N Ö- I») --AIl0. O0 
A; Ay Beau Az Ei 
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=—=n N (6 ) AND, 04) 
A; af "rt ! ’ 
a ve (33). 
DE, Op ® | » 
(A -- 1) Ö ) 
A ur 
A Or 
aAI 
Pi 
Die Fließgrenze nehmen wir laut der Schubspannungstheorie als 
EEE .. 5: 5 eier Bee 
an. 
Subtrahiert man von der ersten der Gl. (33) die zweite und setzt zur Abkürzung 
eher: ee a ee ee 
So wird 
\ \ 
( (O0 n 
(£, Ey) 1 re te re 
of a | of P 


Differenziert man (34) nach der Zeit, so kann man aus der so erhaltenen Gleichung so- 


0) we. 
wie aus der letzten Gl, (33) mit Rücksicht auf (34) und (35) r und * eliminieren und erhält 


‚00 
2A 2 „4 | (37). 
k?— 0° 
Diesen Wert in (36) eingesetzt und integriert wird 
En ee te ee a ee sin 
wo zur Abkürzung 
Ar 3g ee et a a en x I 
gesetzt wurde. 
Aus der Summe der ersten drei Gl. (33) erhält man 
re : 3 en a in nd AS (40) 
und aus der dritten Gl. (33) 
eur the —-D : -» : >» 2 Hr 0 0 8 eine RR: 


Die Integrationskonstanten wurden alle so bestimmt, daß am Anfang des Zugversuches alle 
e und alle o verschwinden sollen. 


Der Operator D auf (35) angewendet, ergibt mit Rücksicht auf (32), da «, von r unab- 
hängige ist, 


en. N a ee 
Addiert man einmal (41) und (42), so wird 
,=la-Dry+i, +PREDB „..- 202 00.8. 805 Ü 
dann (38) und (40) und berücksichtigt (41), so wird 


2 € (2 u P) yp p Ort pP a I A ER (4 t). 
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Setzt man in den letzten beiden Gleichungen 


| 
p 3 (0, +2 Do,+0) 


ein, so folgt 
38, (a+2P)o,+-la—-Pp)(2Do,+o)+3ßPkDR | 
be: 2la- ), +Pa+P)(@Do,+o)+3ßkR | 
2 Do,-+o läßt sich eliminieren, 
3e,=3a0o,+(2a+ß)kEDR—(a pP)kR 
Wendet man auf diese Gleichung den Operator D an, so ergibt sich 
3e.=3aDo,+@a+p)kD’R-(a- PEDR. 


da De.—r,. Setzt man o, und Do, aus (47) bzw. (48) in die zweite 
berücksichtigt, daß laut (39) 


P= k Ig R. 
ist, so erhält man für R die Differentialgleichung 
2(2a+p)D’R (a+2P)R BalgR=VvV 


Aus (47) erhält man die Randbedingune an der inneren und Äußeren 
/ylinders, da dort o,=0 ist, zu 


(45) 


(46). 


(47). 


(48), 


(]. (46) ein. und 
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(0). 
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(39). 


Randbedin- 
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"(da Tr” P) DR (a pP) I 
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Wegen dem Gliede mit Iaq R ist Gl. (50) nicht linear. Wir setzen R 
nr 
reihe nach & = "an. Aus 
R=-a+ta.+a,°+a, + 
folgt 
DR (a,+a)+2la, +a)E+3(la,+a,)&P+... 
D’R - (a,+3a, +2a,+2(a, +5a,+3a,)E-+... 
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. ü a 
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a Q 0 Goj? Ad, 
Bezeichnet r; den inneren, r„ den äußeren Halbmesser des Rohres und wählen wir 
tr ’ r r or 
„=  , „setzen ferner zur Abkürzung d= " . "so geht 
ö N" N", 


gung (Al) an der äußeren Zylinderfläche für &=d in 


IE. 
IR A, 


über und dementsprechend an der inneren Zylinderfläche für 


. 98 
ge=—onm A, —A 


z ö+4,0° 


1 


über, wo die Koeffizienten A nur a. P und die zu bestimmenden Koeffizienten a enthalten. 


Die Summe und Differenz der letzten beiden Gleichungen ergibt 


VA, +...=3la+P)a, +22a+ß) a, 
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Wir unterdrücken von A, an alle höheren Glieder. Zwei weitere Gleichungen erhält man 
zur Bestimmung der unbekannten Koeffizienten a,, a,. a,, a,, indem man die Entwicklungen 
(52). (54) und (55) in die Differentialgeleichung (50) einsetzt und die so erhaltene Entwicklung 
identisch eleieh Null setzt. Die ersten beiden Glieder erzeben 


2(2a-+P)(a,+3a,+2a,)- (a+2P9)a,—- 3alga,V | 


d, . . . . . (AT). 
2a (2a. +5a.+3u,) (a2) a 3a = 0 
/ 1 . + 3 PB 1 Go}? dt, 
Die (Gl. (56) und (57) genügen, um zu Jedem gegebenen e, die dazu gehörenden Koeffizienten 
von R:a,, a, a, a, zu bestimmen. Mit Hilfe der Entwicklung von R erhält man aus (47) 
drei Glieder der Entwieklung von o,, ferner aus 0; = Do, zwei Glieder von o,. Auso=kIgaR 
berechnet man 0.—=0o-+-0; für jede beliebige Stelle im Rohr. (35) und (42) ergeben e; und e,. 


Die durchgeführten Berechnungen sind nur so lange gültig, bis o, tatsächlich die 
mittlere Hauptspannung ist. Da die beiden anderen Hauptspannungen durch 


FE | 
ÖL,» . .) + .) I; 


vegeben sınd, so wird diese Bedineunge an der äußeren Oberfläche 
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Abb. 1. 


In der Abb. 1 sind die Werte von o, an der inneren Oberfläche, in der Mitte der 
Wand und außen, 6; innen und außen, endlich 7 innen, in der Mitte und außen zur Abszisse 
e, aufgetragen. Die Spannung o, ist so klein, daß sie im Maßstabe der Abbildung kaum 
dargestellt werden kann. Dabei wurde a«=0,25 - 10°" cm’/kg, P = 0,625: 10 * em’/kg, 901, 
: = 3000 ke/em? angenommen. Zum Vergleiche ist auch der Wert von o, nach der Fließ- 
potentialtheorie mit Misesscher Fließbedingung gestrichelt aufgetragen worden. 


Aus der Abb. 1 sieht man, daß unsere Gleichungen in unserem Beispiele nur bis zu 
einer Dehnung von rund 0,21 v. H. gelten. Von dieser Dehnung entfallen bloß 0,065 v. H. 
auf die bleibende Dehnung. 

Für eine kleinere Wanddicke wird zwar diese Grenze hinausgeschoben, dem gegen- 
über werden die aus der ungleichmäßigen Wanddicke rührenden Versuchsfehler vergrößert. 
Es scheint also als ziemlich schwierig, experimentell festzustellen, ob sich die Wanddicke 
(des Zylinders bleibend verändert oder nicht, um so zwischen der Fließpotentialtheorie und 
der Gleitebenentheorie zu entscheiden. 
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3. Entscheidung zwischen Fließpotential- und Gleitebenentheorie auf Grund der Rauminhalts- 
änderung. Da eine direkte Bestimmung der Wanddiekenänderung unseres Versuchsrohres sich 
als ziemlich schwierig herausgestellt hat, scheint eine Methode, welche diese Schwierigkeiten 
umgeht und so zur Entscheidung zwischen den in Rede stehenden Theorien führt, mehr Er- 
folg zu verspreehen. Als eine solche Methode bietet sich die genaue Bestimmung der Änderung 
des Rauminhaltes unseres Zylinders. Man denke sich den Zylinder in den Einspannköpfen 
verschlossen und mit einer Flüssigkeit gefüllt, deren Niveauänderung in einer kalıbrierten 
Kapillare beobachtet werden kann. Das Prinzip ist dasselbe wie beim Flüssigkeitsthermometer, 
mit dem Unterschiede, daß jetzt nieht die Änderung des Flüssigkeitsvolumens, sondern die 
des Gefäßinhaltes bestimmt werden muß. Die spezifische Änderung des Rauminhaltes beträgt: 


| 
| Ex er u at „ . ö . . . . . . . . . P . . ° . (SH), 
Mit den Werten von Gl. (25) erhält man für die Fließpotentialtheorie mit Misesscher Fließ- 
bedineune: 
l=a0 . EEE Yo PROBE aE? 7 ESGER | Dr Su ee ie S 
mit Mohrscher Fließbedingung: 
= -Blua—v). . : >.» ee De 


endlieh für die Gleitebenentheorie: 


| 106 5 re nr 


Mit den speziellen Werten (29) bleibt Gl. (60) bestehen, während (61) und (62) in 
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\bb. 2. Abb. 3. 


In Abb. 2 sind zur Ordinate o die Abszissen | aus (60) und (63) aufgetragen. Der Ver- 
lauf der Kurven ıst grundverschieden und eignet sich zur experimentellen Untersuchung. 


Um den Mechanismus der Formänderung ohne Rücksicht auf die Spannungen zu unter- 
suchen, wurden die spezifischen Rauminhaltsänderungen auch als Ordinaten zu den axialen 
Längenänderungen e, in Abb. 3 aufgetragen. e, wurde aus den beiden entsprechenden 
(1. (30) berechnet. Die Abbildung enthält auch die Tangente im Anfangspunkt und die 
Asyınptoten. 2341 
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Achsialsymmetrische, freie Flüssigkeitsstrahlen mit schwacher 
Kontraktion. 
Von H. Reifner, Berlin-Charlottenbureg. 


1. Aufgabestellung und literarischer Rückblick. Das räumliche Problem des freien Strahles 
in einer inkompressiblen Flüssigkeit tritt in der Theorie der Ausflußstrahlen aus Gefäßen und 
in der Theorie der Propellermechanismen auf. Die hauptsächliche mathematische Schwierig- 
keit besteht in der Erfüllung der Randbedingung des Gleichgewichts mit der umgebenden 
Flüssigkeit an einer nicht gegebenen, sondern erst durch die Untersuchung zu bestimmenden 
Oberfläche. 

Dieses Gleichgewicht verlangt für einen Strahl aus ruhendem Ausflußquerschnitt kon- 
stanten Druck an der Begrenzungsfläche, für einen Strahl mit beweztem Ausflußquerschnitt, 
wie z. B. von einem bewegten Fahrzeug aus, streng genommen eine kompliziertere Rand. 
bedingung, da an einem sich kontrahierenden und bewegten Strahle die äußere Atmosphäre 
kleine Radialgeschwindigkeiten aufweisen muß. Jedoch ist die Abweichung von der Be- 
dingung konstanten Druckes bei der im folgenden vorausgesetzten Beschränkung auf kleine 
(Größen erster Ordnung klein von höherer Ordnung. 

Räumliche Potential-Ausflußstrahlen sind in Arbeiten von M. Plancek' und P. Moelen- 
broek’), J. Boussinesq’), Fr. Kötter’), T. Levi-Cıivita°®), E. Trefftz*) und R. Förster’) 
behandelt worden. 

Die beiden Erstgenannten, Planck und Moelenbroek, führen das achsialsymmetrische 
Problem auf die ebene Potentialströmung zurück und geben allgemeine Sätze über Niveau- 
flächen und die freie Oberfläche. Boussinesq, Kötter und Levi-Civita haben ver- 
schiedene Abschätzungsmethoden für die Kontraktionszahlen entwickelt. Trefftz hat eine 
ins einzelne gehende Darstellung des achsıalsymmetrischen schwerefreien, wirbelfreien und 
zirkulationsfreien Strahles aus einer kreisförmigen Öffnung im Boden eines großen Gefäßes 
auf Grund einer von ihm aufgestellten Integralgleichung dureh ein schrittweises Verfahren 
durchgeführt. 

Während sich die vorher genannten Autoren mit bestimmten Ausflußstrahlen aus Ge- 
fäßen beschäftigt haben, hat sich Förster mehr dem allgemeinen Aufbau von Potential- 
strahlen mit Hilfe der Integration der Potentialgleichung durch Reihen und Koeffizienten- 
bestimmung durch die Bedingung der freien Oberfläche zugewandt, wobei er auch die Schwere 
und die Kapillarkräfte berücksichtigt. Seine Arbeit und auch die von Planek und Moelen- 
broek, die mir alle übrigens erst nachträglich durch das schöne Referat von F. Auerbach‘) 
zur Kenntnis gekommen sind, haben in dem Ziel und in einzelnen Ansätzen eine gewisse Ver- 
wandtschaft mit der im folgenden mitgeteilten Formulierung, soweit sie sich auf den zirku- 
lationsfreien Potentialstrahl erstreckt. Jedoch scheint mir auch für diesen letzteren Fall im 
folgenden ein Fortschritt durch eine mathematisch und technisch zweckmäßigere Fragestellung 
erzielt zu sein. 

Das Ziel der hier mitgeteilten Untersuchungen ist nämlich zunächst wesentlich enger, inso- 
fern als sie sich auf Strahlen beschränken, die nur wenig von der Zylinderform abweichen und 
für die im Anfangsquerschnitt die Geschwindiekeitsverteilung oder die Druckverteilung vor- 
geschrieben wird. Für eigentliche Ausflußstrahlen würde es sich also nur entweder um den 
Zustand in großer Entfernung von der Mündung oder um Mündungen von schwach konischer 
Form handeln. 

Andererseits ist aber das Ziel der Arbeit erheblich weiter gestellt, da sogenannte Pro- 
pellerstrahlen betrachtet werden sollen. Tatsächlich handelt es sich hierbei in den meisten 
Fällen zwar um Strahlen mit schwacher Kontraktion, jedoch ist es wesentlich, daß Zirku- 
lation, Rotation und Wirbelverteilung vorschreibbar bzw. durch einen Propellermechanismus 
erzeugbar sein müssen. 

Nun hat allerdings die neuere Propellertheorie ®) die Vorstellung des achsialsymmetrischen, 
scharf abgegrenzten Strahles als Grundform verlassen und arbeitet mit Potentialströmungen 
von sektorförmigem Querschnitt, die von schraubenförmigen Wirbeltlächen (Diskontinuitäts- 

1, Wied. Ann. 21, 409. IS84. 

2, Wied. Ann. 52, 207. 1894. 

3) Comptes Rendus, 70, 33, 117, 1279. 1870. 

1) Arch. d. Math. u. Phys. (2) 5. 1887. 

5), Comptes Rendus, 157, 481. 1013. 

6), Dissert. Straßburg, 1913. Zeitschr. f. Math. u. Phys. 64, 34. 1016, 


”) Zeitschr. f. Math. u. Phys. 62, 319. 1014. 
»), Handb. d. Phys. u. Techn. Mech. (Hort-Anerbach) Bd. V, SS. 22) 


», A. Betz, Schraubenpropeller mit geringstem Energieverlust. Mit einem Zusatz von L. Prandtl, Göttinger 
Nachr. Math.-Phys. Kl. 1919. S. 193. Nendruck in L. Prandtl und A. Betz. Vier Abhandlungen zur Hydrodynamik 


und Aerodynamik, Göttingen 1997. 
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flächen der radialen Geschwindigkeit) eingeschlossen sind und die keinen scharf abgegrenzten 
Strahl bilden, sondern stetig, wenn auch sehr schnell in die umgebende Flüssigkeit übergehen. 
Dabei entspricht die Zahl dieser Diskontinuitätsflächen der Zahl der Schraubenflügel. Jedoch 
hat sich auch in dieser Theorie der Begriff des achsıalsymmetrischen Strahles, der sich er- 
eibt, wenn zu unendlich vielen Flügeln und damit zu einem scharf abgegrenzten Strahl über- 
gegangen wird, als praktisch erwiesen. 


Dieser modernen Theorie ist bisher nur der Aufbau von Strahlen unter Vernachlässigung 


der Kontraktion zelungen. Die eigentliche Veranlassung zu den foleenden Stu- 
dien istnundasZiel. die Umsetzune des Überdruckes undderGeschwindie- 
keitsverteilung. welehe dureh die Strahlkontraktion im Gleichgewicht mit 
der Umgebung bewerkstelliet wird, ım einzelnen, wenn auch im nach- 
stehenden nur für achsiale Symmetrie zu verfolgen. Es zeigt sieh dabei, 
daß man für eine Reihe von gegebenen zylindrischen Ausgangs- oder 
Grundströmungen mit und ohne Geschwindigkeitspotential und mit und 
ohne Zırkulation bzw. Rotation Nachbarlösungen der Geschwindiekeits- 
und Druekverteilune mit vorschreibbarer Strahlkontraktion berechnen 
kann. Dabei ergeben sieh immer dann statt der vollen Strahlquerschnitte 
ringeförmige Querschnitte, wenn infolge der Zirkulation oder Rotation un- 


endliecehe Gesehwindiekeiten an der Achse auseeschlossen werden müssen. 


wobei bestimmte Bedineungen zwischen der inneren und der äußeren Kon- 
traktion erfüllt seın müssen. 

Ührieens erscheint es durchaus nieht ausgeschlossen, daß mit Hilfe des nachstehend 
beschriebenen Verfahrens auch unsymmetrische und unstetige Strahlformen aufgebaut werden 
können. 

Im folgenden ist mehr das Grundsätzliche und sehr viel weniger eine numerische oder 
oraphische Darstellung der Zusammenhänge und Folgerungen gegeben worden. 

Herr K. Hinz, dem ich bei der zuweilen mühevollen Durehfindung bis zu den folgenden 
verhältnisinäßie einfachen Ansätzen wertvolle Unterstützung verdanke, wird in einer weiteren 
Arbeit «diese teehniseh wiehtige numerische Durehführunge und einige der vorstehend er- 
wähnten theoretischen Erweiterungen bringen. 


2. Der zylindrische Strahl der Grundströmung. Die im folgenden behandelte Aufgabe soll 
darin bestehen, zu der bekannten achsialsymmetrischen, gleichförmigen und stationären 


Strömungsform des zylindrischen freien Strahles einer inkompressiblen Flüssigkeit von der 


Diehte v. Nachbarlösungen aufzufinden. 
Jene bekannte zylindrische Strömung wird dargestellt durch die nachstehenden mit den 
willkürliehen Funktionen »,(r) und ww,(r) gebildeten Geschwindigkeitskomponenten: 
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\bb. 1. 


Verlanet man, daß die Strömung wirbelfrei ist, so ist sie mit den Konstanten 


Pr 


und ce, aus dem Geschwindigkeitspotential 
‘ 0 ( 7 _ er U ı) . . . . . . . . . . . . . . . . (2) 
r, 
abzuleiten, hat also die Geschwindigkeitskomponenten 
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° ro» r EEE REN VER RE TEE GE 
Oo Cu 
nn, 
02 7 








Ba 


an 


W« 


od 


eilt 


I 


Id 
VO 


Wi 


re 


wc 
st: 


od 


bi 








Band 12, Heft | ‚ . ö u ra i i „> 
Febrnar 19832 Reißner. Achsialsyvmmetrische. freie Flüssigekeitsstrahlen 27 





ınd die Druekverteilune: 


st. : (>b1. 


W ., 
’o Pa > ev | r: Y. 


Wenn c„=$0 ist, kann die Strömung nicht bis zur Achse gehen, sondern muß innen dureh 
eine nieht wirbelfreie. stationäre und gleiehförmige Strömung oder dureh ruhende Flüssigkeit 
von kleinerem Druck p; begrenzt sein. 


3. Der wirbelfreie, schwach konvergente Strahl. Ks werde als Geschwindigkeitspotential 
angesetzt: 


4 = Un 9 > 
wo unter q, der Ausdruck in Gl. (2) zu verstehen ist und 9, ein Zusatzpotential bedeutet. 


Es muß dann y, der Laplaceschen Gleichung 


04 


() 2 


Ir + N Fr 


oder. wenn man die dimensionslosen Veränderlichen 


einführt. der Gleichung: 


ıo/ dgp,\, %pg, 
o0o ao 3 #2 a PR .. (88) 
genügen. 

Die Ableitung nach dem Winkel # tritt für dieses Zusatzpotential nieht auf, da jede 
mögliche Zirkulation des achsialsymmetrischen Strahles schon dureh das zweite Glied in q, 
vollständig gegeben Ist. 

Es werde nun zunächst als partikuläres Integral von Gl. (3a) eine Lösung angesetzt, 
welche der verlangten Eigenschaft des Strahles, im Unendliehen zylindrisch zu werden, ge- 
recht wird, nämlich: 


p,=e "ja .J, (20) b N,(zo)] I ’ . HH. 


wo J,(20) und N\,(x20) die Besselschen Funktionen nullter Ordnung und «a und 5b Kon- 
stanten bedeuten. 


Aus diesem partikulären Integral kann man offenbar auch die Reihe: 


\' 
Fı u 


l 


en» la, Jo(l&n 0) en Du N, IH o)| . . ’ (ta) 


oder auch das bestimmte Integral: 


TS 
T, \dxe *-la(z)J,(* 0) biz) N,(zo)| a a db) 


“=—=UVU 
bilden. 

Für den im folgenden verlangten Näherungsgrad wird es genügen, mit einem Glied, also 
mit dem Ansatz (4) zu arbeiten. 

Es muß nun gefordert werden, daß die aus Ansatz (4) folgende Strömung ebenfalls wie 
die Ausgangsströmung in Abschnitt 2 freie Stromlinienflächen, d. h. Stromlinienflächen kon- 
stanten Druckes p besitzt, welche durch je einen vorgeschriebenen Kreis, z. B. den Kreis £ 0, 
o — 0, hindurchgehen. 

Zur Erfüllung dieser Forderung ist zuerst die Gleichung der Projektion einer Stromlinie 
auf die Meridianebene # == const mit der Bedingung, daß diese Stromlinie dureh den Punkt 
20, 0 0, hindurchgeht, aufzustellen. 


Die Gleichung einer solchen Meridianlinie lautet: 


Ri) 7 f Ri 47 N 
do u Q 0 Q 0 
dE Mt: 04 0 } Ai 7 f Q 4 N (»,)), 


Q OL Wr 2 
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Begnügt man sieh mit einer Nachbarlösung zur zylindrischen Grundströmung, für die 


Q Fo og, » 
N: e N (6) 
7 dJ,(z 0) dN,(z0) 
ıst, und setzt man Ansatz (4) eın, so erhält man. wenn - «dureh .J,’(z 0) und > 
d(x 0) = d(z0) 
dureh \,(2 0) bezeichnet wird: 
do | Pr N 
- e “la «J, (ko)+bxN, (zo) . . . . ve mal. 
d. Bun . u. 
Ks ergibt sıch daraus: 
| do 
e >» = C, n f r+ . 
H "zlaJ, (ko)+bN, (x o)] 


Die weitere Integration vereinfacht sich, wenn man entsprechend der Näherungsvoraus- 
setzung (6) der schwachen Strahlkontraktion 


. 
1 
) Onso ö, 6 ee ee 


setzt und die Besselsechen Funktionen bıs zu kleinen Gliedern erster Ordnung entwickelt, 
also schreibt: 


\-bN '(xzo 


0 (+0 0 „(1,0 


al ko) + bh N, (2 o)=a., (zo I orlu J (201) +b N. (X 0uo)]- 


Ks entsteht dann die foleende Lösung: 


Kan oax J, (x O,) +hx N. (X 0,) ' 
’ u, \Ö ‚ = aa 4 . 
aJ, (20) +bN, I(z0,)| 2 al, (zo)+bN, (z0,) 


u 0 0 


e 


woher o, abkürzungshalber dureh o, ersetzt ist. 


1,0 


(seht man auch hierin nur bis zu kleinen Größen erster Ordnung, verlangt. daß für (=, 
s=0 wird und nennt man ferner die Strahleinziehung o im Unendlichen (= x)o, , So er- 
oıbt sich: 

7 N (z O,) 5 N, (+ I) | 


Cor 


ai er a  } ii 





oder 7 o_Al1 BT an 


Soll nun «die Stromlinienfläche der Gl. (S) Randstromfläche sein, so muß auf ihr 
die Geschwindiekeit e=e, konstant sein. Es muß also gelten: 


u" —+-rv. + lw Ha, )” 6,” EL EEE ee Ba HM, 


also unter Benutzung der Gl. (2a) bis zu kleinen Größen erster Ordnung: 


IN u .® 
- 0%, Ca N ) er l (9a) 
“ Pr er ur . . ° . . ; e 5 r 6 .«%K 
serz Ö - Cr Cr” 0 
x ‚ Cala ? CH\"[ l 20 
oder: p ‚ze la J.izo) rt l, IV. (zo)| | | | wg 0 3]|- 
erz P 2 C Cu On OÖ, 


Entwickelt man auch «ie linke Seite nach o und schreibt abkürzungsweise J, = .J, (#2 0,): 
N. N,(20,),. so erhält man: 


.) ., “> 1 
-_ X Cala\ C»\ A Cp\ 

es (aJ, + bN,) | | | a” 20| = Bu ey (Oh). 
Car C; Cr Cr 


Ks soll nun erstens 6 0, sein für © x, woraus sich für die Randgeschwindigkeit e„, als 
Bestimmungseleichung eraıbt: 


(Ida). 







































TS. 


W 
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Ferner soll # — 0 sein für ©=0, so daß man erhält: 


PA, 6% 
(aJ, +tbN)=2o | |o, ee EEE FE, ET TE 
Cw Cu - 


Schließlich gelangt man durch Einsetzung der letzten beiden Beziehungen in (9b) wieder zu 
den Stromliniengleichungen (8) zurück. 


Es sind nun zwei Fälle zu unterscheiden. erstens der Strahl ohne Zirkulation 


C» (0, h=V, 


! 


wobei 0,== 0,4 1 gesetzt werden kann, und bei dem alle Geschwindigkeiten bis zur Nabe 
(=) endlich bleiben und zweitens der Strahl mit Zirkulation 
EV, 


bei dem die Strömung vor der Nabe aufhören und durch eine innere freie Stromfläche 
begrenzt werden muß, um das Gebiet unendlicher Geschwindigkeiten auszuschließen. 


(3a). Der Fall ohne Zirkulation ergibt nach (8) 


da ,, 

0,° E a ee er er er Tee. 
und nach den beiden Gl. (10) 
C e.V 

E \—1=®, "AT We ee 53 ©} 
Co } 


Schreibt man die Strahlkontraktion so, im Unendlichen vor, so kann man aus der 
ersten der drei Gleichungen die Amplitude a und aus der zweiten Gleichung die Rand- 


geschwindigkeit e„ berechnen. 

Die transzendente Gl]. (1Vec) für z hat unendlich viele Wurzeln. Zu jeder Wurzel gehört 
eine Strahlform mit um so stärker im Querschnitt oszillierender Zusatzgeschwindigzkeit je höher 
die Wurzel liegt. 


(3b). Der zweite Fall des zirkulierenden Strahles e„$0 macht zunächst die Einführung 
folgender abkürzender Bezeichnungen zweckmäßig: 


I (0) = Ja; I (0a) =Ia; Hal); 0, 6. an der äußeren Stromfläche 
ade Kader, Bee o, =; an der inneren Stromfläche. 


Er ergibt je für die äußere und die innere Stromfläche mit Hilfe der Gl. (10) und (8) 
die folgenden Gleichungen: 


ET A 1). 0 
Ada + I Na u . 9a I 


\C 


7a 
f 6,40, 
aJ; bh N; | O0; u Öj.: 
r \Cın a 


Da das erste der Glieder auf der rechten Seite als klein von höherer Ordnung fortfällt., 
erhält man durch Einsetzen von o„ bzw. o; aus Gl. (8) die folgenden homogenen Gleichungen 
für die Konstanten a und b: 





a\J, - Eu B 1 DIN, | ®) Fl. Az 
(ih). 





dA J; [ | 0; R I un. “ | oO; ’ N; =0 











Aus ihnen gewinnt man dann für # die folgende transzendente Bestimmungsgleichung: 





„el el le Yolle leo. un 
Il u | . j ) ww 


3 Me 
20a \Cw 


3 
X 0; 
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. a . Ca . . . 
Ist nach willkürlicher Festsetzung von E und o; (o„ werde von vornherein gleich eins 
ZE 
aneenommen) aus G]. (12) z bestimmt, so liefert eine der Gl. (11) das Verhältnis: 
L_/8,% L /c,„\° 
U f} 4, 
Z, (2) I. z () 
h 7 0a Cr “0; \Cır 
A | C.\ | e.\2 
N, a N, (=) 
2 0a EC 20; \Cu 


Die Einsetzung «dieses Verhältnisses in die Gl. (S) ergibt dann sowohl a als auch eine 





Bedingung zwischen äußerer Strahlkontraktion o„ und innerer Strahlkontraktion o; nämlich: 
| ©. 2 
? | A 
, 2. 0Og Cr 
(dl Ji NV, m v Cw Ög 
} l Balta.ı 
N | n Na 
2 0a Cyı r 
(13) 
| r 2 
I, „(— vr 
j an “Or Cyı 
(dl ; \ i a rr Cu 0; 
\ | f$ IA ‚ 
u KO C yn/ ER 
und daraus: 
j | ( ” | C ö 
(N, I) vw) rl (2) 9a’) 
cd OR , Cor 20a Cu 
' _ (14). 
On ( N ‚ er C, Ta ‚ 
A | Vn ‚| V. # | Rs _ | fu 
da Cu 2. Og Cy 


Man kann also mit Hilfe der willkürlichen Konstante «a nach der ersten der Gl. (13) 
7. B. die äußere Strahlkontraktion o„ vorschreiben und dann aus Gl]. (14) die mit ihr zwang- 
läufie verbundene innere Strahlkontraktion o; berechnen. 

Hat man so aus €, 04, 9. und o; die Größen x, a, b und 6, berechnet, so fordert die 
erste der Gl. (10) für die Randgeschwindiekeiten e„ bzw. e; die Werte: 


») 
Ca ? Ca u % 
| > | | Oa (da I h V) 
Cu En ser - 
(19). 
Ci r;\ Co\ 2 > 4 
| | if O0; (a.J; I \;) 
{ 1 ’ ( 
! u u 





Es muß sieh demnach, wie schon im Absatz 2 gesagt, an der inneren freien Stromfläche 
dieses Potentialstrahles von ringförmigem (Querschnitt eine andere (größere) konstante Ge- 
schwindiekeit e; und damit ein kleinerer Druck p; ausbilden, was offenbar auch aus der 
Zentrifugalkraft der Zirkulation folgt. 

Die wesentliche numerische Arbeit bei der Durchreehnung soleher ringförmiger mit 
/irkulation behafteter Strahlen besteht in der Diskusion der transzendenten Gl. (12) zur 
Bestimmung von z, die an anderer Stelle gegeben werden soll. Auch hier ergibt sich für jede 
Wurzel der transzendenten Gl. (12) je eine Geschwindigkeitsverteilung im Querschnitt, also 
jedesmal eine mögliche Strahlbildung. 

Für den Ansatz (4a) bedeutet dies Überlagerung, so daß die Faktoren z,„ nicht ganze, 
sondern transzendente Zahlen sınd. 


4. Darstellung allgemeiner wirbelbehafteter, freier achsialsymmetrischer und schwach kon- 
trahierter Strahlen. Die hydrodynamischen Grundeichungen für «die allgemeine, achsialsym- 
metrische und auch wirbelbehaftete Strömung lauten bekanntlich in Zylinderkoordinaten 
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von denen die drei ersten die dynamischen Gleichungen, die letzte die Kontinuitätseleichung 
wiedergeben. 
Man setze wiederum: 
u; v=v,tt; wm tt; P=mT+tP - : :» :.: .: 0.AM, 
wo nach den Gl. (1) », und w, beliebige Funktionen von r und 
v, bzw. mw; P <p 


ee 


l 1 1 


sind und p, durch Gl. (la) angegeben wird. 


Die Gl. (16) lassen sieh dann umsehreiben in: 





Io» dm, om, 
- z m —- M, 
u OL 0oL ’ A 7) 
Iop, Ö, | 
'ZR rn 
u 0 oO 0 A Ü v0 lc 
(Iba). 
or, 1 0 0 
m, Ay ",o Yo 1% 0) 
#) 1,0) = 
s 00 e- OL 


Es werde nun wieder wie im Abschnitt 2, um sich der gesuchten Strömungsform des 
schwach kontrahierten, im Unendlichen zylindrisch werdenden Strahles anzupassen und 
entsprechend dem in z linearen und in den Koeffizienten von z unabhängigen Charakter der 
(Gl. (16) und (16a) für die Zusatzgeschwindigkeiten der Ansatz gemacht: 

u, =U(o) e-*» 
v, V(o)e *- 
w, == Wl(o)e-*-» 
Pı Pio)e ” 


Die Gl. (16a) werden dann: 





U 
| 

P +zw,l! +20,V o"' 
u 
| (6b). 

A E d 

x” W, 0 do” Oo) 
w; ; 
B—g"" (lo) Mu 
- do 


Setzt man die beiden letzten Gleichungen in die beiden ersten ein, so erhält man: 





x d = 
P=#,90° (lo)  w,l 
“ = do ” 
(I6e). 
| p' N i . 4 Vor: d Ä 
“Ww ya br 0) [®) ",o) 
u ur u zm, . Ye 
Die Elimination von P ergibt dann: 
| +.) v4 
"oO d 
‚’ rr Y ‚ ' N “> 2) — 0 -_ 
Ih N ! N VO ı1- N IR U, 0 i. w.IoO “ HH” - ",o) 0) r ; (Ir). 
Ü & 0 CC r 0 ı\G m, do | or | 


Jede wıllkürliche Annahme über die Grundgeschwindigkeiten w, und », erlaubt also, 
aus den Gl. (16b) und (17) die zugehörigen Nachbarlösungen der U’, V, W, P’ zu finden. 


3. Fall der Potentialströmung als Ausgangslösung. Wir fragen zunächst, welche Nachbar- 
lösungen sich ergeben, wenn als Grundlösung eine Potentialströmung, also: 


Cw CyFu J 
=. a . (2a) 
N u 


[I 


nach Gl. (2a) angesetzt wird. 
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Ks ergibt sich dann: 


/ Io (x 0°) () 
also / Ad J (20) f, \,(z0) 
1% 0) 
| d 
| lt o.J, (zo) bo N,(zo)|, 
z0do ’ ” ” “ 


wo ./, und N, wieder die Besselschen Funktionen erster Ordnung bedeuten. 


Nach einer bekannten Differentialformel '’) ıst nun. wenn wie üblich aJ, + N, =Z 


1 I l 
vesetzi wird 


d 
|oZ (zo)| 20 /,(z2 0). 
do 5 z ” 
demnach W= 2,120) = a J, (#20) +-b N,(z0). 


Ks führt also die Potentialströmunge als Grundlösung auch für die Nachbarlösungen 
wieder auf Potentialströmungen, wie sie im vorhergehenden Abschnitt behandelt waren. 


6. Fall der nach einem Potenzgesetz wachsenden achsialen Grundgeschwindigkeit (1%,). 
Für eine zweite Fragestellung fordere man eine Zunahme der achsialen Grundgeschwindig- 
keit »=, nach einem Potenzgesetz und behalte die Voraussetzung, daß die tangentiale Grund- 
eeschwindiekeit dureh eine Zirkulation gegeben sei (oder null sei), bei. 


Man setze z. B. zunächst: 


die Gl. (17) wird dann: 


U+-Ur'+U#=b0, 


also: ! Ad J(zol)—b \„(#0) /,|2 0) 
V 0 
| EEE 
| (0 7,1z0)). 
z0do 2 = 


Da diese Besselschen Funktionen auf der Achse (0-0) durchweg unendlich werden, 
so ist bei diesem Ansatz nur ein Strahl mit rineförmigem Querschnitt (ein Strahlschlauch) 
zu konstruieren. 


Dareeen führt der Ansatz: 





I ie, WE an ee 
zu: Hr —- ee -9I=%, 
also zu: 
| Zn 1180) 
| u 
(18a). 
8 , 
Il 07, (zoll 
0 do he: ie 
Wählt man » 2: (ie, 60°), so ergibt sich: 


u /,\zo) 
W /,\20). 
Dies ist dieselbe Lösung wie für die Potentialgrundströmung (m, = e./r a). 


Ks führt also der Ansatz quadratisch wachsender Achsialgeschwindigkeit », zu denselben 
Strahlformen und zu derselben Verteilung der Radialgeschwindigkeit wie der Ansatz 
konstanter Achsialgeschwindigkeit. 


Für » > erreben sieh dageren andere Strahlformen. 


10, Sjelhn ‚B. Jahnke-Emde. Funktionentafeln, Tenbner 1909, S. 16». 
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7. Fall der nach einem Potenzgesetz wachsenden radialen Geschwindigkeit (1u,). Zu be- 
merkenswerten Ergebnissen führt ferner der ebenfalls physikalisch vernünftige Ansatz einer 
nach einer Potenz von o wachsenden Radialgeschwindigkeit, also der Ansatz: 


"7,5 Sn TE Ca  ;  ° 


(Gl. (17) nimmt dann die foleende Form an: 


w. —w,0" w,[xz” + on?’ —1)|=0 
Wu sr 
Aus Gl. (16b) folgt: 
N a 
x 


8. Strahl mit vorgeschriebener achsialer Wirbelverteilung der Ausgangsströmung. Als 
weiteres Beispiel möge eine Strömungsform gewählt werden, bei der in der Grund- oder 
Ausgangsströmung Wirbellinien in achsialer Richtung vorhanden sind, jedoch so, daß die 
Tangentialgeschwindigkeit in der Achse null ist. 


Man setze entsprechend: 





Us 
v, - o" 
Nu 
(20). 
C in 
= konst 
(d 
(+l. (17) verwandelt sich dann in: 
rmx rs | Y | 07) 15) ) C» £ 157 l 
U’+-lUo u 10° “ ä in Hi)o®’%n \ 0, 
- | fe B - 
| u 


Sa. Die Annahme n» = 0 ergibt eine bequeme Lösung, bei der man allerdings mit Hilfe 
der Besselschen Funktionen zweiter Art den Nabenbereich ausschließen muß, nämlich die 
foleende Gleichung: 





| 1} u) 
is C»\Y" 
U” +-U'o | 1=° O0 11 2 "| | 0 
ö | | Cw/ || 
.. 12} > Ev ‚ ® 
also wenn man zur Abkürzung p? | >| | schreibt: 
u 
= Zy(x0) 
1 © ıy 
—— N iniX 0) 
A Ci” r ” . N . ; ; . (20a). 
‚3 6 
1% d 0 Zn (2 o) 
“oO d 





Sb. Der Wert » 1 entspricht der Drehung des Strahles als starrer Körper mit der 


a ° ‘ . CO» . 5 
Winkelgeschwindiekeit ‚, und liefert: 


i . .; 2, 
woraus, wenn man zur Abkürzung #4,’ »#° + |. -) setzt 





U = 2, (x, 0) 
ee 

7 Mac. iR (206). 
warst | Z, (x, 0) 

jr 0 do [o 1, (%, 0) 0%, A, |\X%,0 


Da die Geschwindigkeiten U und V an der Achse null werden und W dort endlich 
bleibt, wenn man in Z, und Z, nur die Besselschen Funktionen erster Art wählt, so kann 
mit diesen Funktionen ein voller rotierender Strahl gebildet werden. 
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9, Durchführung der Rechnung für einen starr rotierenden Strahl der Ausgangslösung (8b). 
Als Beispiel für die weitere Berechnung der freien Strahlform, d.h. für die Berechnung der- 
jenigen Parameter und Integrationskonstanten, für die es eine freie Strömungsfläche gibt, soll 
der zuletzt angesetzte Fall der Strahldrehung als starrer Körper in der Grundlösung gewählt 
werden. 


Die Gleichung des Meridianschnittes einer Strömungsfläche ist wie in Gl. (5) 


do HM, Ya r 
e Ue", 
sc m 


oder wenn man in Gl. (20b) beachtet, daß nur die Besselschen Funktionen erster Art 
venommen zu werden brauchen: 


do Pi 
j aJ,(x,0)e",*. 
- Fa . 
Setzt man wieder: 

Oo | 6. Wwoo » 


so ereibt die Integration und die Ausrechnung bis zu kleinen Größen erster Ordnung: 


r,d a \ do 0 
en "J,[z,(1-60)] J (z,) 


Die Bedingung, es soll 6 =0 für DO werden, ergibt wieder: 


rd 


% Cw 


Die Kontraktion im Unendliehen wird ähnlich wie in Gl. (Sa) 


i ar, 
o AJ,(,) Ze) ee A VEN 


Die Gleichung einer Stromfläche ist also ähnlich wie in der dritten der Gl. (8): 


Ar, Re ' 
9 — As)il-e"]=0,(1-e") . 0.0.0. . (200). 


l 
m 


Die Bedingung, daß diese Stromfläche freie Stromfläche sein soll, erhält man nun wie 
in Gl. (9) dadurch, daß man auf der freien Stromfläche die Bedingung konstanter Ge- 
schwindiekeit stellt, und zwar: 


I)mr,w e 


"2 -L2r,n.— w A = 


0 
Die Einsetzung der Lösungsfunktion der Gl. (20b) liefert dann: 


. 2 „ |») » 2 2 
Cy J (z,) . ar ed. Cw “ Cr Co R 
AN ed 2 m u 


ae >| 
f 
Ka Ya Ya 


er x 
Entwickelt man bis zu den Gliedern erster Ordnung in o= | -o und beachtet Gl. (20«) 
und (204), so ergibt sich: 


200 Mey d (7,) 
ert20 | 
( 


! 1 | > rn > > 
2, J/,(2,)]| Ca Ya tt ey (1 20). 


2 I. (x,) \Cw) 4 


Damit eine der Stromflächen freie Stromfläche konstanter Geschwindigkeit (e,) Ist, 


müssen also die folgenden beiden Bedingungen erfüllt sein: 


A ET 20,) RE | De BE 0 2: 


Schreibt man nun die Grundgrößen e,„ und ce, der Geschwindigkeiten und die Kon- 


traktion o, vor, so bedeutet die letzte Gleichung eine (transzendente) Bestimmungsgleichung 


für x und x, +4”) die für jede Wurzel z, (z,) einen Strahl mit jedesmal anderer 


( u 


Geschwindigkeitsverteilung liefert, der übrigens den anderen Strahlen in gewisser Weise über- 
lagert werden kann. 


Ist x und x, bestimmt, so liefert Gl. (20«e) die Konstante «a. 


Damit ist dann, wenn man noch den Außenradius vr, im Anfangsquerschnitt vorschreibt, 
die Gesehwindiekeitsverteilung im ganzen Strahl gegeben, nämlich nach den Gl. (20b) dureh: 


















































M- 


di 


dr 
Vet 


be 


dr 
A 


m 


in 


U 


N 
d 








Band 12, Heft I 





Februar 1939 Reißner, Achsialsymmetrische, freie Flüssigkeitsstrahlen 35 
> X Cw 
" m e—-"So, J.(xz,0) 
1 Pr a 1 ‚. 
‚ 7 J, (#,) 
- Cy 
v—®v e-*%:0_D J(z,o) 20f) 
bee zur ee eu ee 
’ Tr dt 
Kl 
ww, -—e"0 0%,J,(x2,0) 





® Ya J; (x) 
Schließlich auch die Druckverteilung 
p=p,-umw, W, 


die sich durch Einsetzung aus Gl. (la) und der letzten der Gl. (20b) umformt in: 


„> 


.(d-0) +; 2 (2,0) e * (20 8) 
P’ sn Ia ha . . 2 0° Tr u . 2 a Ö 0 K, 8 : ” r R ä n , z _. = v 
ra Nu J,(z,) " 


Im Unendlichen (© =») geht, wie es sein muß, die Druckverteilung in die des zylin- 
drischen Strahles über, während im Anfangsquerschnitt (©=0) das letzte Glied die Überdruck- 
verteilung angibt, die sich bei der Kontraktion des Strahles in Geschwindigkeit umsetzt. 

In dieser Weise kann die Rechnung auch für die andern nur bis zur Konstanten- 
bestimmung durchgeführten Beispiele weitergeführt werden. 


10. Zusammenfassung. Es wird die allgemeinste zylindrische Strömung mit freier zylin- 
drischer Oberfläche sowohl für den wirbelfreien als auch für den wirbelbehafteten Fall als 
Ausgangszustand zugrunde gelegt. 

Durch partikuläre Integrale der Potentialgleichung bzw. der allgemeinen hydrodyna- 
mischen Gleichungen werden Nachbarlösungen aufgestellt, die im Unendlichen asymptotisch 
in die zylindrische Lösung übergehen und es wird ihre sich schwach kontrahierende (unter 
Umständen auch divergierende) freie Mantelfläche bestimmt. 

In den Fällen, in denen die Zirkulation bzw. Rotation des Strahles unendliche Ge- 
schwindigkeiten an der Nabe bedingen würde, wird die Geschwindigkeitsverteilung, die 
Druckverteilung, und die äußere und innere freie Mantelfläche der zugehörigen Ring- 
strahlen bestimmt. 

Wenn es nur darauf ankommt, Beispiele von solchen Strahlen zu konstruieren, genügt die 
Ansetzung eines eingliedrigen partikulären Integrals mit einer Besselschen Funktion nullter 
Ordnung, in der der Argumentfaktor einer transzendenten Gleichung zu genügen hat. Mit diesem 
partikulären Integral ist dann für jede Wurzel dieser Gleichung je eine Geschwindigkeits- und 
Druckverteilung insbesondere im Anfangsquerschnitt verknüpft, welche man sieh durch einen 
Beschleunigungsmeehanismus (einen Propeller) im Anfangsquerschnitt erzeugt denken kann. 

Wenn es darauf ankommt, eine vorgeschriebene Verteilung der Geschwindigkeit und 
des Druckes im Anfangsquerschnitt wiederzugeben, ist dies durch Überlagerung beliebig 
vieler partikulärer Integrale zu erreichen. 

Die partikulären Integrale enthalten bei vollem Strahlquerschnitt nur die Besselschen 
Funktionen erster Art, bei Ringquerschnitt sind es Binome aus den Besselschen Funktionen 
erster und zweiter Art. 

Für einen wirbelbehafteten Strahl kann man auch bei eingliedriger Lösung der Geschwin- 
digkeitsverteilung oder der Druckverteilung oder der Wirbelverteilung im Querschnitt eine 
willkürliche Bedingung auferlegen. In allen durchgerechneten Fällen ist die Form der 
Mantelfläche bestimmt durch die folgende Beziehung zwischen spezifischer Kontraktion o in 
der achsialen Entfernung © und der spezifischen Kontraktion o, für [= «: 

oO O5 (1 e *>), 
wo x der aus einer transzendenten Gleichung zu bestimmende Argumentfaktor einer Bessel- 
schen Funktion ist. 

Von den 10 Bestimmungestücken bzw. Funktionen des Strahles, nämlich den beiden 
Komponentenverteilungen w,(r) und ®,(r) der Grundgeschwindigkeit und den drei Kom- 
ponentenverteilungen U, V, W der Zusatzgeschwindigkeiten, ferner den spezifischen Strahl- 
kontraktionen (im Unendlichen) innen und außen o,,a und o,,; und den Begrenzungsradien 
r„ und r; im Anfangsquerschnitt und der Druckverteilung können fünf willkürlich vorge- 
schrieben und die übrigen aus ihnen berechnet werden, wobei es immer erreicht werden 
kann, daß die Voraussetzung der Kleinheit der Zusatzgeschwindigkeiten und der Kontraktionen 
genügend gewahrt bleibt. 

Es erscheint möglich, von einer achsialunsymmetrischen, aber zylindrischen wirbelfreien 
oder wirbelbehafteten Strahlströmung ausgehend, auch freie achsialunsymmetrische, sich 
schwach kontrahierende Strahlen von vorgeschriebenen Eigenschaften aufzubauen. 218 
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Über Einflußlinien und Einflußdiagramme. 
Von Emil Miller, Wien. 


err P. Nemenvi hat in seiner bedeutsamen Arbeit „Über die Singularitäten der Elasti- 
H zitätstheorie“ ') unter anderm die Bedeutung der Differentialquotienten von Einflußlinien 
infolge wandernder Einzellasten für den Fall gerader Stäbe mit konstantem Trägheitsmoment 
dargelegt. Im folgenden mögen die entsprechenden Betrachtungen für gekrümmte Stäbe mit 
überdies veränderlichem Trägheitsmoment angestellt und im Anschluß daran einige Be- 
ziehuneen zwischen den Einflußlinien für wandernde Einzellasten verschiedener Richtung 
besprochen werden. 


Die Ableitungen der Einflußlinien für eine Einzellast als Einflußlinien für besondere An- 
griffe. Kine lotreeht gedachte Einzellast P=1 wandere über einen beliebigen Stabzug und 
bewirke, wenn sie im Abstand a von einer gewählten, gleichfalls lotrechten Achse A angreift, 
daß eine VEWISSE (röße ıdD den Wert 7 annehme: dann Ist 7 —= f(a) die Gleichung der Kinfluß- 


.. | 
linie von P für eine Einzellast (Abb. 1). Läßt man an der Stelle a P ‚an der Stelle 
re 
a--e P angreifen, so bewirkt der Angriff dieses Kräftepaares ein 
e 
| | 
D (ae la 
—f — | 


Dieser Wert geht beim Grenzübergang e>V0 in P=f’(a) über, während gleichzeitig die 
Kräfte immer näher aneinander rücken. Da hierbei jedoch das Moment des Kräftepaares 
unverändert eleieh 1 bleibt, ist die Belastung im Grenzfall mit dem Angriff eines Momentes 
\/=1 an der Stelle a zleiehbedeutend. Es folgt demnach der 

Satz 1: »„’ gibt die Einflußlinie von 2 für ein über den Stabzug wanderndes Moment 
M==1. das im Sinne des Uhrzeigers dreht, an. 





Abb. 1. Abb. 2. 

Dieser meines Wissens zuerst von Herrn Schwätzer?) ausgesprochene und hier nur 
der Vollständigkeit wegen nochmals ausgeführte Satz gilt wie ersichtlich unabhängig von der 
Form des Stabzuges und dem Verlauf seiner Trägheitsmomente. 


läßt man nun bei a ein Moment M = ‚ bei a+e ein Moment JM: angreifen 
e rt’ 
Abb. 2), so liefert dieses „Momentenpaar* („Doppelmoment* nach Nemenvi) ein 
es gr 
P= f (a+e) f (a 
e e 


während der Winkel, den die Tangenten an die Biegelinie des Stabzuges an den Stellen a und 
« --e miteinander einschließen, 

l eseeg 

e EJ 

wird. Hierin bedeutet g die Neigung einer der Tangenten an den Stabzug im Intervall e 
und .J das Trärheitsmoment an einer Stelle des Intervalls.. Für e>0U wird nun einerseits 


«ıb . f” A 
anderseits 
sec 
= — ii). 
F M 


1), Zeitschr. f. angew. Math. u. Mech., Bd. 1, Nr. 4, S. 383. 
?\ Schwätzer, Statische Berechnung der eingespannten Bogenträger bei beliebigen Lastangriffe. Österr. 
Woechensehr. f. d. öffentl. Baudienst, Jahrg. 1912, S. 184. 



































Bee 





jand 12, Heft 1 Z ” er A 22 
u Bi 1932 Müller, Über Einflußlinien und Eintlußdiagramme 37 





worin i (a) eine nur von der Lage und den Abmessungen des Stabzuges abhängize Funktion 
ist und p und J sich jetzt auf die Stelle a beziehen. Denkt man sich also bei a in den Stab- 
zug ein Gelenk eingeschaltet und die dort zusammenstoßenden Stabteile um den Winkel z=i 
derart gegeneinander verdreht, daß der Stabzug beim Fortschreiten im Sinne wachsender a 
einen Kniek nach rechts erfährt, so bewirkt dieser Vorgang, den wir kurz einen „Knick x“ 
nennen wollen, ein 

Da f” (a), 


Demnach besteht der 


Satz 2: 7” (a) gibt die Einflußlinie von P für einen wandernden Kniek von der im 
allgemeinen veränderlichen Größe i an. 

Für J=0 (Gelenkstelle) wird <=, also auch der anzubringende Kniek und damit im 
allgemeinen „’=®x; d.h. 7’ macht an einer solchen Stelle einen Sprung. Um dessen Größe 
zu finden, ist an der kritischen Stelle ein Momentenpaar M„:.,=1.M,_. | anzunehmen. 


) 








Seine Wirkung auf ® ist dann durch f" (a +0) - f’(a - 0) gegeben. 
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Abb. 3. - \bb. 4. 
i ' s i FR (a) 
In Fortsetzung unserer Betrachtungen nehmen wir nun wieder beı a einen Knick 
ev 
r : IR: N (Ad - e) r E ) z ’ 
bei a--e einen Knick 5 an (Abb. 3). Nach dem vorhergehenden Satze ruft dieses 
„Kniekpaar“ 
Z} [Z, 
l (dt 0) ) (Ad) 
D- / | / 
e e 
’ u (a) . 
hervor. Anderseits bewirkt der kleine bei m angebrachte Knick eine Verschiebung 
Pe 


(a) ’ Br 
des Punktes » senkrecht zur Sehne mn um eseep=tl(a)seep und einen Knick des 
e 


Stabteiles II gegen Stabteıl I von 


t (a) (ae) 
w— -Ap+ 


e @ 


Bei dem Grenzübergang e>0 wird P= 17" (a), die Verschiebung fällt in die Richtung der 
Stabzugnormale in m, und der Knick wird, da Ig mit e verschwindet, 


.r 
H r (Aa) 


Die Größe der Verschiebung werde durch ihre lotrechte Komponente » gekennzeichnet. Man 
erhält 


5 (a): 


ein positives » bedeutet eine Verschiebung nach links beim Weiterschreiten in der Richtung 
der wachsenden a. Das Anbringen eines Knickes z und einer Verschiebung senkrecht zur 
Stabtangente. deren lotrechte Komponente » ist, wollen wir das Einfügen eines „Knieksprunges 
nennen. Wir können dann den Satz aussprechen: 


., 


(x.v) 


Satz 3: 7” (a) gibt die Einflußlinie von P für einen wandernden Knicksprung von der 
im allgemeinen veränderlichen Größe (, ü) an. 
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Wir bemerken, daß (#2, 0) das Anbringen eines Knicks ohne Verschiebung, (0o,v) das An- 
brinzen einer Verschiebune ohne Kniek, d. h. einer derartigen Verschiebung bedeuten, daß 
die Tangenten an die beiden Stabteile in a zueinander parallel bleiben. Ferner ist unmittel- 
bar klar, daß 


und alleemeın 


%,, V, %, 0,)= (8%, + %., td, +; 


.y . .r . y 
ta), ice] i (a) (a) 
(Am; ") e ( e | 
und in @a—+e einen Knicksprung 
(a-+te),iia-+e dia-+te) i(a-+te)| 
(Ans Un) | i | 
P P E 


an, so wird hierdurch eın 


«db / T / 


hervoreerufen. Der hierbei entstehende Knick des Stabteiles II gegen den Stabteil I beträgt 
Abb. 4 


Ä ii) Ä i (ae 

x u + Io + x. = -Aor + 

n / ) r / r 

und die Verschiebung senkrecht zur Sehne mn... Um Sec 9 — AZmesecpg—+r„secgy, wobei 


der Einfluß von Ig, der ohnehin mit e verschwindet, von vornherein vernachlässigt ist. Für 
e>V wird nun 


D 7 IV) % ta 
und 
(a) 2 t(a-te) f 
v = lim (u, — %m e+v,) =lim — + (+ 24V (a) 
> (l e—!|l f € 


Demnach gılt der 


Satz 4: Via) gibt die Einflußlinie von P für einen wandernden Knicksprung von der 
) ;t 


im allgemeinen veränderlichen Größe (i”, 27) an. 
Durch sinngemäße Wiederholung der letzten Betrachtung wird leicht gefunden der 


Satz 5: 7)’ (a) gibt die Eintlußlinie von ® für einen wandernden, im allgemeinen ver- 
änderliehen Knieksprung 


%,.v0) DL > m — din 3]. 
Für eerade Stabzugteile konstanten Trärheitsmoments verschwinden laut den Gäisen 4 


und 5 alle höheren Differentialquotienten als »’”. Daher bestehen die Einflußlinien längs 


solcher Stäbe wıe bekannt aus Parabeln 3. Ordnung. 

Beispiele. Die Richtigkeit der abgeleiteten Beziehungen möge noch an zwei Beispielen 
bestätigt werden: 

I. Einflußlinie der Durehbiegung in der Mitte eines freiaufliegenden 
Balkens mit veränderlichem Trägheitsmoment. 


Sie ist eleich der Biegungslinie für die in Abb. 5 dargestellte Belastung. Die Differential- 
eleichung des linken Astes dieser Biegelinie lautet 


und demnach 
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Die Anbringung eines Knickes 1 bei a gibt, wie Abb. 6 zeigt, eine Durchbiegunge der Träger- 
mitte um 


2) 3 e Zul 2 


Die Anbrinzunz einer bloßen Verschiebune "=1 nach Abb. 7 











| 
rn) : 
I» / BZ 
| N 
. 1 | =7 en 
on Mar re u _ | U d ER 
7 Tl L-a D 
_ ( | a er 2 
m 5 o nr) 
(A —— 
I 5, _ Q a l-a - = a | 
2 2 2 2 
Ar RAı62 Ze AA162 7 
Abh. ». Abh. 6. Abb. 7: 


Die negativen Vorzeichen ergeben sich, da der positive Knick und die positive Verschiebung 
eine Hebung der Trägermitte hervorrufen. Demnach entsteht 


. : £ Ad . [zZ 
infolge (x,0)=(,0)... 0 a | 
( ”) j’ ) . : ae, 
H ı,ı ( 5 ı > Y) 
i dd 
x.n) ” 2) i N f 
») 
.) 00 
Ad rt 
/. h Ta .„ did V) 
z.v0) >=(1 r 31 a ) .) ! .) 1] . 


2. Einflußlinie des Horizontalschubes eines halbkreisförmigen Zwei- 
eelenkbogens mit unveränderlicehem Trägheitsmoment (Abb. 8). 
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Abh. 8. Abh. 9. 
Mit bekannten Bezeichnungen wird 
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Die Anbringung eines Knickes x dag (Abb. 9 


„ 
Ns 07 da yda 


demnach die Anbrinzung eines Knickes x 1 eın 9; 


rıbt eın 


er 


J 


y; folglich wird 


ÖOy; 2E.J 
H, n „.c084 
OH on 
Die Anbringung einer Verschiebung » - 1 (Abb. 10) bewirkt ein Op, tgq 
2) » 2 E43 
H,=z- tg 
OH vn 
. . Sec 
Infolgeex ; wJ entsteht demnach 
) 
Hı-;=4#H, z I" 
a 





\bb. iv. 


Zur Ermittlung von //”” wird ö gebraucht. Man findet aus 


dag see q 
rda COSgQ da a 
da r 
und damit oda di dg tgysec’g 
! . - 
dc da da A FE J 


Infolge x =?’ wird nun 


AR Ai Thai 
en 
und ınfolge > 
u—Hh Bi m 
also infolge e,%). I __ H,+iNH, BE, 


Zusammenhang zwischen Einflußlinien für lotrechte und waagrechte Belastung. 


Ztschr. f. angew. 
Math. und Meeh. 





somit eın 


Als Anhang 


zu den vorausgehenden Betrachtungen möge aus Satz 1 noch ein Zusammenhang zwischen den 
Kıntlußlinien für eine lotrechte und waagrecehte Einzellast behandelt werden. 


Es sei wiederum », (a) die Einflußlinie einer Größe ® für eine über den Stabzug wandernde 


lotrechte Einzellast VI und Sta) die entsprechende Einflußlinie für eine waagrechte Last 
I" 1. Dann ist der Einfluß einer unter dem Winkel 9 gegen die Lotrechte wirkenden Last 


PP’ 1 Abb. Il) gezeben dureh 


«D |: cos + -| - sın q 
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Ordnet man nun jedem Punkt a einer waagrechten Geraden einen Punkt m zu, indem man von 
a um (a) aufwärts und um &(a) nach links geht, so ist die Wirkung von P=1 durch die 
Strecke a m’ gegeben; hierin ist m’ die Projektion von m auf die durch a parallel P ge- 
zogene Gerade. Dreht sich die Kraft P um ihren Angriffspunkt, so bewegt sich m’ auf dem 
Kreis. dessen Durchmesser a m ist. Die Länge der durch «a parallel P gezogenen Kreissehne 
oejbt dann den Einfluß von P=1 auf ®P an. Bei der gewählten Anordnung ist der Einfluß 
positiv, wenn die Pfeilspitze innerhalb, negativ, wenn sie außerhalb des Kreises zu liegen 
kommt: das letztere ist bei der Kraft in Abb. 11 der Fall; dabei ist vorausgesetzt, daß die 
Kraft immer so eingetragen wird, daß die Pfeilspitze nach «a fällt. Wirkt P in a parallel a m, 
so wird ® ein Maximum; wirkt P senkrecht a m, so verschwindet ®. 


In Abb. 12 sind zu einer Reihe von Ängriffspunkten a die zugehörigen m dargestellt; 
‚lie Zuordnung der a und m ist durch die Verbindungsstrecken a m ersichtlich gemacht. Das 
Einzeichnen der Kurve der m ist für das Interpolieren vorteilhaft. Schließlich sind noch die 
oben erwähnten Kreise eingetragen. Dieses „Einflußdiagramm* gibt über den Einfluß jeder 
beliebigen in irgendeiner Richtung an dem Stabzug angreifenden Kraft erschöpfend Aufschluß. 
Die Strecken am mögen „Gräten“, die Winkel «a, die sie mit der Lotrechten einschließen, 
„Grätenwinkel“ genannt werden. 





Abb: 9. Abh. 13. 


Es ist nun bemerkenswert, daß, wenn der Stabzug und etwa die 7-Linie gegeben vor- 
liegen, hiermit auch die &Linie im wesentlichen bestimmt ist; ja daß allgemein immer dureh 
zwei der drei Kurven (Stabzug, »)- und &-Linie) die dritte gegeben ist. Trägt man nämlich 
(Abb. 13) die & als Funktion von „, die „, als Funktion von x auf, so geben nach Satz 1 so- 


- 


än_. dE 


wohl j. wie die Einflußlinie für ein über den Stabzug wanderndes Moment M 1 an. 
X (dd 


D. h. es muß sein 
dı ds de 
! ee ee ee . ib 


dr dy ydı 
oder de yydı, 


worin die Striche Differentiation nach .r bedeuten. Durch Integration folgt hieraus 
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Die Beziehungen an geraden Stäben. Beispiel. Handelt es sich um einen geraden Stah- 
zugteil, so ist y’  const und es wird 


=yn+t 


2 
oder wem a. .. u. 0% 4 ee 
7 7 


Ir 


(sehört überdies zu & 0... VO und wird vorausgesetzt y' $# x. so wird © 0. Man findet 


- 


14 


dann tga=y; 


d. h.: gibt es in einem geraden Stab eine Stelle, an der & und „ gleichzeitig verschwinden, 
so stehen die Gräten des zu dem Stab gehörigen Diagramms auf dem Stab senkrecht. 


Gibt es längs eines Stabes keine gemeinsame Nullstelle für & und »,, so ist nach (3) 


ee Caygy. 
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liegt also beispielsweise die „-Kurve für einen geraden Stab gezeichnet vor (Abb. 14), so 
schneiden die dureh die Kurvenpunkte senkrecht zum Stab gezogenen Geraden auf der Grund- 
linie die Streeken E— (' ab. Trägt man diese von einer Waagrechten e nach aufwärts auf, 
so erhält man hierdurch bereits die Form der &-Kurve. Die Abstände dieser Kurvenpunkte 
von einer um € tiefer als e gelegenen Waagrechten d stellen die £ dar. Die £ werden natür- 
lich auch dureh eine Parallelverschiebung der (£— O)-Kurve um € nach aufwärts erhalten. 


Da nun nieht mehr konstant, sind die Gräten in diesem Falle auch nieht mehr zueinander 
/ 


parallel. 


Die Einflußlinie des waagrechten Auflagerwiderstandes H am rechten Gelenk des in 
Abb. 15 dargestellten Zweigelenkrahmens für eine über den Rahmen wandernde lotrechte Last 
> 


besteht aus zwei Parabeln 3. Ordnung. Am linken Gelenk ist &= „7° 0; demnach sind die &, 


a In | h | | 
nach (3) längs des Stabes I proportional den 77; und zwar ist, datga= 5, &° „9. Im Stab I 


dagegen fällt die Nullstelle von 77 (am rechten Gelenk) nicht mit einer Nullstelle von &7 zu- 
sammen; denn am rechten Gelenk wird Sp 1. Folglich schneiden die durch die Punkte 
der »-Kurve senkrecht zum Stab gezogenen Strahlen auf der Grundlinie Strecken &y7- (ab, 
die ın diesem Falle, weil sie links der Fußpunkte der Ordinaten zu liegen kommen, negativ 
sind. Sie werden daher in Abb. 15 von einer Waagrechten nach abwärts (gestrichelte Linie) 
aufgetragen. Hiermit ist die Form der &pKurve gefunden; sie muß aber jetzt noch parallel 
zu sich selbst so weit verschoben werden, bis am rechten Gelenk die Ordinate 1 auftritt. 





\bb. 1. Abb. 15. 


Die Beziehungen an krummen Stäben. Beispiel. Ist der betrachtete Stab gekrümmt, so 
muß auf Gl. (2) zurückgegriffen werden. Aus dieser sieht man zunächst, daß man wegen der 
Vertauschbarkeit von y und » dieselbe &Linie erhält, gleichgültig, ob man die y-Linie als 
Stabachse und die »-Linie als Einflußlinie für lotrechte Einzellasten oder umgekehrt die y-Linie 
als Eintlußlinie für lotrechte Lasten und die »-Linie als Stabachse auffaßt. Obwohl der 
natürlichste Weg, die &Linie aus der Stabachse und der n-Linie zu konstruieren, schon in 
Abh. 13 angedeutet ist, sei doch noch eine für die Ausführung vielleicht etwas bequemere 
Konstruktion angegeben. Sie ist gleichfalls aus Gl. (1) unmittelbar abzulesen. wenn man diese 
als Differenzengleiehung in der Form 


schreibt. 


In Abb. 16 ist das Intervallöi kin eine Anzahl Unterabschnitte I. geteilt. Zu jedem 
Ir gehört ein Iy, z.B. zu I, ein I, rp. Zieht man nun durch p die Parallele pg zu 
n,. wobei », die mittlere Normalenriehtung der »-Linie im Abschnitt 1.r, bedeutet, also etwa 
die Normale zur Sehne des Abschnittes 2, so schneidet p g auf der Waagrechten durch r die 
Streeke rg I, ab. Mit Hilfe der so gefundenen IF erhält man die Form der &-Kurve, die 
man sowohl als Funktion der y wie in Abb. 16 oder als Funktion der x auftragen kann. 
Zur Bestimmung der & selbst ist im allgemeinen wie vor noch eine Parallelverschiebung der 
vefundenen Kurve oder die Angabe der Grundlinie notwendig, von der aus die &-Werte zu 
messen sind. Zu diesem Zwecke muß wieder Ü aus zwei zusammengehörigen Werten & und 
ermittelt werden. 


In ähnlicher Weise läßt sich jede der drei Kurven aus den zwei andern finden. 


rr 
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Als Beispiel für die Anwendung der erörterten Beziehungen sei noch die Einflußlinie 
des Momentes im Punkt m eines frei aufliegenden, kreisförmigen Bogenträgers für waag- 
rechte Lasten aus jener für lotrechte Lasten und aus der Stabachse abgeleitet (Abb. 17). Zu 
beachten ist hierbei, daß die Vorzeichen der 1£ sowohl von y’ wie von Ay abhängen. A&, 
bis 1, sind danach positiv, da von 1 bis 3 die 7 und Ay positiv sind; JE, und A£, sind negativ, 
da die zugehörigen n’ negativ, die Ay aber positiv sind, und IE, bis A&,, sind wieder positiv, 
da hier 7’ und Ay negativ sind. In Abb. 17 sind die zu einem Punkt p gehörigen 15, durch 
die Strecken bestimmt worden, die die durch p gezogenen Normalen zur n-Linie auf der 
Waagrechten durch p — 1 abschneiden. Bei diesem Vorgang sind die von der Lotrechten 
durch p nach links fallenden Abschnitte A& positiv, die nach rechts fallenden negativ in 
Rechnung zu stellen. Eine Parallelverschiebung der so erhaltenen £-Linie ist hier nicht mehr 
erforderlich, da sie bereits so aufgetragen wurde, daß am linken Auflager, wie es sein muß, 
E=( ist. 


Trägt man » als Abszisse und das zugehörige & als Ordinate in einem rechtwinkligen 
Koordinatensystem auf, so erhält man eine Kurve K (Abb. 17), deren Punkte je einem be- 
stimmten Punkt der Stabachse, nämlich jenem, zu dem das betreffende »; und & gehören, zu- 
geordnet sind. K hat nun die Eigenschaft, daß ihre Tangenten parallel den Tangenten in 
den entsprechenden Punkten der Stabachse sind. Denn der Riehtungskoeffizient einer Tan- 





de 

u: dx i ’ , 

gente von K ist rg und dieser Wert ıst nach (1) gleich . 
/ 
dr 
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Die Darstellung der entsprechenden Beziehungen bei räumlichem Kraftangrifi. Zum 
Schlusse sei noch darauf hingewiesen, daß sich dieselben Betrachtungen, die bier über 
Einflußdiagramme von ebenen Traggebilden angestellt wurden, ohne Schwierigkeit auf 
räumliche Gebilde übertragen lassen. Anstatt der Einflußlinien treten dort Einflußflächen 
auf. So läßt sich z. B. bei einer Kuppel jeder Größe P je eine Einflußfläche für Lasten in 
der x-, y- und z-Richtung eines rechtwinkligen Koordinatensystems zuordnen. Um die Wirkung 
einer beliebigen an der Kuppel angreifenden Einzellast zu beschreiben, hätte man dann jedem 
Punkt a der Grundrißebene (x, y-Ebene) einen bestimmten Punkt m mit den Koordinaten 
&, n, £ (Einflußwerte von ®P für Lasten in der .-, y-, z-Richtung) zuzuordnen, dessen Pro- 
jektion ın’ auf eine Parallele zu P durch a in der Strecke am’ wieder die Größe ® für die 
betreffende Lastrichtung darstellte. Entsprechend der bei den ebenen Gebilden verwendeten 
Kreisschar hätte man hier durch jeden Punkt der Grundrißebene eine Kugel mit dem Durch- 
messer am zu legen; die Größe der von der Kugel auf einer durch «a parallel P° gelegten 
Geraden abgeschnittenen Sehne gibt dann die Größe von P für P=1 an. Die Erläuterung 
der Beziehungen zwischen den &, n- und £-Flächen möge einer allfälligen künftigen Arbeit 
vorbehalten bleiben. 162 
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Interpolationsverfahren zur numerischen Integration 
gewöhnlicher Differentialgleichungen. 
Von @ünther Schulz. 
\us dem Institut für angewandte Mathematik der Universität Berlin.) 


ür das Adamssche Verfahren fortlaufender Extrapolation zur rechnerischen Ermittlung 
F des Integrals einer gewöhnlichen Differentialgleichung hat R. v. Mises vor einiger Zeit!) 
eine exakte, bequem anzuwendende Abschätzungsformel für den begangenen Fehler angegeben. 
Im foleenden soll zunächst eine Abänderung dieses Verfahrens, die die Kxtrapolation durch 
eine Interpolation ersetzt und daher eine noch bessere Approximation erwarten läßt, geschildert 
und für sie eine Abschätzung des Fehlers durchgeführt werden. Es wird sich zeigen, daß 
man in der Tat zu einer wesentlich niedrigeren Fehlerschranke kommt, als wenn man nach 
dem ursprünglichen Rechenverfahren in der Form, wie es M, S. 85 zugrunde gelegt wurde, 
rechnet. Zur praktischen Durchführung der Rechnung, die nur wenig mehr Zeit als beim 
oenannten Extrapolationsverfahren kostet, werden vorteilhaft Iterationen verwendet. Be- 
inerkenswert ist, daß das sog. Adamssche Verfahren samt der hier behandelten Abänderung 
in nicht wesentlich verschiedener Gestalt, aber ohne jede exakte Fehlerabschätzung, seit 


langer Zeit vielleicht auch schon vor der Anwendung durch J. ©. Adams im Jahre 1883 
auf die Bereehnung der Meridiankurven von Tropfen?) — in der rechnenden Astronomie zur 
Bahnbestimmung der Himmelskörper benutzt wurde’). Schließlich wird im Teil II die 


Fehlerabschätzung auch noch auf den Fall ausgedehnt, daß man sich die Punkte eines Stückes 
der Interralkurve, das zu Beginn der Anwendung des Adamsschen Verfahrens schon bekannt 
sein ınuß, nicht dureh Reihenentwicklung, sondern ebenfalls mit Hilfe der gleichen Interpolations- 
formeln und Anwendung von Iterationen verschafft. Für die Konvergenz werden hinreichende 


Bedineungzen angegeben Soweit in den Abschnitten 1, 3 und 4 einige der nötigen Entwick- 
luneen analoge denen in der Arbeit von v. Mises verlaufen, sind sie hier nur kurz an- 


eedeutet worden. 


l. Das Verfahren fortlaufender Interpolation. 


I. Allgemeiner Ansatz. Ist die numerische Integration einer Differentialgleichung erster 


„ f’ | 
A A 


deren rechte Seite der Lipscehitzbedingung 


fir. fl, my) x»|y--y" a er 


venüge, so weit vorgesechritten, daß zu gewissen mit der Schrittweite h äquidistanten Abszissen 
Eos Ks... 27, anstatt der genauen zugehörigen Ordinaten y(#,), la)... ., 9(2,) der gesuchten 
Integralkurve tr) die Näherungswerte %,, Yı3 +: :, Mu bereits ge fanden sind, so liefert das 
\damssche Verfahren fortlaufender Extrapolation einen weiteren Punkt (#„.;,, Y„.4,) in 
foleender Weise. In der exakten aus (1) durch Integration entstehenden Formel 
Yn+i 
Yln+ı) Y\ın) \ fie, y()| er. N RE ARE (9) 
Xn 
wird der Integrand «durch das Polynom r-ten Grades ersetzt, das an den Stellen #,, @n-ı3 ++» 


r,_r die Werte f(a,, 4). Sn» Mn) =: +» Flan-r: Yun r) annımmt, und die linke Seite dureh 
Yyıı Mu. Die entstehende Formel, die für die Fehle rbestimmung in der Gestalt 


X 


Hn | Ur . h = Ar Lu fan ry Yı = ‚) . . . . . . . . . ( | ) 


benutzt wird, dient dann zur Bestimmung von Y+, 


Dieses en kann nun dahin abgeändert werden, daß man gleichsam als wäre 
die gesuchte Ordinate 7,.;, schon bekannt den Integranden durch das Polynom r-ten Grades 
ersetzt, das an den Stellen «,. inet Daattı ME EEE Filanız Mazı)s Fin, Mu), 

Ih Zeitsehr. f. angeew. Math. u. Mech. 10, 1030, 8. 81 bis 92. Diese Arbeit soll der Kürze halber hier stets mit „M“ 


tiert werden. 
°», Fr. Bashforth and J. C. Adams, An attempt to test the theories of eapillary action, Cambridge 1883. 
), Verl. B. C. L. Charlier,. Die Mechanik des Himmels, Bd. ?, 2. Aufl., Berlin und Leipzig 1927, insbesondere 
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f(@n-rtı> Yn_r;,) annimmt. Auch hier wird zur Fehlerbestimmung das Polynom in der Form 
der Lagrangeschen Interpolationsformel angesetzt, für die eigentliche Bestimmung der 
Integralkurve in der Newtonschen Form, die später angegeben wird. In analoger Weise, 
wie man zu Gl. (4) gelangte, erhält man jetzt aus (3) die Formel, die das „Verfahren fort- 
laufender Interpolation“ kennzeichnet: 
ä 
Yntı Inh are Fantı-e Intı-0) >» > nn.) 
0=V 


Dabei sind die a,." wohlbestimmte, nur vom Näherungsgrad r abhängige Koeffizienten, deren 
Werte noch angegeben werden. Der gesuchte Wert %y„;, kommt in (5) zweimal vor, sowohl 
auf der linken Seite als auch auf der rechten Seite als Argument der Funktion f. Es wird 
allerdings nur in Sonderfällen möglich sein, Gl. (5) exakt nach „+, aufzulösen. Der Bau 
dieser Gleichung legt es aber nahe, y,„;, durch Iteration?) zu berechnen. Man verschafft sich 
dazu zunächst einen „vorläufigen* Wert %,„.;," von Y+,. etwa mittels Gl. (4) mit dem 
Näherungsgrad s 


(v0) 


S 
Y ’ . 
Int In h > ao lan. Be 5 re Ne (b). 


[2] 0 
® 


(0) 


Wird dieses „+, 
man aus 


auf der rechten Seite von Gl. (5) an Stelle von y„;, eingesetzt, so findet 


ae en % f/n [0 Bein *f/. Pan 
Anti Ynzhı a, ; en+ı, Yntı ) Far, Jan Yu) t--- + Gr antı-r Yntir)] (“) 
einen neuen im allgemeinen besseren Wert ya4+,ı, kann fly: Yu4 1) berechnen und 
in gleicher Weise ein %,+,” finden. Allgemein gilt für die Iteration 
(41) 
| ER Ef(. a & 
Y, 14 - h lür. AG, n+1» Int ) r Arı f(4 n» In) | 


TE Mu (9). 
w +a,, flanı; r» Yaztı ‚)] (vr v, = .. | 


Für die praktische Rechnung werden an Stelle von (6), (7) und (8) die entsprechenden 
A-Formeln benutzt. Konvergiert die so gefundene Folge y,:," , Ya  » Musi +... gegen einen 
Wert y„+,, so genügt dieser der Gl. (5). Dies ist der gesuchte Näherungswert für (a4). 
der zur Bestimmung weiterer Ordinaten Y,:13; Yun 4s.... der Integral-Näherungskurve benutzt 
wird. Führt man noch die Fehler 


De a (+) 


ein, so gelangt man in analoger Weise wie beim Verfahren fortlaufender Extrapolation 
(M, S. 83) zu 





? 
y Mn | 
Ey | Ey h 2. A, Flau+; 0 Yatı 0) fau+, 0 Kurt o)) 
‚—l , 
(10). 
/ Aura 
HL ar Flrntı-, Yan) | Fa, yaa))diu 
0=V Ku 


Wird die obere Schranke des Absolutbetrages des zweiten Gliedes auf der rechten Seite für 
alle in Frage kommenden Intervalle — der „Quadraturfehler* mit c,* bezeichnet und die 
Lipschitzbedingung (2) benutzt, so folgt 


a 
eu, — &,|Snh >’ |." |’ |en4, etc" Dr BE ES 
ren 
u 
oder 
1 
eu, 1 —»hla, |]es || I +shla,,*||Hah N are | |enrı-el+c*. (12). 
Damit (12) brauchbar bleibt. ist zu verlangen, daß 
a 5 ne ehe ao 


eine Bedingung, die sich später erneut ergeben wird. 


4) Zum Iterationsverfahren vgl. die Arbeit von R. v. Mises und H. Pollaezek-Geiringer, Praktische Ver 
fahren der Gleiehungsauflösung, diese Ztsehr. 9, 1929, S. 58 bis 77 und S. 152 bis 164. 
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Wie ein Vergleich der Gl. ®t, (S) und MW, (10°) einerseits mit (5) und (12) andererseits 
zeigt, wurde beim ursprünglichen Adamsschen Verfahren beim Rechnen mit dem Näherungs- 
gerad r, d.h. bei Benutzung von Parabeln r-ter Ordnung, aus den r-+1 Ordinaten %,, Yn-ı; 

Yu r die nächste Ordinate y„;, bzw. aus den r—+1 Werten von |e az 
eine obere Schranke für |e„;,| gewonnen. Hier ist dagegen beim gleichen Näherungsgrad r 
nur die Kenntnis der r {nicht #-+1) vorangehenden Ordinaten Y,, Yn-ıs == :+ Mn+ır bzw. der 
r Werte von | 8,1. |&yıls 2+.. | Pntıor)| nötig, man braucht also auch vom Ausgangskurven- 


stück einen Punkt weniger zu kennen. 


€ € 


„ n I 


Läßt man in (12) das Gleichheitszeichen gelten, so ist durch diese Gleichung zusammen 
mit den gegebenen r Fehlern |, |, |&, |: »-:, |&r,_, | der Anfangswerte eine Folge Y, von oberen 
Schranken für die |e,| bestimmt. Wir schreiben also 
ha. Net itahla., Fa I Far tr 


[n. 
+) 


Mit der Anzabe «dieser Formel ist die Fehlerabschätzung erundsätzlich erledigt. 


2. Konvergenz der Iterationen. Die Berechnung von Y,„,;, mit Hilfe der Gl. (6) und (8) 
nd » ® .. ® a .. * . . } 
erscheint auf den ersten Blick mühevoll zu sein. Indessen genügt es meist, in (6) s=0 zu 
wählen, und da «a, nach WM, S. 85 den Wert 1 haben muß, zu setzen 


"ou 


(0) 


Hr ee 


Dadureh kann die Konvergenz der anschließenden Iterationen höchstens etwas verzögert werden. 
Aber auch dann sind in der Regel nur noch 2, höchstens 3 Schritte nötig, bis sich der jeweils 
aus (S) berechnete Wert von y,„;,” innerhalb der gewählten Stellenzahl nicht mehr ändert. 


Um zu hinreiehenden Bedingungen für die Konvergenz der Iterationen zu gelangen, 
sehreiben wir Gl. (8) einmal für vr und einmal für » -1I und ziehen die entstehenden 
(leiehunzen voneinander ab: man erhält dann’) 


AT (V 1) 
(cv) * . 5 (vr) ® M . 
Yan Anti ha,, (u n+ı» In+tı ) Man+ı 4,4, . . . . . . (16) 
oder mıt Benutzung von (2) 
vr: ev) = cv) „ . m 
Y.ı, Hz; „Aa, Wu; Y,,, 8 


(a,, wird sich nämlich für alle r als positiv erweisen.) Macht man nun die Voraussetzung 


a an tr er A 


die sich sehon einmal als nötig erwies, dann konvergiert die rechte Seite von 


(cv) 1) (ı) ’ 


da Ya Fe llaaı  — Yarı It las Yaaı)t > . . (8 


"—>% 


absolut gegen einen Wert 9,.,. der die Gl. (5) befriedigt. Da gezeigt werden wird, daß 


RE re A en A a rar 
—>I 
so folet, daß dureh Wahl eines größeren r der Faktor zha,," beliebig verkleinert und damit 
die Konvereenz des Iterationsverfahrens verbessert werden kann. 


3. Der Quadraturfehler. Der (Juadraturfehler e,* wird in gleicher Weise wie beim Extra- 
polationsverfahren [M, S. 853°)] berechnet, nur daß statt der Polynome @s(w) und der Ko- 
effizienten a,o und ?, hier andere Polynome @;*" (u) bzw. a," und P,“ auftreten. Wir führen 
« statt .r als unabhängige Veränderliche durch 

] Bu nu . 
ve = &.,t+hu, u a Ge eG a 
7 ) 
ein und setzen 
2, + 5er ne rn 


Die Verwendbarkeit von Gl. (16) wird in 5 erörtert. 
ersetzt. 


0434 


Die dort eingeführten Bezeichnungen ©, (#) und 3, werden vorteilhaft durch %, |, und > 
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Ferner sei 
a. (wert oa (u)=u-—I1 
9 (u)=(un — 1I)u (un) (un Iulu+|) 29) 
9 (Uu)= (un Yuw-+l)...(wm+o 2) 





Dann gilt nach der Newtonschen Interpolationsformel 


r 

y \ “5 , Pr | . “ 

Fin) Fi): ou) + iu CE Ed ı )  ; -; \ 
A o! " (r -])! J I 


0 —( 


mit ‚+1 " I, da « nur im Intervall ®©...1 für uns in Frage kommt. Sämtliche hier 
auftretende, mit | bezeichnete Differenzen sind „rück wärtsgenommene*“ *) und durch die Gleichung 


1? Fu) le! Fin) je" Fu 1) a le 00 eV 


definiert. Schließlich folgt für den Quadraturfehler, da ,,,"(w) zwischen VO und I sein 
negatives Vorzeichen nicht ändert, 


a Ze ER te. 5 > ABER 
Dabei ist gesetzt 
Po’ = \ 90 (u) du u Be ee er I (26). 
“m 
Aus Gl. (26) findet man für die Po (o=V,1,..., 6) die ın der folgenden Zahlentafel angegebenen 


Werte. Die beim Extrapolationsverfahren auftretenden Koeffizienten 5. sind zum Vergleich 
daneben gesetzt°). 





1 | 2 
I’o I’ | 
| 5 | ” 
Pı 5) 0.» P, ) 0). 
Rn | ERBEN " n 
P, 1» 0,0833 D, 1» 041664 
e | . > BR 
Ps 34 004165 Ps re VD ’ 
(4) 
19 2 
RR 4 +) - 0 )4j* 
P, 0 002659 I rt 0.345861 
A s 001875 | 032986 
5 160 . a Ep Pr ISS due St) 
al S63 001427 i 19087 031559 
de 0480 ee Peg 


Wie aus (22) und (26) folgt, haben alle Po*(o=1,2,...) negatives Vorzeichen; sie nehmen 
dem Betrage nach mit wachsendem o ständig ab und zwar, wie ein Vergleich mit den Po 
zeigt, wesentlich schneller als diese. Leicht läßt sich beweisen, daß die ?o und Po* füro—a& 
gegen Null konvergieren. Wir schreiben 


Bo \l1 | t hi | er | - au E (28) 
und i 
Bo* (1 4 ! ‘)( 2 z n | ir 2 = *) a, 
0 - / 


’, Vielfach, z. B. bei J. F. Steffensen, Interpolation (Baltimore 1927), werden sie auch mit 7 bezeichnet. 
>, Bei Bashforth und Adams 8. 18 findet man die Koeffizienten 7, und ?,* bis o=4. 
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Nun berücksiehtirge man, daß das unendliche Produkt li eeeen Null geht und zwar 
N () 


gleichmäßig für alle positiven » in jedem endlichen Intervall, das den Nullpunkt ausschließt. 
Zerspaltet man das Integrationsintervall in (25) in zwei Teile 0...1—- y„ und 1—n...1, so 


ist bei beliebig klein vorgegebenem & für alleo - P der Integrand des ersten Integrals < 5 : 
€ . 
das Integral selbst also „(1 1); ersetzt man im zweiten Integral alle Faktoren durch 1, 
so erkennt man, daß es n ist. Man findet so 
> € . 
0 .) (1 1) T 1, . . ; . . . . ’ ‘ . . . ö (30). 
Wählt man schließlich noch „= „. so folgt #0 8. Analog verläuft der Beweis dafür, 
2 u 


daß lim 80° VO ist. Hier ist aber die Zerspaltung des Integrals in (29) unnötig, da /J(1 

> « oe \ G 
im Intervall 0 u I gleichmäßig geren Null geht. Gleichzeitie wird damit verständlich, 
warum die 80* so viel schneller gegen Null konvergieren als die do. Schließlich ist jetzt auch 
die früher aufgestellte Behauptung 


ı) 


ee a ae ri 
> 


bewiesen, «da aus der Darstellung des Interpolationspolynoms in der Lagrangeschen Form sıch 
l 


| 
| Pi)... de+r—1 | 1 — u\ 
a, |", ) (n } 2. (fi “ 1 "au ß 


nm 
— 
.. 
ww 
-— 
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= en r | Bu r 
ergibt. Aus der Identität der Interpolationspolynome r-ten Grades in der Lagrangeschen 
und Newtonschen Form folgt nach Integration 

r “ r a, 
y FU—- der HEN). . : +. 8. RR. 
u" 0: 0 
Durch Vergleich der Koeffizienten rechts und links erhält man den gleichen Zusammenhang 
zwischen den «a,.* und den %0*, wie er zwischen den a,. und Po bestand, nämlich 


RR SE | 5.2 2: 5 


Die Zahlenwerte sind: 





PB: | 
| E l 
} hi: = ü; = 
.) R S r | 
) + Ga Pa Ass r j 
12 ” 12 e 12 (34). 
> ” ) er | ‘ ) . . ie | 
‚ >. Ayo I ’ Az, If ’ Ass 4 ’ Ag 4 
251 646 264 106 N 19 
; - A F _ . l -_ . Bi _ . ee -A, 
"4: de ng Ar = 790: Cs 720 u 720 “ 720 
! er} . * n4 ” 
Auch hier gilt %'a,.* 1. Die Tafel der a,.* verleitet leicht zu der Vermutung, daß, im 


IE 
Gerensatz zu den a,.. die Beträge aller Koeffizienten unter I bleiben. Dies ist jedoch nur 
. „ 2713 - > 
für r 6 richtig; dagegen ist z. B. a," = 5-59: also I. Für die Summe der absoluten 
| A Ye 
Beträge 


Y 67 >. .)- 
\ Aro A, . . . . . . j . . . . (39) 
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ergibt sich: 
57 1286 





A mei, 2° . — 1,1667, A,*: 4 1,1167, 4A,* 0° ‚861 (86) 
und für die Quotienten |ß,:,*|:4,*=y,*, die später gebraucht werden, 
y,*— 0,083 33 v,—0,2083  ] 
y,* — 0,085 71 y,— 0,1023 r 
y,* 0,018 63 y,— 0,052 29 EN. 
y,* — 0,010 50 y, 0.026 94 





Zum Vergleich sind die Werte von f,+,:4A,=- y, danebengestellt, wie sie das Extrapolations- 
verfahren benutzt. 


Zum Schluß sei noch die Formel angeführt, nach der die eigentliche Rechnung (nicht 
die Fehlerabschätzung) durchgeführt wird. Sie folgt aus (5), (20), (21) und (26) und lautet: 


+1 a r or . (v) 
J FF a In h Ri Po MR flan- ı» Yn4ı ) h Ent Yn+ı ) 
ee | \ >. (88). 
y >= Eu ) er | 
„ AfG..) po rTG..) 34 Te. .) og A a 





Während des zu einem Schritte gehörigen Iterationsprozesses, der Yu,” , Ynı" »... und 
damit auch f(&, Yu+1" )» F (&, Yusı ), --. liefert, ist das Differenzenschema immer wieder von 
neuem durch Berechnung der Arf(@u+1, Yn+ı )v=0,1,2,...) zu verbessern, bis innerhalb 
der der Rechnung zugrunde gelegten Stellenzahl keine Änderungen mehr eintreten (vgl. das 
Beispiel in 5). Da, wie bereits erwähnt, bei der Iteration meistens 2 bis 3 Schritte genügen 
und die Änderungen, die das Differenzenschema erfährt, dann nur in den letzten 1 bis 2 Dezimal- 


stellen liegen, ist die Rechenarbeit geringer als es auf den ersten Blick scheint. Um aus der 
v-F1) = .. 2 3 v-+1) 
) die Anderung von 9y,,,” in y zu be- 


Anderung von f(#4+1, In. ) in flausıs Y n-+i 


rechnen dient Gl. (16) (vgl. dazu auch 5). 


n-Hi 


4. Die Fehlerfortpflanzung. Anstatt aus der Differenzengleichung (14) die einzelnen Fehler- 
schranken Y, sukzessive zu berechnen, kann man auch das von v. Misesa.a.(). benutzte 
„Verfahren der majoranten Partikularlösung* auf (14) übertragen, um eine Schranke für |e,| 
ohne Kenntnis der vorhergehenden zu gewinnen. Wir haben also eine solche Partikularlösung 
Y, von (14) aufzusuchen, die sich aus einer partikulären Lösung der inhomogenen Gleichung 
und einer partikulären Lösung der zu (14) gehörigen homogenen Gleichung additiv zusammen- 
setzt und für v=0,1,...,r—1 Werte ergibt, die größer oder gleich den vorgegebenen Werten 
von |&,|, |&ils »..,|&r_,| sind. Wegen der Positivität der Koeffizienten von (14) majorisiert 
dann die so bestimmte Lösung jene Lösung, die für v—0,1,...,r—1 genau die Werte |e,|, 
l&,|,.- .,|er_,| annimmt. Die beiden fraglichen Partikularlösungen sind — e,*:#h A,* und (z”, 
die gesuchte majorante Partikularlösung mithin 


H 


( 


Y ,— Oz” r dr . . . . D . . . D . 0 
’ “hA,“ : 


(39). 


Dabei ist z die kleinste reelle Wurzel der Gleichung 
a — e-"—=xhlla,,"|-2r+le,, "2" +... + ll - - -» 0.080. (40), 


die >1 ist. Daß es stets mindestens eine solche Wurzel gibt, folgt daraus, daß für 2-1 die 
linke Seite von (40) gleich 0, die rechte gleich <h A, ist, und für große z die linke Seite 
wegen xha,*<-1 stärker unendlich (allerdings von gleicher Größenordnung) wird. Die Kurven 
der Abb. 1 gestatten zu gegebenem xh für r—1,2,3,4 diesen z-Wert aufzufinden, der dann 
nötigenfalls noch mit Hilfe von (40) genauer zu bestimmen ist. Die Abbildung zeigt zum 
Vergleich in einem x h-z-Koordinatensystem die Kurven für das ursprüngliche Extrapolations- 
verfahren (gestrichelt) und das abgeänderte Verfahren (Interpolationsverfahren). Man sieht, 
daß dieses bei gleichem #h und r einen wesentlich kleineren, also günstigeren z-Wert liefert 
als jenes. Unglücklicherweise nimmt bei beiden Verfahren, wenn man bei festem xh zu einem 
größeren Näherungsgrad r übergeht, der Wert von z auch zu, so daß die Fehlerabschätzung 
(bei genügend großem n) schlechter wird, obwohl in der Regel die Rechnung genauer ist. 
Dieses Anwachsen des z-Wertes bei zunehmendem r ist jedoch bei den ausgezogenen Kurven 
wesentlich langsamer als bei den gestrichelten. 
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Weiß man nur, daß |e,|,|& |,..., &r_,, alle unterhalb e liegen, so ist die günstigste 
Wahl für © 
(' Wr. 
ie. 0 ee re (41). 
Somit ıst also 
v1 e nv ») 
P, e2 +7, „5 @ re | 
oder weren (25) 
£, gz’ 7 in ] / - we R I) (43) 
T4r T a T ie | 
/ r=:h/ IT r 
| Y | | 7 A 
| / | Pa s 
RR tt a Er 
zZ | | 4 I 
/ | wa Pe 
/ | + | or 
13 -—,., > 4—r 
# Pa Fr | 
/ ri | ‚7 
91:5 
/ _ | vr 7=f,,- 
/ / | wei a 
12 / # en ni > 1 | 
/ v ‚ | Be | | 
/ Pi ”ü | et, | J 
I L > M et. I, T Eee 
/ 
v ' 
, 4 | 
77 rar 4 — 
/\/ | 
/ | 
N san; 
// + P— | | 
NV | 
„es EA IE 
/ l l En u Re = N rn BRUNS zn, 
u: 002 004 006 008 0n ar UN 


Abb. 11. 


Kin Vergleich «dieser Formel mit jener für das Adamssche Extrapolationsverfahren (Gl. Dt, 
(32)). unter Heranziehung der Zahlentafel (37) und der Abbildung zeigt, daß sie sich nur in 
den auftretenden Konstanten unterscheiden. Bei gleichem r ist das z und die Konstante 
Pr, Ar" in (43) kleiner als der entsprechende z-Wert und die Konstante P,;,:4A, ın 
»t, (32), mithin bei gleichem &, »v, h, f (r, y). 2 die aus (45) gewonnene Fehlerschranke für | e, | 
eünstieer. Um welche hohen Beträge sich bei großem » die Fehlerschranken unterscheiden 
können, zeigt am besten ein Beispiel. 


5. Beispiel. Als Beispiel sei auch hier die Integration der Differentialgleichung 


dy „ x 


d.ı BE „ (44) 


mit dem Anfangswert #0, 9-1 mit A=002 und r t gewählt (vel. Mt, S. S9 bis 92). Das 


folgende Rechenschema zeigt, wie bei Anwendung des Iterationsverfahrens %, gerechnet wird, 
wenn YoNYır ++; y, oder, was dasselbe ist, wenn f(x, , 21): If,» 9): I’ f(@,. 9,) bekannt 
sind. Die Formel 
’ | | FA a 
Ur >, Hu I (tn; Hu) rs‘ l/ (Ms Yun) T 12 Fflau Yn)+ ha "ilzu, Yu) 
(49), 
ee 
T mo | Kan; In) 
v0 
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nach der beim Adamsschen Verfahren in seiner ursprünglichen Form gerechnet wird, liefert 
mit dem Näherungsgrad 4 für n=5, «0,10 den Wert 7,” = 1,091126. Mit Hilfe dieses 
Wertes wird sodann der Funktionswert hf(xr,,»,'”) berechnet und das Differenzenschema 








N r Y h } hAf- 10%  AA?F- 10%  AA®T- 10% | AArF- 10% 
4 0,08 1.074195 = Y, 0.017 228 623 4» 7 0 
5 0.10 1.091 126 — y,(") 0,016 642 586 37 -5 > 
D 0,10 1,091 127 ey, ) = y, 0,016 642 586 37 5 > 


dureh Bildung von hAfla,, y,"),..., AN f(r,,y,”) vervollständigt. Die Gl. (35) des Inter- 
polationsverfahrens liefert, ebenfalls beim Näherungsgrad 4, 9,’ 1,091 127. Sodann wird ein 
verbesserter Funktionswert hf(@r,,9y,"’) und ein verbessertes Differenzenschema h Af(r,, y,”). 
.., hNfi#,,9y,”) gebildet. Man darf „,“’ schon als endgültigen Wert „, betrachten, da Be- 
nutzung von Gl. (16) zeigt, daß sich hf(a,,y,”) von hf(#,,y,"’) innerhalb der gewählten 
Rechengenauigkeit nicht mehr unterscheidet. 

Die Rechenarbeit ist noch etwas geringer, wenn man den vorläufigen Wert ," nur mit 
dem Näherungsgrad 0 berechnet, also Gl. (15) benutzt. In diesem Fall erhält man für 
n—=5, «0,10 den Wert „, = 1,091423. Nach Berechnung des zugehörigen Funktionswertes 








N r y hj hAf- 106  nA?f- 106 m A87 10% nA} - 10% 





* 4 0,08 1.074 195 — y, 0,017 228 623 42 7 0 
5 0.10 1.091 423 — y,(") 0,016 648 585 38 4 3 
5 0.10 1.091 27 = y,)—y, 0.016642 586 37 5 2 


(ı) 


und des Differenzenschemas ergibt Gl. (35) mit dem Näherungsgrad 4 den Wert , 
— 1,091 127. Daß auch hier ,“’ schon für den endgültigen Wert , gesetzt werden darf, 
>51] 


_ 


zu 349 ro] : 
20 0,349 ergibt: 


zeigt Gl. (16), die hier mit a," — 
a Y — () 349 |0,016642 0,016643] I — 10°. 
.) 


Durch * sind im Rechenschema stets die endgültigen Werte der Punktfolge (.“,, ,), die man 
für die gesuchte Integralkurve erhält, kenntlich gemacht. 
DO ! 


Es werde nun für dieses Zahlenbeispiel |e,,| nach Gl. (43) abgeschätzt. Da r=4 und 
“h=0,W06, so liefert uns das zh-z-Diagramm den z2-Wert 1,01, den man mit Hilfe der 
Gl. (40) in 2= 1,011 (etwas zu groß) verbessern mag. Daraus folgt z'’= 1,116. lerner ist 
ft = 8540, h= 0,02, Ai? 64- 10°, 1,419: 4,°= 0,0105. Mit e=0,5- 10% ergibt sich 


0.5:10°%.1.116 0.0105 10T ago 0.116 
Es is. 98» 1, rV, .) 0.006 h U), ) 


0558-10 %+0,115 - 10% — 0,673: 10%, 





Dagegen hatte man DW, S. 91 gefunden: | 8,, 3,27 10%, also eine 5mal so große Schranke®). 





Die Gl. (43) für die Fehlerschranke ist deswegen besonders wertvoll, weil sie auch für 
große Schrittzahl » Schranken ergibt, die den Fehler, der im ungünstigsten Fall begangen 
wird, nicht allzu sehr übertreffen. Beim Adamsschen Extrapolationsverfahren ohne Iteration 
brachte der Übergang von Gl. M, (24), nach der jeder Fehler sukzessive einzeln berechnet 


9, Es sei bemerkt, daß dureh Abrundungsfehler, die ja nur in den Anfangswerten (£e=,+10 6), aber nieht in 


der Rechnung berücksichtigt sind, die eben erhaltene Schranke 0,673 +» 1076 etwas liberschritten werden kann. Bei der 
l * 
D 
erößerer Stellenzahl durehgeführt werden, so daß die Zahlenwerte nieht dureh Abrundungsfehler beeinflußt werden. 

4* 


vorliegenden Untersuchung des „Formelfehlers“* müßte eigentlich unter Beibehaltung von : 106 die Rechnung mit 
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werden muß, zu Gl. ®t, (32), die die Kenntnis der früheren Fehler nicht ‚voraussetzt, eine so 
erhebliche Verschlechterung der Fehlerschranke, daß M, (32) für große v, wie das folgende 
Beispiel zeigt, kaum angewendet werden kann. Der gleiche Übergang beim Interpolations- 
verfahren von Gl. (14) zu (43) bringt zwar für die Fehlerschranke auch eine Verschlechterung, 
die aber wegen der dort auftretenden kleineren Koeffizienten kaum ins Gewicht fällt, so daß 
(43) auch für große » praktisch anwendbar bleıbt, wo ein schrittweises Berechnen der Einzel- 
fehler zu mühevoll und zeitraubend wäre. Das Adamssche Verfahren wird nämlich oft für 
Scehrittzahlen von der Größenordnung v = 100 und darüber benutzt. 


Denken wir uns nun das eben behandelte Zahlenbeispiel nieht nur bis » = 10, sondern 
bis » = 100 fortgeführt und berechnen wir die Schranke für |e,,, aus (43). Man hat z= 1,011, 


z'° — 2,99, also mit den oben angegebenen Werten: 
” | „64-10 
€ 05:10 °°%.2,99-+- 0,0105 SS40 . 1,99. 


100 


0,006 


1.50: 10% +1,97 - 10° = 3,47 - 10°%, 





Die Anwendung der Fehlerschrankengleichung W, (32) für das ursprüngliche Verfahren würde 


für dasselbe Beispiel folgende Zahlenwerte liefern: z = 1,067, 2" 655, 


+ 12 
e . 


+1 
0,5: 10°°%.659 + 0,0269 SS40 : 654 


i 0.006 


100 


(3285 + 1662) - 10 = 19% - 10°, 





also einen 574 mal so großen Wert wie beim Interpolationsverfahren. 


II. Die Bestimmung der Ausgangswerte. 


6. Hilfssätze. Um den Gedankengang der nun folgenden Überlegungen nicht unterbrechen 
zu müssen, seien zunächst zwei Hilfssätze vorausgeschickt. 


Satz I: Bestehen zwischen den Gliedern von k unendlichen Reihen 


U,o 7 U, T:+:-:TU,y 

Prien u U PT SEE nl Pe (46) 
) 

UkotUkıt::- TUkyr tr: 





mit positiven Gliedern (v,, 0) für» 0,1,2,... (oder wenigstens für alle » > N) die k 
linearen Ungleiehungen 


U,,r+ı > CC Wr TC Us» Tr HC KR ÜR» 


Us,y+rı > Caı Up T Caa Us . TC, k Uky 


(47), 





Un, ’+1 EC} 1 ", Ib + Cr 2 U, 4 + ae + Ck k tx y 


und besitzt deren Matrix & — (e, „) nur nichtnegative Elemente (e,„ = 0) sowie charakteristische 
Zahlen «u |[d. h. Wurzeln der Gleichung det (CE  «&)=0], die dem Betrage nach kleiner als 
| sind, dann konvergieren die k Reihen sämtlich. Dieser Satz ist offenbar eine Verall- 
gemeinerung eines bekannten Konvergenzkriteriums. 


Zum Beweis lasse man in (47) das Gleichheitszeichen gelten, schreibe @,, an Stelle von 
t,, und verlange noch ü,,= u,, (2° 1,2,...,k). Dann sind durch 


une rt ertren, Bell... vr, MA, 
wegen €,„-.0, k Zahlenfolgen “,, mit positiven Gliedern definiert, für die 
Be, EEE HERE ne a. en 
[7 
eilt. Für die n-ten Teilsummen s, '4t,,„ dieser Folgen läßt sich ein ähnlich wie (48) 


ı u 
zebautes, aber inhomogenes System von linearen Differenzengleichungen aufstellen, dessen 
explizite Lösung die Konvergenz der s,, für n>x in Evidenz setzt. Oder einfacher 
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man faßt ü, ,„, %, ,,...,@;„ als Komponenten eines k-dimensionalen Vektors u auf; dann schreibt 
sich (48) kurz in der Form 


0 te ee 


und man hat die Konvergenz der Vektorreihe 


wru,rw,+...=(E + +E’+..)W: : » 2.2.0.0... (Öl) 
zu untersuchen. Die Reihe der Matrizen auf der rechten Seite konvergiert aber nach be- 
kannten Sätzen der Algebra, sofern &E nur charakteristische Zahlen vom Absolutbetrag | 


T. 
hat, was hier auch vorausgesetzt war. Die Reihe \'u, konvergiert also und mit ihr die Reihe 
van 
Rn . 
Y’u,, die durch jene majorisiert wird. Die Konvergenz der Vektorreihen zieht die Konvergenz 
0 
L 
der aus den entsprechenden Komponenten gebildeten Reihen \'’«,, nach sich, womit die Be- 
hauptung bewiesen ist. 


Satz 11: Bestehen für je k reelle Zahlen &,,&,,...,&; die 2% Ungleiehungen 


a, ta: &+::.: ta, &—b, 0 


(l,, &, 4 (A, E. = en + d,R &; b, . 0) FIO 
N et a A ee La Der 





a, +tar,&+.. + bu <0 


0,80 ee DAR, Fe VE 


be 
-i 5) > ser. 6C 


wobei die Determinante der a,, nicht verschwindet und a,, >0,a,,„<0(=*x) und b, >0 ist, 
dann erhält man obere Schranken X,,X,,..., X; für alle diese Zahlengruppen, falls das zu 
(52) gehörige Gleichungssystem 


4, Sı Aa St... TAık Sk b, 


d,, &, + (,, 7 + ..r + Ay,R Er b, _ 
ee 





Ar, di + Aka &,-+. ; . + Ad;% Er. h,. 


eine positive Lösung hat. Ist also X, >0,X, >0,..., X; >0 diese Lösung, dann gilt 


I<E LEX HS, 


1: (55) 


für alle Zahlengruppen &,,&,,..., &, die den Ungleichungen (52), (53) genügen. 
































Abb. 2. Abh. >». 


Die behaupteten Aussagen veranschaulicht man sich vorteilhaft geometrisch, indem man 
die £, als rechtwinklige Koordinaten in einem k-dimensionalen Raum deutet. Jede lineare 
Gleichung zwischen den £, stellt eine (k# —- 1)-dimensionale Hyperebene dar; auf ihr liegen 
alle Punkte, deren Koordinaten £, die Gleichung befriedigen. Die Hyperebene selbst teilt 
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den Raum in zwei k-dimensionale Halbräume, in denen die Punkte liegen, deren Koordinaten 
£, der Ungleichung genügen, die aus der Hvperebenengleichung entsteht, wenn man das 
Gleichheitszeichen dureh bzw. ersetzt. Sind für die Koordinaten £, mehrere Unglei- 
ehungen vorgeschrieben und gibt es £,-Werte, die diesen genügen, so bedeutet dies, daß die 
zugehörigen Punkte in einem (endlichen oder unendlichen) konvexen Gebiet liegen müssen, 
dessen „Stützebenen* die zu den Ungleichungen gehörigen Gleichungen (oder wenigstens ein 
Teil von ihnen) sind. 

Der Satz behauptet nun, daß die %# (sich im Nullpunkt o schneidenden) Koordinaten- 
hyperebenen &, 0, &,=0,...&,;,=0 und die k übrigen sich in einem Punkte p mit den Ko- 
ordinaten X,,X,,..., Xx schneide nden Hyperebenen einen endlichen, nicht verschwindenden, 
konvexen Te ı T des k -dimensionalen Raumes (konvexes Polytop) bestimmen, dessen Punkte 
sämtlieh Koordinaten £, = X, haben, der also ganz in dem durch O=Z2,<XN, (z=1,2,...,k) 
verebenen Quader O liegt, mit dem er die beiden Punkte o und p und Stücke der Koordinaten- 
hyperebenen gemeinsam hat. Die Abb. 2 und 3 zeigen das Aussehen des Raumteils 7 im 
Falle k=2 und k=3. Aus den Vorzeichenbedingungen für die a,, folgt übrigens, daß jede 
der # sich im Punkte p schneidenden Hyperebenen nur einen positiven Abschnitt auf den 
Koordinatenachsen ausschneidet, und zwar die #-te Ebene auf der £,-Achse. 

Der Beweis wird nun, wie folgt, geführt. Zieht man vom System (52) das System (54), 
das ja für &,— X, erfüllt ist, ab, so folgt, wenn zur Abkürzung £, X, = 7; gesetzt wird, 


4, 9 #9, N. +-.-- FR NY SV 

Os N, rt GN, rt. TEN SU du 
(6). 

Ar + Art. RN ZU 





Das Ungleiehungssystem (52) oder (56) bestimmt einen gewissen unendlichen Raumteil. (Im 
Fall #=3 ist es einer der S Teilräume, in die der Raum durch die 3 sich im Punkte „7,0 
schneidenden Ebenen zerlegt wird.) Um nun zu zeigen, daß alle Punkte des Raumteils (56) 
im Raumteil 

7 MEN te EN ae a ee 


71 72 


liegen, genügt es nachzuweisen, daß alle Punkte der Begrenzung von (57), also gewisse Punkte 
der Hyperebenen „©, außerhalb des Raumteiles (56) liegen und daß es einen Punkt gibt, 
der zugleich innerer Punkt von (56) und (57) ist. Die fraglichen Randpunkte des Raumteils 


(57) haben im »,„-Koordinatensystem die Koordinaten 


Wi Br rn 
oder N.,; 0, Nagy +; Nok 
oder 
oder Ku. ee 


wobei die »,, irgendwelche nichtpositive Zahlen sind. Wegen der Vorzeichenvoraussetzungen 
für die a,, hiegen alle diese Punkte er des Raumteils (56), da stets mindestens eine 


der U ngle ichunge n (56) verletzt ist. Da alle X, >0 und alle b, > 0, gehört der Punkt £, =0, 
das heißt 7, \, sowohl dem Raumteil (56) wie (57) an: aus N„=QV folgt aber &,=X,, 


womit der Satz bewiesen Ist. 


7. Allgemeiner Ansatz. Wie bekannt, können die Rechenformeln (45) bzw. (38) des 
Adamsschen Verfahrens erst dann angewendet werden, wenn bereits einige Punkte aus der 
Punktfolge, durch die die gesuchte Integralkurve approximiert werden soll, bekannt sind. 
Meistens verschafft man sich dazu die Tavlorsche Entwicklung von y (#) für den vorgegebenen 
Anfangspunkt (os Yo) und setzt dort zo. +h, a, +2h, @,+3h,... (Beispiel: M,Ss. S)). 
Hierbei können folge nde Nachteile auftreten: Die Potenzreihe konvergiert möglicherweise sehr 
langsam, so daß sehr viele Glieder herangezogen werden müssen, was besonders dann, wenn 
die Berechnung der Koeffizienten mühevoll ist, die Brauchbarkeit des Verfahrens zur schnellen 
Bestimmung einer Integralkurve sehr herabsetzt. Außerdem kann der Konvergenzradius der 
Reihe so klein sein, daß aus der einen Entwicklung nicht die erforderliche Zahl von Aus- 


gangswerten berechnet werden kann. Schließlich kann — insbesondere bei langsamer Kon- 
vergenz die Beurteilung schwer sein, wie gut die aus der Potenzreihe berechneten Werte 


die genauen Werte approximieren. Für eine exakte Fehlerbestimmung ist jedoch die Kenntnis 
des e in Gl. (41) unerläßlich. 

Eine andere Möglichkeit zur Bestimmung einer genügenden Anzahl von Ausgangswerten 
Yys Ms, .. ohne die eben geschilderten Nachteile bietet sich, wenn man auch hier schon die 
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Formeln des Verfahrens fortlaufender Interpolation, bei der praktischen Durchführung in Ver- 
bindung mit Iterationen, benutzt. Die unten angegebene Fehlerabschätzung gestattet die 
Beurteilung, um wieviel die so erhaltenen Werte y,,%,.... von den auf der genauen Integral- 
kurve gelegenen Werten y(#,), 9(#,),... höchstens abweichen; für die Iterationen werden 
hinreichende Konvergenzbedingungen aufgestellt. 

Wir wollen annehmen, daß zu dem vorgegebenen Anfangspunkt „, noch weitere k Punkte 
Yır Yas«+ +, Yr, die den Abszissen #,+h,x,+2h,...,: r,_+kh entsprechen, bestimmt werden 
sollen. Will man später mit dem Näherungsgrad r weiterrechnen, so ist bei Benutzung von 
(45) (fortlaufende Extrapolation) k=r, bei Benutzung von (35) (fortlaufende Interpolation) 
kr —1 zu setzen. Werden nun die gesuchten Werte %,.9s,..., 9, Schon als bekannt an- 
eesehen, dann kann man durch die k+1 Punkte (#,,f,): (@,:f})»-:-, (r, fr) die Parabel k-ten 
Grades legen (zur Abkürzung ist hier f(“#,,9,)=f„ geschrieben worden). Setzt man das Inter- 
polationspolynom in der Lagrangeschen Form an, so liefert Integration über die Intervalle 
(ll. er; Jeweils Ausdrücke, die in f,, f,: fa; - - -, fj, Jinear und homogen 


sind. Diese werden gleich 9, — 9, Ps — Has: Yr MY, gesetzt. So entsteht analoz zu (5) 
das folgende Gleichungssystem 


Hi Y%= h ls.fo+ Rn 5 +.:..T70% fi] 


Ur Yo h [930 fh +0,, f; F...40;% fi| (58) 





Y: ur h 57; 0 E ü r OR 1 f\ + + Ök «fi] 


Die o sind dabei wohlbestimmte Zahlenkoeffizienten, die sich ähnlich wie die « und a* be- 
rechnen; z. B. ist 


nl I)(m 3) 
ER \ Ida 


I:-(—-N:(-2) 5 für k=3. 


NE m 


u" 


Die Rechnung ergibt folgende Zahlenwerte für die rechteckigen Matrizen 





























Go 91 9,1 
Oso Gpı Il, 
Öko Or; Se Ö]; l; 
| | 
für k=1: 1; - | 
> S | 
12 12 12 
für k=3: 
| 4 | 
) 3 3 
4) 19 » | 
24 24 24 24 (39). 
is: 2 I - | 0 
u ur 3 3 3 
; ‘ı) ) | 
S S S S 
351 646 264 106 19 
120 120 120 120 120 
29 124 24 | l 
/ wo u) % ww %w 
für l = t: L PR. 
2 102 (2 42 ’, 
Ss0 Ss0 Ss0 Ss0ı Ss0 
14 64 24 64 14 
15 45 15 15 5) 
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Die Summe der Elemente der ersten, zweiten, ... Zeile jeder Matrix beträgt 1, 2, .... 


Die unbekannten Werte %,, %s, -.., 47, kommen sowohl auf den linken Seiten der 
Gleichungen des Systems (58) vor wie auf den rechten Seiten als Argumente der Funktion f. 
Auflösung nach dem Iterationsverfahren ist der praktisch allein mögliche, dem Bau des 
Systems auch durchaus entsprechende Weg. Dazu werden zunächst k „vorläufige Werte“ 
MY, Yo", 2... 9% bestimmt, entweder durch wiederholte Anwendung von (45) mit dem 
Näherungsgrad 0, d. h. von (15): 


= yHhh, ty Hl”, men th" 2. (60) 
oder mit von Schritt zu Schritt wachsendem Näherungsgrad 


P 1} (0) | r 0) I 1 , 0) 
a Yywthf: . ya =y,; ei Tr .) IA . 


| » | (0) | > 200 
4 y, Ri 1" di > Af, ) FR 19 ?f, . 





P (0) l y (0) N ) Zu (0) 
Ya Hk + h D 1 7 ‘) ıfı. 1 au 2 Be + Pr 1 1 fi: 1 | 


(Zur Abkürzung ist dabei f(#,, 1’) 


dann nach den Gleichungen 


- f,”” geschrieben worden.) Die Iterationen verlaufen 
ua 2. un hluhtrufi” +0: f” +:...+0;% 7 


. Fa a Y nn hi [930 f, + Oz, A + Öas Fake + eo.“ + Oak fe] (62) 
In 





Yo Pr” Un a h lo%, fi + Oki ya r Okas aa + nn + Okk r,. 


mit »—=0, 1,2, .... Bei der praktischen Durchführung der Rechnung werden natürlich an 
Stelle der Gl. (62) die entsprechenden J-Formeln benutzt, die später in (76) angegeben 
sind. Konvergieren die Iterationen, so liefern sie k Werte %,, Ys, -.., Yx, die den Gl. (58) 
genügen. 


8. Konvergenz der Iterationen. Denkt man (62) einmal für » und einmal für » —1 ge- 
schrieben und subtrahiert entsprechende Gleichungen voneinander, so entsteht 


yo: 1) Ara hilo,if,” -f’)+...+0xr (fr nr »)] 
Yo di 47° h [9,, A A . + ... + 05% "iu Br u 








. (63) 
ya? = hlon, fi” Fr) +... + on fe” — 9] 
und weiter nach Berücksichtigung der Lipschitzbedingung (2) 
wer yr|sahllo.l)y” N r9l+...+loRl 9” —yr®|] 
Y 7 — y,Prisahlo,))y” —-yri+...+loxl|%a” ua” "|] 6 
. . . +). 
Y% 0 — a |\sahlo,)|y” —-y’ri+... tor |y — Yu” |] 


Y 


Nach dem Hilfssatz I in 6 konvergieren aber die k auf der rechten Seite stehenden Reihen in 


y=imy,” =” +? u) Hy ur)+... Kl... .. . (66 
r—>% 


dann absolut gegen Werte %,. Ys, -:-, 4x, Wenn sämtliche charakteristischen Zahlen der 
Matrix xAh(o,,) dem absoluten Betrage nach < 1 sind, oder, wenn sämtliche Wurzeln der 
Gleichung 


Ö u OÖ R we O,k 


(66) 





































(61). 





d 
di 


m 
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dem absoluten Betrage nach - —E sind. Damit ist aber folgende hinreichende Bedingung für 


die Konvergenz des Verfahrens gefunden. Wird die Wurzel « vom größten Absolutbetrag 
mit 1: A bezeichnet, dann muß sein 


„h > Bee: 
A ist also der absolut kleinste „Eigenwert“ der Matrix (o,,), wenn man die reziproken 


charakteristischen Zahlen einer Matrix Eigenwerte nennt. Mit den in (59) angegebenen Zahlen- 
werten lautet Gl. (66) 


für k=1: Su I=0 05 

für k=2: Ye — Ju-+l 0 0,373 
fürk=3: 72. — 108 + 2361-30 1,255 
für k=4: 121500 u.” — 243000 u + 66345 1° —- 1226 u — 2220 — 0 1.6853 


In der letzten Spalte ist die Wurzel vom größten Absolutbetrag angegeben. Daraus folgt 
für Gl. (67 

h< 2 für k=1 | zh<0S10 für k=3 

h< 1146 fürk=2 zh< 059% für k=4, 


Je größer also die Zahl der Ausgangswerte ist, die man nach dem geschilderten Verfahren 
gleichzeitig berechnen will, um so niedriger ist die Schranke, unter der #h liegen muß. 


9. Abschätzung der Fehler. Nachdem nun für ein unterhalb der angegebenen Schranken 


liegendes #h die Konvergenz der y,”, 9”, ..., 95° gegen gewisse Werte %,, Mas +... Hk 
sichergestellt ist, sind schließlich noch für die k Fehler |»,| = y,  v(r,), z=1,2,..., k) 





obere Schranken anzugeben, deren man zur Abschätzung von |&;;,|, |&r+s|,; -.. bedarf und 
aus denen insbesondere die Konstante € in (39) bestimmt wird. Genau in gleicher Weise, 
wie (11) aus (3), (5) und (10) gewonnen wurde, subtrahiert man von der ı-ten Gl. (58) «= 1,2, ..., k) 
die Gleichung 


X 
wa yet : » ee 
xy 
fügt h N’ o,„fl@z,(@,)] einmal mit positivem und einmal mit negativem Vorzeichen hinzu 
%==() 
und erhält schließlich 
k ’ 
€, Eu, ah er O,2\'\f% + d,’ . . r s . : . . : F i ; (69), 
“== () 


wobei d,' eine obere Schranke des «ten „Quadraturfehlers“ 


2 x, 
[29 30.79 1 EBETTE 79) a #262) E F \ /ı) 
x==() Xo 
ist. Setzt man noch d,’+|\r,\=d,, so folgt für |e,|, |e,|, ..., |*;| folgendes System linearer 
Ungleichungen 
IR 
Y 
le,|z&xÜMÄ zZ 0, |-Je,|+d, 
==() 
i k 
I. = xh > 19% 8214, 
%=—=0O . . . . . . . . . . . . . (71). 
k 
y 
\alSxh Mlor.|-|ex|+ di 
“0 





Die k Größen |e;|, |e.|, . . ., [ex] genügen also 2% Ungleichungen, nämlich den Ungleichungen (71) 
und den k Bedingungen |*,|=0. Man überzeugt sich leicht, daß sämtliche Voraussetzungen, 
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unter denen der Hilfssatz Il in 6 abgeleitet wurde, erfüllt sind. (Daß in (59) für k— 


+) ® 
o ein 
o gleich Null ist, ist belanglos.) Durch Auflösung des Gleiehunessvstems 





I; 
E „h >’|o,xl- En d, A Er 


(12) 
z ==) 

erhält man also % positive Zahlen Y,, für die Je, | =Y,, |s,|=<Y,, ..., |. |< Y, gilt. Die 

Determinante des Systems (72) ist von Null verschieden, solange <h kleiner ist als der 

setrag des absolut kleinsten Eigenwertes, solange also (67) erfüllt ist. Die Auflösung des 


Systems (72) geschieht praktisch ebenfalls durch Iterationen; diese konvergieren bei genügend 
kleinem #%h sehr schnell. 





Schließlich sind noch Formeln für die Berechnung der k Quadraturfehler anzugeben. 
In Analogie zu den Gl. (23) und (25) erhält man zunächst für die Ausdrücke (70) 


I 


(vu — k)u— k+1)... (uw —1)u 
DHkr1) . -. 
v) 
Hierbei ist », eine (von u abhängige) Zwischenstelle im Intervall (0....). Berücksichtigt 
man, daß AUFD — „hr fin) ist und bezeichnet man mit |f**V|, „, das Maximum des Absolut- 
betrages, das die (k®-+1)-te nach ‚x genommene Ableitung von f(x, ,y) auf der Integralkurve 


im Intervall (@#,....r,) annimmt, so darf man jedenfalls setzen 
(u- Ku k+1l)... um —1)u 
ya k+r2., (k+1) k 
} h f Dr \ (ki)! du 
"(u k)...(u  1)u 
=d’+her2.|fkri|. "I j 
Te ee a 
I: 
(w—k)...um—1)u 
d,. d. 4’ BAr2.| far) , / 
I} kr th / max \ (kLD! du 
k-1 





10. Praktische Durchführung der Rechnung und Zahlenbeispiel. Die Gl. (62), an denen 
die Konvergenz der Iterationen untersucht und die Abschätzung des Fehlers durchgeführt 
wurde, sind für die praktische Rechnung durchaus ungeeignet. Daher muß wieder der Über- 
gang zu den entsprechenden Differenzenformeln mit Hilfe der Beziehungen 


«=: -DAkted) It —:::+-V dA k=1,2,..,%&). . . (5 
zemacht werden. Man erhält so zunächst aus (62) %» Formeln für jedes k, nämlich je eine 
für 4, Mo» Ms Mor =: Mr Mo. MPureh Subtraktion je zweier aufeinander folgender erhält 
man daraus Formeln für 1, Mas Ms Mir +++ Mk Ya; unter diesen müssen nun gewisse 


in den Zahlenkoeftizienten übereinstimmen, dürfen sich nur in der Anzahl der vorkommenden 
Differenzen unterscheiden, da offenbar für die Koeffizienten nur der Abstand der zur Ver- 
wendung kommenden Differenzenzeile im Schema von dem zu berechnenden y-Wert maß- 
eebend ist. Daher lassen sich die an Stelle von (62) tretenden Rechenformeln in folgender 
Weise schreiben 





’r1) »r1) <y) ' Hi) 4 2 [? .«y) 212) E (9) 251 j? vv) y 
Yr-3 Y:-4 Äh Iı .) Ifı 719 1 ar I "0 I er 
- .)+) ») 19 | 
(+1) v+1) 4 .) EPORG +) er +.) ar . u 
AT ANT AT Fan (76) 
| £ „ | 
y- l ’-+-]) u. Pa .) .ıy) +) ep ()) | 3 f ’) | 11 R I’) | 
IR -1 YE-3 7 hf > ıfı "19 Pi" 75 "ff, + 0 "fr. +.. | 
,—] +] . ıy) | . pr | or On | . » ıy) 19 » yı 
N, UFR l I 1 ) ıfı > "fi. 4 "fi 20 "fı Tee 
Sollen #,. Ys+ ».-. 4% berechnet werden, so sind die k letzten dieser Gleichungen zu benutzen, 


wobei überall nach der A-ten Differenz abzubreehen ist. Tritt im Subtrahenden auf der 
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linken Seite der untere Index O0 auf, so kann natürlich dort der obere Iterationsindex fort- 
gelassen werden, da , fest gegeben ist. 

An und für sich wäre auch das System (63) zur praktischen Rechnung nicht ungeeignet, 
sobald y,”, f”, 9", @=1, 2, ...., k) einmal bekannt sind, es liefert aber im allge- 
meinen etwas ungünstigere Ergebnisse als (76), da zur Berechnung von y,”*" noch Funktions- 
werte vom (» — 1)-ten Schritt Verwendung finden. Man wird (63) ähnlich wie Gl. (16) 


(v) 


aber trotzdem dann vorteilhaft benutzen, wenn man nach Berechnung der y,"” bei dieser 
Iteration stehenbleiben, dabei aber noch abschätzen will, um wieviel sich die y,"”*" von den 
y,„”” unterscheiden; der Vorteil liegt darin, daß man das Differenzenschema der f,”” nicht zu 
berechnen braucht und man es nur mit sehr kleinen Zahlen zu tun hat. 

Als Beispiel möge wieder die Dt, S. S9 und hier in 5 behandelte Differentialgleichung (44) 
mit dem Anfangswert „,=0, y(#,)= 91 dienen. Die Ausgangswerte mögen für k=3 
berechnet werden; man hat also später die Möglichkeit, bei Verwendung des Extrapolations- 
verfahrens mit dem Näherungsgrad 3, bei Verwendung des Interpolationsverfahrens sogar 
mit dem Näherungsgrad 4 die Integralkurve weiter zu berechnen. 

Zunächst wird Af,, dann mit Hilfe der ersten Gl. (61) y,'” berechnet, sodann hf”, 

If,” h(f,” - f,), dann mit Hilfe der zweiten Gl. (61) 9,” und weiter der Reihe nach 
nf." HAL®, NEL y® hf" hAf®, hf”, hA®f®. Damit ist ein „vorläufiges“ Diffe- 
renzenschema aufgebaut worden. Nun finden die drei letzten hinter den Gliedern mit 4° 
abgebrochenen Gl. (76) Verwendung; sie liefern mit » = 0 Werte für y,°, 9,”, 9,°’, mit denen 
verbesserte Funktionswerte und ein verbessertes Differenzenschema gebildet werden. Um den 
Einfluß der Abrundungsfehler zurückzudrängen, ist dabei die Genauigkeit der Rechnung um 
eine Dezimalstelle erhöht worden, was auch in der folgenden Fehlerabschätzung seine Be- 
rechtigung finden wird. Nochmalige Anwendung der Gl. (76) mit » 1 liefert nun Werte 
y,Q 9," A N” » deren erste 6 Dezimalstellen sich durch weitere Iterationen nieht mehr ändern 
und, wie ein Vergleich mit Wt, S. 89 zeigt, mit den dort berechneten Werten genau überein- 
stimmen. Die Zahlenwerte sind in folgender Tafel gegeben: 











z y'® hi” IhAFPIhA?FP| RA y hi hafı’ IhA?FPIRAFFV y' 
0.00 1.000 000 0.020 000 1.000 000 0.020 000 0 1.000 000 0 
769 7695 
0.02 1,020 000. 0.019 231 4-55 1.019611 0.019 2305 554 1.019 6103 

714 6 7141 70 
0.04 1,038 846 0.018517 r 49 1.058 481 0,0185164 + 484 1,038 479 4 
— 665 — 6657 
0.06 1.057029 0,017 852 1.056 662 0.017 850 7 1.056 659 2 


Bei der Fehlerabschätzung sind folgende Werte zu benutzen: h=0,02, 2 = 053, #h= 0,006, 
E,| 0), Aus 


et y? 
w_ -q (5) Suse 40 . _ 
Di. (et 


folgt ferner f” „.x 640. Damit ergeben sich aus (74) die Quadraturfehler 


d,—5,41:10%, d,=853:10%, d,— 13,94: 10%. 


(32° 10a’ y + 13a y? — 16°) 


3 
19 11 19 
720° 720° 720° 
Eine einzige Iteration mit den Gl. (72) gibt, ausgehend von Y," = d,, Y,"”-d,, Y,"”=d,, 
bereits mit hinreichender Genauigkeit die oberen Fehlerschranken für '&,, &,, 8, 


Y,=5,45:10*, Y,=859:10*, Y,= 14,07: 10%, 


Die drei in (74) auftretenden Integrale haben der Reihe nach die Werte 


(2) 


so daß, wenn von den Abrundungsfehlern abgesehen wird, bei den y“ in der obigen Zahlen- 
tafel die siebente Dezimalstelle um höchstens 1 bis 2 Einheiten unsicher ist. 

Man könnte einen Mangel dieser Bestimmung der Ausgangswerte darin vermuten, daß, 
während die er bei Verwendung von genügend vielen Gliedern die 
Werte y(#,), y (x), beliebig genau zu approximieren gestattet, hier die Größe der Fehler 
a ee 2 die ( Genauigkeit der Annäherung, außer durch die Eigenschaften von f(x, Y), 
auch dure hAhund k festgelegt ist. In den meisten Fällen ist aber, wie im Beispiel, die auf 
dem hier geschilderten Wege erreichte Genauigkeit sogar noch größer als erforderlich, so daß 
die oft mühsame Potenzreihenentwicklung durch ein Verfahren ersetzt ist, das auch die für 
die Untersuchung der Fehlerfortpflanzung nötige Kenntnis der Fehler der Ausgangswerte 
vermittelt. | 191 
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KLEINE MITTEILUNGEN 


Zur elektrischen Festigkeit am Rande 
eines Plattenkondensators. Bei einem 
Plattenkondensator, der aus zwei unendlich 
dünnen Platten besteht (Abb. 1). die sieh nach 
einer Seite ins Unendliche erstrecken und zwi- 
schen denen eine zewisse  Potentialdifferenz 
herrscht, treten, wie Maxwell zuerst zezeigt hat. 
an den Kanten unendlich hohe Feldstärken auf. 
Kirchhoff hat nachzewiesen,. daß auch für den Fall 
endlieher Dieke der Platten die Feldstärke ın den 
Keken unberrenzt hohe Werte annimmt. Um ein 
Durchschlagen des Kondensators zu vermeiden, 
muß man danach trachten. daß die Feldstärke in 
der Nähe der Kanten nieht höher ist als im Innern 
des Kondensators. Praktisch wird dies durch Ab- 
runden der Kanten erreicht. Um das Maß der er- 
forderliehen Abrundunz zu ermitteln und um den 








\bb. 1. 


Wert der Feldstärke längs der Abrundune an- 
reben zu können. kann man nach einem Vorschlag 
von Rogzowski!) so vorgehen, dab man eine 
Niveaufläche herausgreift. z. B. die schraffierte 
Niveaufläche der Abb. 1. Diese Fläche denkt man 
sieh erstarrt und betrachtet sie als neue Konden- 
satorplatte an Stelle der ursprünglichen ebenen 
Platte. Länzes der neuen Fläche läßt sieh die Feld- 
stärke berechnen. Man hat nun die Aufgabe, durch 
Probieren eine solche Niveaufläche auszuwählen, 
länes welcher die Feldstärke keinen Wert erreicht. 
der größer ist als im Innern des Kondensators. 
Rozrowski hat diese Untersuchungen durehgzeführt 
und auf Grund seiner Ergebnisse eine geeienete 
Form für die Kondensatorplatte angegeben. 








C Ä Br 
6 ı — 
[77 
£ er tt 
fd Y fc. 
RAz3 D 


\bh. ?. 


Im folgenden soll ein anderer Weg aufgezeigt 
werden, der es ermöglicht, auf direktem Wege die- 
Jenige Form der Kondensatorplatte zu finden, auf 
welcher die Feldstärke einen konstanten Wert be- 
sitzt. Wir kommen dazu, wenn wir weitzehend 
(ebraueh machen von der Analorie,. die zwischen 
elektrischen und hydrodynamischen Feldern be- 
steht. Wir können ja beispielsweise Abb. 1 auch 
als die Strömung einer reibungslosen Flüssirkeit 
auffassen, die aus dem äußeren Raum in den Spalt 

1) W, Rogro wski r Die elektrische Festirkeit am Rande 


eines Plattenkondensators. \rehiv für Elektrotechnik. 


XII (1923) S. 1. 


zwischen den beiden Platten einströmt. Es ent- 
sprechen sich folgende Begriffe: 


Niveaulinien Stromlinien 
Kraftlinien - Potentiallinien 
Feldstärke Geschwindigkeit. 


Der unendlieh großen Feldstärke an den Kanten 
entspricht demnach unendlich große Geschwindig- 
keit der Strömung. Da aber eine unendlich große 
(seschwindiekeit nieht möglich ist, so löst sich in 
Wirklichkeit die Strömunz an der Kante ab und 
tritt in Form eines Strahles in den Raum zwischen 
den beiden Platten hinein (Abb. 2). Diese Strömung, 
die in der Hydrodynamik als Einströmung in eine 
Bordasche Mündung bekannt ist. ist von Helm- 
holtz berechnet worden. Sie hat die besondere 
Kirenschaft, daß weren des konstanten Druckes 
in der Umzebunz des freien Strahles nach der 
>ernouillischen Gleichung auch die Gescehwindig- 
keit längs des freien Strahles konstant ist. Daran 
ändert sich nichts. wenn der freie Strahl etw: 
durch eine feste Wand beerenzt würde, die genau 
der Form des Strahles folet (Abb. 2 unten). da wir 
ja eine reibungslose Flüssirkeit angenommen 
haben. 

Deuten wir nun umgekehrt dieses Strömungs- 
feld als elektrostatisches Feld. so ergibt sich, daß 
bei einem zeladenen Kondensator von der in 
Abb. 2 gezeichneten Form die Feldstärke längs der 
Perrenzung A—B (bzw. D—E) konstant ist, von 
4 nach € (bzw. D nach F) hin nimmt sie geeen 
den Wert VO ab. Legt man den Anfangspunkt eines 
Koordinatensystems in den Punkt A, und die 
negative N-Achse in die Richtung A--C, die posi- 
tive Y-Achse senkrecht dazu nach unten. so sind 
die Koordinaten der Beerenzuneskurve durch die 
folzenden Gleichungen zereben: 


IP „9, 2 
7 sin’ — In c085 


“5 pr 


Y: > (9 — sin O). 
"T 

wobei © von O0 bis z geht. Dabei bedeutet 45 den 
Abstand der beiden geraden Begrenzungsstücke 
AC—DF. In analoger Weise läßt sich der Ausfluß 
eines Strahles aus dem Spalt einer ebenen Wanil 
als eine Kondensatorform deuten, welche ebenfalls 
länes der beiden Kurven, die die Strahlgrenzen 
darstellen, konstante Feldstärke ergibt. Die Me- 
thode. nach welcher die Begrenzungskurve sowie 
das hydrodynamische (bzw. elektrostatische) Feld 
ermittelt werden kann, ist an verschieilenen Stel- 
len ?2) ausführlich geschildert, so daß es nicht nötig 
ist, darauf einzugehen. Es sollte hier nur gezeigt 
werden, daß die Methoden, die in der Hydrodyna- 
mik zur Bereehnung von Flüssigkeitsstrahlen ge- 
bräuchlieh sind, auch zur Lösung gewisser elektro- 
statischer Aufgaben heranzezoren werden können. 

, Wieselsbergzer 237 


ist Schrägüberschreiten der Fahrbahn 
gefährlich?! Bei Gerichtsverhandlungen über 
Kraftfahrzeuzunfälle wird sogar von 
Sachverständigen die Meinung geäußert, daß einen 
Verletzten die Schuld träfe, weil er vor dem heran- 
nahenden Kraftfahrzeug die Fahrbahn schräg über- 
sehritten habe. 

Das vorgeschriebene reehtwinklige Über- 
aueren des Straßendammes ist zweifellos zweck- 


0) 


?\ Vel. z.B. H. Lamb, Lehrbuch der Hydrodynamik. 
R.v. Mises, Berechnung von Ausfluß- und Überfall- 
zahlen. Z.d. V.D.I.61 (1917) S. 47. 





A er I 


. 


”.. 





Band 12, Heft 1 
Februar 1932 





Buchbesprechungen 61 





mäßige, da es die gerinzste Wahrscheinlichkeit be- 
dinet, daß der Fußgänger - überhaupt überfahren 
wird. Dareren ist die Bedrohunz dureh ein be- 
stimmtes Fahrzeug in diesem Falle nicht am klein- 
sten. Das folgt schon aus der Betrachtung des 
Grenzfalles: Ein mit Fahrzeugzgeschwindigkeit lau- 
fender Fußzänger wäre zerade bei der Bewegung 
in der Fahrbahnlängsrichtung sicher. 

In der Tat gibt es einen günstigsten 
Schräzwinkel für das Überschreiten der 
Fahrbahn. In der Abbildung stellt das größere 
Rechteck ein Kraftfahrzeug von der Breite b, das 
kleinere einen Gegenstand (Handwazen, Fuß- 
gänger usw.) von der Länge p-b dar, der sich mit 
der k-fachen Geschwindigkeit des Kraftfahrzeuges 
(k <1) unter dem Winkel « gegen die vorschrifts- 
mäßige Überquerung nach rechts über die Fahr- 
bahn bewegt. Der Kraftfahrer erblickt die Spitze 
des Gegenstandes im Abstand gb vor seinen 
linken Rädern. 

Gefährdet ist alsdann das Ende des Gegen- 
standes und, da dessen verhältnismäßige Ge- 
schwindigkeitskomponente in Richtung der Fahr- 
bahn, k-sin a, stets kleiner als 1 ist, das rechte 
Vorderrad des Kraftfahrzeuzes. Vor diesem geht 
das Ende des Hindernisses vorbei, wenn es die 
Strecke 


b p-+ 


Cosa 


zurückzelert hat: das Kraftfahrzeug mit seiner 
er uf Se WR 
] -fachen Geschwindigkeit ist dann von diesem 
” 

Punkte noch um die Sicherheitsstrecke 
fernt, und es ist 


s-:b ent- 


1 | E %, . 
s=( I e- . 
{ I P: COS u) “ 
Für das Bremsen steht also nicht beim senk- 
rechten Überqueren (@«=- 0), sondern, wie man 


dureh Nullsetzen der Ableitung nach « erkennt, für 


sina=k 


die größte Länge zur Verfügung. Sie hat den Be- 
trag 


pr] 1 k? 
Ss —— ( 
{ I 
Nr . . . ur) » 
Für die mindestens erforderliche .Sicht- 
strecke“ g:b ergibt sich unter der verein- 


fachenden Annahme » =1 (Hindernislänge gleich 


Wagenbreite): 


1-HV1—R? a 
= h :cig —. 


u a : . 
noch für /y k hat die mindest er- 


Aber > 
forderliche Sichtstreeke den zleiehen Wert wie für 


das vorschriftsmäßize Querüberschreiten, nämlich 


ö p+1 
IR 

Zwei Beispiele mözen die Verhältnisse er- 
läutern: Der Fußgänger hat die Geschwindirkeit 
6 km/h, das Kraftfahrzeuzr 30 km/h. Dann ist das 
Geschwindigkeitsverhältnis %# = 0,2 und der gün- 
stigste Schrägwinkel 11.5° Die Mindestsicht 
strecke, auf der der Führer das Hindernis vor 
seinem linken Rade erbliecken müßte, um sicher 
hinter ihm vorbeizukommen. betrüge (dessen 
Länge wieder gleich der Wazrenbreite angenom 
men) 99 Wagenbreiten zewen 10.0 beim senk 
rechten Überqueren. Aber noch bei einem Scehräg- 
winkel von 22,6° genügen diese 10,0 Warenbreiten. 


IN 














RZ) 
Abh. 1. 
Weiter sei angenommen, der Kraftwagenführer 


bremse beim Erblieken des Hindernisses nach Ab 
lauf der Schreceksekunde so weit. daß seine Dureh- 


schnittsgeschwindiekeit 20 km/h betrage. Der 
Fußgänger aber verdoppele seine Gangart auf 


12 km/h. Dann ist das Geschwindigrkeitsverhältnis 
k = 0,6, der günstieste Schrärwinkel schon 36.9° 
und der Mindestwert der Sichtstrecke 3.0 Waren 
breiten gezen 3.3 beim senkrechten Überschreiten 
der Fahrbahn. Noch bei dem großen Schrärwinkel 
von 619° genügen «liese 3.3 Warenbreiten zum 
sicheren Vorbeikommen hinter dem Fußgänger 
(wobei dessen Länge ungünstige eleich der Wagen- 
breite angesetzt ist). Hierbei ist noch nicht be- 
rücksichtigt, daß bei schrägem Überschreiten der 


Fahrbahn auch die für den Entschluß und die 
Wirksamkeit des Bremsens erforderliche Zeit 
erößer wird. K. Everline. 244 


BUCHBESPRECHUNGEN 


(Die hier angezeigten Bücher sind durch die VDI-Buchhandlung Berlin NW 7, Ingenieurhaus, zu beziehen.) 


Prof. Dr. teehn. e. h. PHILIPP FORCHHEIMER, 
wirkl. Mitglied der Akademie der Wissenschaften 


in Wien. Die Berechnung ebener und 
sekrümmter Behälterböden. 3. Aull. 


Verlag Wilhelm Ernst & Sohn, Berlin 1931. 49 S. 
mit 31 Abb. im Text. Preis geh. 4.80 M. 

Diese Schrift erscheint jetzt in der dritten Auf- 
lage und hieraus ist deutlich zu sehen, daß sie den 
praktischen Bedürfnissen entgegenkommt und 
dem Konstrukteur die gewünschten Behelfe für 
den Entwurf eiserner Wasserbehälter an die Hand 


eibt. Behandelt sind Behälterböden, und zwar an 
zwei Seiten aufliegende ebene und flach ge- 


krümmte Streifen. rechteekiee und kreisförmige 
allseits aufliegende «dünne Platten, flache Kugel- 
hauben,. Umdrehungstlächen, Halbkugeln. Die Be- 
lastune ist dabei eine gleichförmie verteilte, 
Ferner sind z.. T. referatartig besprochen (lie 
Aufhängung und Auflagerunz von Behältern, die 
Untersuchung der Wandungen, das Ausbeulen von 
Boden und Wandungen,. sowie «lie Einwirkung des 
Winddruckes. wobei neuere Ereebnisse Berück- 
sichtieung gefunden haben. Ein bis zur Gegen- 
wart ergänztes Literaturverzeichnis wird dem 
l‚eser willkommen sein. 


Ratzersdorfer. 
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Dr.-Inz. ROBERT SCHÖNHÖFER, orıd. Prof. an 
der Teehnisehen Hochschule in Braunschweig. 
Wirtschaftliche Stützuner von Trage 
eebilden. Zeiehnerische und rechnerische Ver- 
fahren zum wirtschaftlichen Vorentwerfen von 
Bauwerken (Deeken. Dächer. Hallen, Brücken usw. 
unter Berücksichtigung der von seiten der Praxis 
eeeebenen Bedingungen und Wünsche, Univ.-Ver- 
lae Robert Noske, Leipzig 1931. VII + 185 >. 
Preis 15 M. 


Das vorliegende Werk, das aus einigen Abhanıl 
lungen des Verfassers entstanden ist. beschäftigt 
sieh mit den Aufgraben, wie bei einem Tragrgebille. 
z. B. bei einer Brücke, die Stützen anzuordnen 
sind. «damit «die Kosten des Bauwerks einen 
Kleinstwert erreichen. Man braucht zur Lösung 
(dieser Fragen «die Angaben über ılie Kosten des 
Tragwerks in Abhängigkeit von seiner Stützweite 
und über die Kosten «der Stützen. bei ihrer An- 
lage an «den einzelnen Orten. Diese Zusammen- 
hänge werden graphisch dureh die Tragwerk- 
kostenlinie und die Stützenkosten- 
linie dargestellt, wobei «lie Kosten als Ordinaten 
zu den Stützweiten. bzw. über «len Stützenstellen 
im Längsschnitt durch das Bauwerk aufgetragen 
sind. Ist die Öffnung L einer Brücke mit gege- 
benen Emdstützen z. B. «dureh einen Zwischen- 
pfeiler in die Längen 1 und L-—l unterteilt, so 
bildet die Summe aus den Ordinaten der Trag- 
werkkostenlinie in I und L-—1l, sowie der Stützen- 
kostenlinie am Vfeilerort. «lie Kosten von Trag- 
werk umd Pfeiler. Für verschiedene 1] bekommt 
man eine Kurve der Kosten, «deren kleinste 
Ordinate der günstigsten Stellung (des Pfeilers ent- 
sprieht. Die Lösung wird jedoch sehr kompliziert. 
wenn «die Lage der Endstützen nicht feststeht. 
wenn eine Unterteilung in viele Öffnungen vorzu- 
nehmen ist und wenn in den einzelnen Öffnungen 
andere Tragwerkkostenlinien zu berücksichtigen 
sind. Der Verfasser führt «das Verfahren in allen 
Fällen zeiehnerisch elegant umd  übersiehtlich 
dureh umd mit den Annahmen ezeradliniger uni 
parabelförmiger Ersatz-Kostenlinien lassen sich 
die Aufgaben auch rechnerisch erledigen. 

Zur praktischen Anwendung ist zu sagen, dab 
eine Stützenkostenlinie verhältnismäßig leicht zu 
finden ist. Das Aufsuchen einer Tragwerkkosten- 
linie ist hinzeren umständlich. denn es bedeutet 
ja die Bereehnung «des Überbaues mit wechselnden 
Spannweiten. Aber bei einem großen Brückenbau 
wind man immer Varianten ausarbeiten und daher 
in der Lage sein. einige Punkte dieser Linie an- 
zugeben, Schließlich sind bei bauenden Behörden 
umd Unternehmungen auch schon Aufstellungen 
vorhanden, die verwertet werden können. Und so 
erhält man die Grundlagen. mit denen sich die 
Ergebnisse dieser interessanten Studien in die 
Praxis umsetzen lassen. Ratzersdorfer 28. 


Dr. KARL DOEHLEMANN, weil. Professor an 
der Technischen Hochschule München. Geome- 
trische Transformationen. Zweite, um- 
vearbeitete Auflage, herauszereben von Dr. Wil- 
helm OIlbrieh. Professor an der Hochschule für 
Bodenkultur in Wien. Mit 89 Fire. u. 4 Abb. i. 
Text. Göschens Lehrbücherei. 1. Gruppe. Bd. 15. 
Verlag Walter de Gruyter & Co. Berlin und Leip- 
zie 1930. 254 S. Preis 14,50 M. 


Die vorliewende 2. Auflare ist aus den beiden in 
der Sammlunz Schubert erschienenen Bänden .Geo- 
metrische Transformationen"“ von Doehlemann 
nach einiger Kürzunz hervorgerangen. Die Her- 
ausgabe besorgte W. Olbrieh, nachdem der Ver- 
fasser zestorben ist. Die Kürzungen beziehen sich 
hauptsächlich auf die Anwendungen in der dar- 
stellenden Geometrie. Trotzdem bietet das Buch 


auch in der vorliegenden Form manches dem an- 
vewandten Mathematiker. Es wird darin die pro- 
jektive Geometrie vom Standpunkt der Transfor- 
mationen aufgebaut. Es kommen natürlich im 
erster Linie lineare und quadratische 'Transforma- 
tionen in Betracht. Die gruppentheoretischen Ge- 
sichtspunkte gelangen nur in geringem Mabe zur 
Geltung. Nach Einführung der projektiven Koor- 
dinaten werden die projektiven Abbildungen auf 
der Geraden und im Strahlenbüschel, ferner in der 
Ebene und schließlich im Raume behandelt, sodann 
auch kurz quadratische Transformationen im 
ebenen und räumlichen Falle gebracht. Bemerkens- 
wert sind die zahlreichen eingefügten Konstruk- 
tionsaufgaben sowie die Beschreibung von Kon- 
struktionsapparaten. So findet man im S$ 18 die 
mechanischen Methoden zur Herstellung kolline- 
arer Figuren (Storehschnabel in verschiedener Aus- 
führung, Perspektograph von H. Ritter), im S 24 
die Beschreibung eines Relief-Modellierapparates, 
im $ 31 den Peaucellierschen und Hartschen In- 
versor sowie den von Sylvester und Kempe. Die 
Anwendungen auf die darstellende Geometrie (Pa- 
rallelprojektionen) werden besonders zusammen- 
vestellt. Auch die Anwendungesmöglichkeiten der 
projektiven Geometrie in anderen Gebieten( Kristal- 
lographie, Potentialtheorie) werden kurz ange- 
deutet. Ein reichhaltiges Literaturverzeichnis be- 
schließt das Werk. 


G,.Szeeö, 


OSWALD VEBLEN, Prof. of Mathematics, 
Prineetown University, Analysis Situs. 
American Mathematical Society. Colloguium Publi- 
eations, Vulume V, Part. II. Second. Edition. 1931. 
X + 1978. Preis 2 8. 


Ein dem Wesen nach unveränderter Nachdruck 
der im Jahre 1921 erschienenen ersten Auflage. 


G. Szegö. 241. 


Dr.-Ing. K. RAUH, Privatdozent für Getriebe- 
lehre an der Teehnischen Hochschule Aachen. 
Praktische Getriebelehre. 1. Banll. 
Verlag Julius Springer, Berlin 1931. VII + 1398. 
mit 196 Textabb. und 19 mehrfarb. Abb. auf 
S Tafeln. Preis geb. 22.75 M. 

Wenn Verfasser im vorliegenden ersten Banıl 
seiner Getriebelehre,. der sieh nur auf das Gelenk- 
viereek und jene Getriebe besehränkt. die in be- 
kannter Weise (daraus abgeleitet werden können. 
ach seinem Vorwort die „Zusammenhänge und 
Übergänge zwischen den einzelnen Getrieben“ an- 
sehaulich zeigen will. so dürfte ihm das, insbe- 
sondere auch dureh gute, z. T. farbige Abbildungen, 
durchaus gelungen sein. Auch die Sehilderung 
der „Veränderung in den Bewegungsäußerungen 
bei der Wahl anderer Getriebeabmessungen“ 
dürfte, soweit die Geschwindiekeit und rein an- 
schaulieh die Bahnen in Frage kommen, gelungen 
sein. Vermißt wird aber dabei die Betonung «ler 
Beschleunieung, die bei schnellaufenden Getrieben 
besonders wichtige ist und nieht vergessen werden 
darf, 

Wenn jedoeh der Verf. von einer „weit verbrei- 
teten Riehtung“ spricht. „für die die Getriebelehre 
eine interessante, zum Teil an bedeutende prak- 
tische Beispiele angelehnte Aufgabensammlung für 
die Mathematik und Mechanik ist “so ist das 
erstens sehr einseitige und zweitens möchte man 
allerdings wünschen, daß er sich der Mathematik 
etwas und der dieser Wissenschaft eigentümliehen 
Ausdrucksweise besonders bedient hätte. Warum 
wird (z, B.) der Krümmungsmittelpunkt bei 
Koppelkurven roh dureh den „Schnittpunkt“ dreier 
endlich benachbarter Bahnnormalen statt mit den 
bekannten. gerade beim Gelenkviereck einfachen 
Verfahren bestimmt? Verwundern muß auch die 
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„Schwinge“ von „Schwingen-* bzw. „Kurbel-“ 





Band 12, Heft 1 
Februar 1932 
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Zitierung einer Diplom-Arbeit (5. 85) als Gewährs- 
stelle für die Tatsache. dab ein Punkt der Gang 
polbahn in dem Augenblick, wo er zum Momentan- 
pol wird, eine Spitze beschreibt. Ebenso der bald 
foleende Satz: „Die Länge des Stillstandes eines 
solchen Punktes ist natürlich je nach der Form der 
aufeinander rollenden Polbahnen verschieden, 
immer aber. auch bei nur angenähertem Stillstand. 
verhältnismäßie kurz.“ Durch die ungewohnte 
Festsetzune, daß das kleinste Glied beim Gelenk- 
viereecek mit Kurbel. das gerenüberliegende mit 
Schwinge, die benachbarten mit Koppel bzw. Steg 
bezeichnet werden, treten Unklarheiten auf, insbe- 
sondere hört man bei festgestellter „Kurbel“ oder 
statt 
Koppelkurven. Weitere Eigentümlichkeiten an- 
zugeben verbietet der Raummangel. 

In den Abschnitten über Koppelkurven ist man 
erfreut. ihren eroßen Formenreichtum in vielen 
Abbildungen vor sich zu sehen (.„Brot“-, „Tropfen“-. 
„Bretzel“formen, „Kreiskurven“ u. ä.!). jedoch 
scheint dem Verfasser enteangen zu sein. daß es 
rerade über Koppelkurven mannigfache Literatur 
rjbt (Burmester, R. Müller. Allievi u. a.), welche 
besonders bei der Ermittlung von Rastgetrieben 
aus Gelenkvierecken hätte herangezogen werden 
können. So hat doch schon Alt (auf Burmester 
und R. Müller fußend) in mehreren Arbeiten, auch 


‘ 
m) 
‘ 
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auf einer Getriebetagunge darüber berichtet. 
Warum werden (diese Tatsachen nieht vom Verf. 
in die vielleicht für die Praxis „anschauliche“ 
Form gebracht? Wenn Verf. versuchen würde, 
diese Kluft zwischen den in z. T. rein mathema- 
tischen Arbeiten niedergeleeten Ideen und den 
Weren des systematischen Probierens zu über- 


brücken (s. z. B. Beyer, Technische Kinematik), so 
würde das Kapitel über Rastgetriebe wesentlich an 
Wert. Klarheit und Exaktheit gewinnen. Dab 
keine Angaben über den Verlauf der Beschleuni- 
eung bei Rastgetrieben gemacht werden, ist be- 
dauerlich. 
Wenn 
Wertvolle 


das Werk auch manches anschaulich 
enthält. so wird die Getriebeliteratur 
durch die mit Nachdruck vorgetragenen, primi- 
tiven und mit Unexaktheit gepaarten Methoden 
nieht bereichert, und auch beim Praktiker, für den 
ddas Buch ja besonders geschrieben ist, wird eher 
Verwirrung als Klarheit geschaffen. 


Kobern-Mosel. W,.Meyerzur Capellen. 247 


Rendieconti del Seminario Mate- 
matieo e Fisieo di Milano. Nol. II (1928 

- VD. 1929. X + 199 8... Vol. II (1929 — VID. 
19030. XV + 267 8. Libreria Editriee Politeeniea. 
Milano. 

Der Band II enthält nach einem Nachruf auf 
Aldo Pontremoli u. a. folgende für die Leser 


dieser Zeitschrift wichtige Beiträge: Einen knappen 
Überblick von Fubini über die Ziele der projek- 
tiven  Differentialgeometrie. Eine  zusammen- 
fassende Darstellung des Problems der Änderung 
der geographischen Breite durch Ceeehini. Einen 
Aufsatz von Levi-Civita über die maximale dyna- 
mische Beanspruchung elastischer Systeme. 
(‚anassini handelt von der mathematischen Formu- 
lierung einiger hydraulischer Probleme. Gabba 
schreibt über Schiaparellis Bedeutung für die Ge- 


sehiehte der Astronomie. PBottani behandelt die 
Mathematik der dreiphasigen Wechselströme. 
Amerio bespricht die neuere Entwicklung der 
Unterwasserakustik. Danusso schreibt über die 
Statik erdbebensicherer Konstruktionen. Band III 
enthält folgende Beiträge: Belardinelli, Endliche 


Deformationen des Gummi: Maggi, Über die Glei- 


chungen der Dynamik:  Finzi,. Statistische 
Mechanik: Cisotti. Grundlagen der Elastizitäts- 


theorie; Marcolongo, Die Dynamik Leonardo da 


Vineis; Lori. Ohr, Telephon und Grammophon: 
Volterra, Biologische Fluktuationen; Bozza, Ther- 
modynamische Potentiale; Masotti,  Himmels- 
mechanik und allzemeine  Relativitätstheorie; 
Bianchi, Über Sternparallaxen; Amerio, Über die 
Messung der Sonnenstrahlung. 


Bieberbach. 240 


sind bei 
einzeraneen 
vorbehalten): 


Ferner 
Bücher 
bleibt 


der Sehriftleitung folgende 
ausführliche Besprechung 


Technisch-wissenschaftliche 
dem Osram-Konzern. ?., 
der Hauptstelle für 
erstattung 
Springer. 


Abhandlungen aus 
Band. Herausgeg. von 
wissenschaftliche  Bericht- 
Osram-Konzerns. Verlag Julius 
1951. VI 365 S. mit 398 Abb. 


ddes 
erlin 


l.. M. BERKELEY, Addition-Subtrae- 
tion Logarithms to five deeimal places. 
White Book and Supply Co. 36 West 91 Street, 


New York City. AH + 135 8. Preis geh. 3.25 $ 
franko, 
G. GAMOW, Der Bau des Atomkerns 


unıd die Radioaktivität. Ins Deutsche 
übertragen von Ü. u. F. Houtermans. Bd. I der 
Sammlung „Neue Probleme der Physik und 


Chemie“, herauszee. v. Dr. Euren Rabinowitsch. 


S Hirzel Verlag. Leipzie 1952. X + 147 8. mit 
| Fieuren im Text u. auf einer 


. 


377 Tabellen u. 41 
lafel. Preis kart. 10 M. 

Dr.-Ing, P. WERKMEISTER, ord. Prof. a. d. 
Techn. Hochschule Dresden, Vermessungs- 
kunde. I. Teil: Stückmessung und Nivellieren. 
>. Aufl. Sammlung Göschen, 468. Bid. Verlag 
Walter de Gruyter & Co., Berlin u. Leipzig 1932. 
163 S. m. 146 Abb. Preis 1.80 M. 


Dr. KONRAD KNOPP, o. Prof. an der Univer- 
sität Tübingen. Funktionentheorie. 2. Tl.: 


Anwendung und Weiterführung der allgemeinen 
Theorie. 4. verb. Aufl. Sammlung Göschen, 703. 


Bd. 


Leipzig 1931. 136 S. mit 7 Fig. 


Dr. KONRAD KNOPP, o. 


Verlag Walter de Gruyter & Co., Berlin 
Preis 1,80 M. 


umel 


Prof. an (der Univer- 


sität Tübingen. Aufgabensammlune zur 
Funktionentheorie. 1. Tl.: Aufgaben zur 
elementaren Funktionentheorie. 2. verb. Aufl. 


Walter «le 


135 8. 


Sammlung Göschen. 877. Bd. Verlag 
Gruyter & Co., Berlin und Leipzig 1931. 
Preis 1.80 M. 


ALFRED MANES, Versicherungswesen, 
System ‚der Versicherungswirtschaft. 5. völlig 
veränderte und erweiterte Auflage in drei Bänden 
mit Einschluß der Sozialversicherung. 2. Bd.: 
Güterversicherung. B. G. Teubner, Leipzig uni 
Berlin 1981. X +B316 S. Preis geb. 20 M. 


Dr. ARTHUR HAAS, Prof. für Physik an der 
Universität Wien, Das Naturbild deı 
neuenPhvsik. 3. verm. u. verb. Aufl. Verlag 
Walter de Gruyter & Co. Berlin und Leipzig 1932. 
129 S. mit 8 Fig. im Text. Preis geb. 6 M. 


Handwörterbuch. 
von ARNOLD BERLINER und KARL 
2. Aufl. Julius Springer Verlag, Berlin 
VI1+1428 S. mit 1114 Textfiguren. Preis 
03,60 M. 


Physikalisches Herausgeg. 
SCHEEL. 
1932. 


geh. 


Allge- 
1930. 


Jahrbuch des Forschungsinstituts der 
meinen Elektrieitäts-Gesellschaft. 2. Bd. 


332 S, mit vielen Abb. Verlae Julius Springer. 
- | s 
Berlin 1931. 
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Jahrbuch für die Gewässerkunde Norddeutsch- 
lands. Besondere Mitteilungen. 7. Bd. Nr. 1. 
Über Gescehwindiekeitsformeln Von 
W.Soldan. Verlae E. S. Mittler & Sohn. Berlin 
1931. 30 8. mit 5 Zahlentafeln und 13 Abb. Preis 
eeh. 6.50 M. 


Dr.-Ine.. KARL SCHÄFER, Zur Behani 
lunze baustatischer Aufraben als 
kandwertprobleme. Verlag R. Oldenbourg. 
München umd Berlin 1931. 77 S. mit 47 Textabb. 
Preis eeh. 5M. 


Molekülstruktur. Leipziger Vorträge 1931. Her 
auseee, von Prof. Dr, P. DEBYE, Direktor (es 
’hysikalischen Instituts der Universität Leipzig. 
S. Hirzel Verlag. Leipzig 1931. VIIE + 197 S. mit 
14 Fig. u.5 Tafeln. Preis kart. 10 M. 


Dr. W. LIETZMANN, Oberstudienidirektor in 
(Göttingen. und Dr. J. JAROSCH, Direktor ıl. Bun- 
desrealschule in Wien XL. Arithmetik für 
die I. und Il. Klasse. 4. Aufl. II+139 S. 
mit 22 Fie. im Text. Aus der Sammlung ..Mathe- 
matisches Unterriehtswerk für Mittelschulen“, 
Verlag Franz Deuticke, Wien 1932. Preis 3.20 M. 


Dr. J. JAROSCH, Direktor der Bunidesrealschule 
in Wien XII. und K. PILIZOTTI, Professor an der 
Bundesrealschule in Wier XII. Darstellende 
Geometrie für die VI. und VII. Klasse ı(er 
Realeymnasien. Aus dem Mathematischen Unter- 
riehtswerk für Mittelschulen v. W. Lietzmann uni 
J. Jarosch. Verla@e Franz Deuticke, Wien 1931. 
128 S, mit 217 Abb. im Text. Preis öst. Seh. 6.60. 


Jahrbuch über die Fortschritte der Mathematik. 
Herauszeee,. von der Preuß. Akademie (der Wissen- 
schaften. Sehriftleiter Geore Feigl. Verlag 
Walter de Gruyter & Co. Berlin-Leipzig 1931. 
Bd. 54, Jahrg. 1928. Sonderheft III: Analvsis. 
386 S, Preis 42 \M. 


Bd. 55r. Jahrg, 1929. Sonderheft Il: Ana- 
Ivsis. 190 8. Preis 21.60 M. 


dl. 54. Jahre, 1928. Sonderheft IV: Geo- 
metrie 224 S. Preis 25.40 M. 


sd. 55. Jahrg. 1929, Sonderheft IV: Geo- 
metrie. 1388. Preis 15.80 M. 


Professor TULLIO LEVI-CIVITA, Caratte- 
ristiche dei sistemi differenziali e 
propagazione ondosa. Lezioni raccolte 
dal Dott. G. Lampariello. Verlag Nieola Zanichelli. 
Boloena 1931. VII-+ 108 S. 


R. COURANT, 0. Professor an der Universität 
Göttingen, Vorlesungen über Differen- 
tial- und Integralreehnune. 2. Bad. 
2, Aufl.: Funktionen mehrerer Veränderlicher, Ver- 
lag Julius Springer, Berlin 1931. VIIT+412 S. 
mit 106 Textfiguren. Preis geb. 19,60 M. 


Pietro BURGATTI, Teoria matematica 
della elastieitä. 3.Bd. der Analisi vettoriale 
generale e applieazioni. Verlag Nicola Zanichelli, 
Bologna 1931. VIII +37 S. 


Dr.-Ing. JOHANNES FISCHER, Techn. Hoch- 
schule Karlsruhe. Theorie der thermi- 
schen Meßzeräte derElektrotechnik. 
Grundlagen zu ihrer Berechnung. Ferdinand Enke 
Verlag. Stuttgart 1931. XI+ 147 S. mit 30 Abb. 
Preis zeb. 14.50 M. 


Dr. LUDWIG PRANDTL, Professor an der Uni- 
versität Göttingen und Direktor des Kaiser Wil- 
helm-Instituts für Strömungsforsehung. Abriß 
der Strömungslehre,. Verlag Friedr. Vie- 
were & Sohn Akt.-Ges., Braunschweig 1931. VI+ 
223 S. mit 221 Abb. Preis geb. 15.40 M. 


Dr. ALEXANDER OSTROWSKI, Prof, an der 
Universität Basel. Studien über den Schott- 
kyvschen Satz. Verlag B. Wepf & Üie., Basel 
1931. 112 S. Preis brosch. Schweiz. Frk. 5. 


BERNHARD HAURWITZ, Zur Theorieder 
Wellenbewegungen in Luft und Wasser. 
Aus den Veröffentlichungen des Geophvsikali- 
schen Instituts der Universität Leipzig. Heraus- 
geg,. v. L. Weickmann, 2. Serie: Spezialarbeiten 
aus dem Geophysikalischen Institut. Bd. V, Heft 1. 
106 S. mit 5 Fig. und 14 Tabellen. Leipzig 1931. 


Dr.-Ing. e. h. R. W. POHL, Professor der Physik 
an der Universität Göttingen. Einführungin 
dieMechanik und Akustik. 2. verb. Aufl. 
Verlag Julius Springer. Berlin 1931. VIII + 251. 
mit 440 Abb. Preis geb. 15.80 M. 


NACHRICHTEN 


Gesellschaft für angewandte Mathematik 
und Mechanik. 


Ende Januar ist der Jahresbericht für 1930/31 
und ein ausführlicher Bericht über «lie auf der 
Hauptversammlung in Bad Elster 1931 gehaltenen 
\orträre an alle Mitelieder versandt worden. Die 
Mitelieder werden ersucht. ihren ‚Jahresbeitrag 
für 1931/32 umeehend einzuseniden. 


Ortsgruppe Berlin. 


Am 22, Januar sprach Hr, Priv.-Doz. Dr. A. 
Weinstein-PBreslau zur Theorie der Grund- 
wässerströmung., 


Der bereits in Heft 6 des vorigen Jahrganges 
aneezeiete Vortrae von Hrn. Prof. A, Korn- 


Charlottenburg über „Probleme, die sich auf mehr- 
lagige Spulen in der Elektrotechnik beziehen“ 
wird am 19. Februar stattfinden. 


Prager Mitglieder. 


Am 19 November hielt Hr. Prof. Dr. Harry 
Schmidt-Köthen in Prag einen Vortrag über 
„Stäbe und rechteckige Platten unter dem Einfluß 
beweglicher Lasten“, 


Persönliches. 


Hr. Prof. Dr. W. Hort ist zum beamteten a. 0. 
Professor für Maschinenbau an der Technischen 
Hochschule Charlottenburg ernannt worden, 





Für die Sehriftleitung verantwortlich: Professor Dr. von Mises, Berlin NW 87, Siegmundshof 9; 


für den Anzeigenteil: Peter Valerius, Berlin NW 4. 


Druck von A. W. Ziekfeldt, Osterwieck am Harz. 
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